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Preface 


This book aims to provide a comprehensive review of the latest advances in 
applying the microfluidic platform and technologies for addressing the problems 
and challenges arising from environmental and energy areas. Nine chapters are 
included in the book, covering a wide range of topics and establishing a basic 
framework that will be helpful for readers to appreciate the important role that has 
been played by microfluidics and envision the new areas and problems that deserve 
more efforts to explore in future endeavors. 


It is the outstanding contribution from the academic researchers working in the 
fields that make this book possible. The authors of the book chapters are greatly 
acknowledged, and the sincere gratitude also goes to the publisher. 


Yong Ren 
University of Nottingham Ningbo China, 
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Chapter 1 


Droplet Microfluidic Device 
for Rapid and Efficient Metals 
Separation Using Host-Guest 
Chemistry 


Yehezkiel Steven Kurniawan, Ramachandra Rao Sathuluri 
and Keisuke Ohto 


Abstract 


Metals are pivotal elements in our daily life and industrial processes, to produce 
electronic devices, catalysts, smart materials and so on. However, they are mostly 
present as a mixture in the environment that makes their separation challenging 
over the past decade. Host-guest chemistry principle thoroughly has been used to 
design and synthesize thousands of organic receptors with high complexation abil- 
ity and selectivity to certain metal ions. On the other hand, the droplet microfluidic 
device is well-known for its unique characteristics of fluid dynamics, such as large 
specific surface area and short diffusion distance making the process robust and 
efficient. Therefore, many reports of research employ host-guest chemistry of the 
droplet microfluidic system for the effective metal separation process. This chapter 
deals with up-to-date examples of the droplet microfluidic system application for 
separation of base and alkali metals, recovery of rare-earth and precious metals and 
removal of heavy metals either from the competitive metal system or from the real 
waste solution sample through solvent extraction techniques utilizing host-guest 
chemistry principle. 


Keywords: droplet microfluidic device, host-guest chemistry, solvent extraction, 
recovery, metal separation 


1. Introduction 


Nowadays, smart and advanced technology researches have infiltrated on all 
aspects of human daily activities [1]. For example, smart devices containing metal 
elements have massively been applied in pharmaceutical and electronic fields [2]. 

In general, based on their properties and usage, metal elements can be divided into 
four groups, i.e. alkali and base metals (lithium, iron, cobalt, nickel, zinc, etc.), 
precious metals (silver, gold, palladium, platinum, rhodium, iridium, osmium and 
ruthenium), rare-earth metals (lanthanum, neodymium, europium, holmium, etc.) 
and heavy metals (cadmium, lead, mercury, chromium, etc.). Silver materials have 
been applied as an antimicrobial agent while lithium ones as well-known energy 
storage materials [3]. The demands for these smart materials have kept increasing, 
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therefore, their supply has been depleting day-by-day. Global maldistribution of 
metals, furthermore, is worsened by the metal supplies in several countries. Over 
the past several years, researchers are putting their continuous efforts on metal 
recycling and recovery process from wastewater and other secondary resources [4]. 

Metal recycling is a process to recover metals either from generated solid waste/ 
liquid metal waste or spent home appliances using several techniques for a constant 
supply to meet the demand. Many kinds of research have been conducted to evalu- 
ate the metal recycling process; however, the metal recycling is quite complicated 
due to the presence of other metals with similar physicochemical properties in vary- 
ing concentrations [1]. Additionally, heavy metals removal from the environment is 
difficult because they usually exist in trace amounts compared with other elements 
[5]. Therefore, versatile techniques which exhibit high metal selectivity and high 
recovery percentage are highly required [6]. 

Solvent extraction or liquid-liquid extraction is one of the metal separation 
techniques where the metal extraction efficiency mainly depends on the nature 
of the extraction reagents used. For example, di(2-ethylhexyl) phosphoric acid 
(D2EHPA) extracted lead(II) from chloride aqueous solution, however, the selectiv- 
ity against iron(III), zinc(II) and copper(II) ions is unsatisfied [7]. Meanwhile, a 
complete separation of lead(II) ions from iron(III), zinc(II) and copper(II) success- 
fully achieved using an amide derivative of calix[4]arene due to its high and suitable 
host-guest interaction [8]. 

On the other hand, droplet microfluidic devices are popular in biology and 
analytical chemistry areas because of their unique characteristics and advantages 
[9-15]. Numerous research papers related to rapid and efficient metal separations 
reported recently, and their numbers keep increasing with time [16-20]. This chap- 
ter serves an up-to-date review on the application of the droplet microfluidic system 
for separation and recovery of metal elements using a combination of host-guest 
chemistry principle and solvent extraction technique. 


2. Host-guest chemistry 


Host-guest chemistry is one of the chemistry fields studies various non-covalent 
interactions between host and guest molecules [21]. Host molecule defined as 
a molecular entity that forms complexes with organic or inorganic guests, or 
a chemical species that can accommodate guests within cavities of its crystal 
structure. The guest molecule is a molecule having a complementary structure 
to the binding site of the host compound, thus generating a selective and strong 
interaction between them [22]. The non-covalent interactions between host and 
guest molecules are including size effect, hydrophobic, cation-z, n-z, electrostatic, 
and hard-soft acid-base interactions as well as hydrogen bonding. Because of these 
specific interactions, a selective complexation between host and guest molecules 
could be achieved [23]. 

Thousands of host molecules are being prepared and evaluated to find out the 
best suitable host candidate for metal separations. Among them, calixarenes are 
found to be interesting host molecules and they attracted a lot of attention for 
research works due to their special inclusion phenomena [24]. Calixarenes are mac- 
rocyclic compounds that are prepared from formaldehyde and p-alkylphenol as the 
repeating units. A single step preparation of calixarenes was firstly introduced by 
Prof. C.D. Gutsche in 1978, subsequently, it has grown rapidly due to some merits, 
such as large-scale synthesis, rigid structure, easy modification and stable host 
compound [25]. It is also interesting to note that, the ring size and conformation of 
calixarenes are controlled by using template and reaction temperature. 
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(a) The structures of calix[n]arenes and (b) the conformations of calix[4]arene. 


So far, there are four main subfamilies of calixarenes based on their ring size, 
i.e. calix[4]arene, calix[5]arene, calix[6]arene and calix[8]arene as displayed in 
Figure 1 (a). By increasing the ring size of calixarenes, their conformation could 
not easily be adjusted because of weaker hydrogen bonding on the calixarene 
structure [26]. 

Calix[4]arene is the most simple and special subfamily of calixarenes because it 
serves the most stable conformations. In general, there are four available conforma- 
tions of calix[4]arene, i.e. the cone, partial cone, 1,2-alternate and 1,3-alternate 
conformations as shown in Figure 1(b). Calix[4]arenes offer ideal frameworks 
for metal separation because of their size discrimination and chelating effects that 
contributed to high selectivity and high complexation ability with certain metal 
ions [27]. Furthermore, it is possible to change its sequence of metal selectivity by 
modifying their upper and lower rims [28]. Therefore, it is reasonable to under- 
stand that a lot of research works utilized calix[4]arenes for metal extraction, metal 
adsorption, and as chemosensor of metal ions, etc. [29-31]. 


3. Droplet microfluidic device 


Microfluidic system is a miniaturized device that manipulates and controls 
fluids, usually in the range of microliters to femtoliters, in the channels with dimen- 
sions from tens to hundreds of micrometers [14]. The earliest microfluidic system 
was reported by Hodgson and Charles in 1963 while they investigated droplet 
generation when oil and water phases were introduced into the microchannel [32]. 
In 2001, Burns and Ramshaw firstly reported that the mass transfer of acetic acid 
from kerosene to aqueous phase was significantly enhanced by using droplet micro- 
fluidic device [33]. As of today, there are many types of microfluidic devices, such 
as continuous microfluidic devices, droplet microfluidic devices, paper microfluidic 
devices, and digital microfluidic devices have been fabricated and applied for many 
applications [34, 35]. Each of them serves as a unique platform for certain applica- 
tions. Continuous microfluidic devices have been widely applied to particles and 
cell separation due to density differences. Droplet microfluidic devices serve as an 
ideal platform for metal separation, nanoparticle fabrication, drug delivery, and 
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cell assay. Paper microfluidic devices have also been reported for routine chemical 
analysis and biochemical assays, while digital microfluidic devices are employed for 
immunoassays and chemical synthesis, etc. [36]. 

As the name indicates, the droplet microfluidic device generates droplets of 
two immiscible liquids such as organic and aqueous phases in the microchan- 
nel. Based on the liquids’ intersection, there are 2 types of droplet microfluidic 
devices, i.e. T-type and Y-type of droplet microfluidic devices [9, 12]. A schematic 
representation of a droplet microfluidic device is shown in Figure 2. When the 
liquid is introduced into the microchannel through an inlet which is perpendicular 
to the main microchannel of a T-type droplet microfluidic device, it was defined 
as the dispersed phase, while, the continuous phase liquid that has introduced 
through an inlet which is in line with the main microchannel of a T-type droplet 
microfluidic device. However, both phases are equal to the Y-type of droplet 
microfluidic devices [13]. 

Kashid et al. have observed the droplets generation on both T-type and Y-type 
droplet microfluidic devices in 2010 [37]. For example, on the T-type droplet 
microfluidic device, first, the dispersed phase is entered the main microchannel 
and starting to block the flow of the continuous phase by exerting the pressure. The 
generated pressure squeezes the neck of the dispersed phase and the droplet of the 
dispersed phase is formed. Since both phases are immiscible, each droplet has a role 
as an individual reactor with a short diffusion distance. Furthermore, each droplet 
serves a large specific surface area. The combination of a short diffusion distance 
and a large specific surface area contribute to boosting the mass transfer process of 
chemicals from one phase to another phase [38]. 

Even though solvent extraction offers a selective and efficient technique for 
metal separation, the extraction kinetics rate in some processes is slow and time- 
consuming [38]. For comparison, Figure 3 shows the obtained parameters and 
experimental data for silver(I) extraction using a tetramethylketonic derivative of 
calix[4]arene in batch-wise and microfluidic systems. In the batch-wise system, the 
extraction was carried out in a small glass vial with a diameter of 4 cm and height 
of each phase is about 1 cm. The extraction using a continuous microfluidic device 
was performed in a microchannel with 20 mm of length and 200 um in depth and 
width while the extraction in a droplet microfluidic device was carried out ina 
microchannel with 73 mm of length and 200 um of depth and widths. As shown in 
Figure 3, the diffusion distance of chemicals using batch-wise (20,000 um) was a 
hundred times farther than both the microfluidic devices (200 pm). Furthermore, 
the total specific surface area of the droplet microfluidic device employed was ~550 
times larger than that of the continuous microfluidic system while ~55,000 times 
larger than that of the batch-wise system. Because of that, it is reasonable to know 
that the required extraction time to reach an equilibrium state is 259,200 s (72 h), 
15.0 sand 1.00 s for the batch-wise system, continuous microfluidic system, and 
droplet microfluidic system, respectively. This result demonstrates the remarkable 
advantage of using the droplet microfluidic device for metal extraction and separa- 
tion process [39]. 
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(a) T-type and (b) Y-type of droplet microfluidic devices. 
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Figure 3. 
Comparison of batch-wise, continuous and droplet microfluidic systems for silver(I) extraction using a 
tetramethylketonic derivative of calix[4]arene [38]. 


4. Metal separations using the droplet microfluidic device 


Due to the significant enhancement of the extraction process, droplet microfluidic 
devices have been employed for metal ion separations through a solvent extraction 
technique [18-20]. The chemical structures of several extraction reagents employed 
in metal separations are shown in Figure 4. In general, metal ion concentrations were 
measured using either atomic absorption spectroscopy (AAS) or induced coupled 
plasma-atomic emission spectroscopy (ICP-AES) or induced coupled plasma-mass 
spectrometry (ICP-MS). The measured metal ion concentrations were further used to 
calculate extraction percentages (96E), distribution ratio (D) and separation factor (f) 
as mentioned in Eqs. (1)-(3) whereas [metal ion] aqueous,o is initial metal ion concentration 
in the aqueous phase before solvent extraction process, [metal ion] aqueoust is the metal 
ion concentration in the aqueous phase after t extraction time, and [metal ion] organics is 
the metal ion concentration in the organic phase after t extraction time [26]. The metal 
extraction and separation processes are summarized and listed in Table 1. 


[metal ion] aqueous,o - metal ion] aqueous, 


%Extraction = Dosral ion] ine € 
aquoeus,0 

[metal ion] organics (2) 

p - D metal ion A (3) 


D metal ion B 


4.1 Alkali and base metals 


Alkali and base metals are non-precious metals which are abundantly found in 
nature [40]. In this part, the example of utilization of droplet microfluidic system 
for extraction and separation of copper(II), iron(II), indium (II), calcium (II), 
cesium and lithium ions were discussed. 

Evaluation of copper(II) extraction using Y-type microreactor was carried 
out by Yang et al. The commercially available extractant, AD-100 which con- 
sists of 70-75% (w/w) of 2-hydroxy-5-nonylbenzaldehyde oxime as the active 
reagent in 260# solvent oil as the organic diluent was used. The concentrations 
of AD-100 used influenced the copper(II) extraction percentage. Higher AD-100 
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Figure 4. 
Host compounds and extraction reagents used for metal separations in this chapter. 


concentration gave higher extraction percentage while a higher concentration of 
copper(II) gave lower extraction percentage. The extraction percentages were 
also increased at longer extraction time, lower Reynolds number, higher pH of the 
aqueous phase, and higher temperature. Using the T-type of droplet microfluidic 
device, 96.0% copper(II) ions were extracted in 32.0 s extraction time [41]. On 
the other hand, Jiang et al. employed a Y-type microfluidic device to achieve an 
effective separation of copper(II) over iron(II) and zinc(II) with a mixture 

of 2-hydroxyl-5-nonyl acetophenoneoxime and 5-nonylsalicyladoxime in 1:1 
volume ratio, called as DZ988N, as the extraction reagent. Using the microflu- 
idic device, as much as 80.6% of copper(II) extracted at 0.78 s extraction time. 
While as much as 98.9% of copper(II) was stripped at 3.92 s by using 1.84 M of 
sulfuric acid solution in 2:1 organic to aqueous phase ratio. In a comparison of the 
microfluidic device to the batch-wise system, the separation factors of copper(II) 
over iron(III) and zinc(II) were significantly enhanced from 86.2 to 466 and 123 
to 1089, respectively [42]. 

Iron(III) extraction kinetics using D2EHPA was investigated in a Y-type micro- 
fluidic device. The experimental results showed that the extraction of iron(III) 
was controlled by mixed reaction-diffusion resistance regime on a second-order 
extraction rate equation. Moreover, since the required extraction times for iron (III) 
and cobalt(II) are different, iron(III) could be separated from the cobalt(II) in 
the metal mixture solution, which may contribute to the real hydrometallurgical 
process [43]. 

Indium(III) separation from a metal mixture solution has been successfully 
achieved using a T-type of droplet microfluidic device. The metal mixture contains 
23.0 mM of indium(III), 35.8 mM of iron(III), 44.2 mM of zinc(II), 153.8 mM 
of aluminum(III) and 141.5 mM of magnesium (II) ions [44]. It was found that 
indium(III) was extracted in 98.896, while the other metal ions were hardly 
extracted (<5.00%). The separation process yielding that the separation factors of 
indium(III) over iron(II), zinc(II), aluminum (III) and magnesium (II) ions were 
18,720, 266,373, 41,955 and 11,809, respectively, demonstrating that a high selective 
indium (III) separation was achieved within a minute. 
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Jahromi et al. separated calcium(II) isotopes, ^5calcium(II) and *°calcium(II) 
ions using dicyclohexano-18-crown-6 in a Y-type microfluidic device. Using 
a combination of the dicyclohexano-18-crown-6 and microfluidic device, 

“°calcium (II) was selectively extracted in 98.2% within 5.70 s. By using a continu- 
ous cycle extraction process, the amount of ^'calcium(II) becomes lower in the 
aqueous phase, therefore, the concentration of ^5calcium(II) ion was enriched 
[45]. Another experiment on calcium (II) extraction was carried out by using 
ionic liquids as the extraction reagent, named as 1-ethyl-3-methylimidazolium 
bis(trifluoromethylsulfonyl)imide (EMIM NTf2) and 1-butyl-3-methylimid- 
azolium bis(trifluoromethylsulfonyl)imide (BMIM NTf2) in a Y-type of droplet 
microfluidic device. The pH value of the aqueous phase did not affect the extrac- 
tion percentage because the main extraction mechanism is not an ion-exchange 
mechanism. Even though some parameters were optimized, the highest extraction 
efficiency was only 52.096, which was still unsatisfied [46]. 

On the other hand, cesium(I) extraction in a Y-type of droplet microfluidic 
system was studied by Tamagawa et al. Stable droplet generation was observed 
when the total flow rate of both phases is less than 100 mL h™. By using D2EHPA 
as the extraction reagent, the required times to reach the quantitative amount of 
extraction percentage (100%) were shortened 25 times from 1000 to 40.0 s using 
batch-wise and droplet microfluidic system, respectively [47]. 

Lithium metals are widely used as the main part of lithium-ion batteries and air- 
craft alloys, the continuous supply of lithium may become critical over time. Urban 
waste, especially from lithium-ion batteries and seawater, may serve as the potential 
secondary resource in the future. However, lithium-ion batteries contain lithium, 
cobalt(II) and manganese(II) ions, while seawater contains lithium and other metal 
ions, such as sodium, potassium, magnesium, and calcium ions in different concen- 
trations, the lithium recovery from these resources poses to be a tough work [48]. 

Muto and co-workers investigated the extraction process of lithium using 
D2EHPA in a droplet microfluidic device. They found that a glass substrate gave 
higher extraction efficiency as well as the mass transfer coefficient due to stronger 
convection rate and larger specific surface area of the organic droplets. When 2.0 M 
D2EHPA was used as the extraction reagent for 1.0 mM of lithium ions in the aque- 
ous phase, around 0.3 mM of lithium ions were extracted in 5.00 s as the extraction 
time. Furthermore, by the addition of 0.18 M tributyl phosphate (TBP), the amount 
of the extracted lithium can be improved to 0.5 mM at 5.00 s extraction time. Even 
though the extraction efficiency was only around 50%, the droplet microfluidic 
device serves as a potential platform for quick metal extraction compared with the 
conventional extraction process [49]. 

Hirayama et al. evaluated the separation of cobalt(II) over lithium from a model 
solution of lithium-ion battery waste. The mixture of 0.20 M of D2EHPA and 
0.01 M of NaOH, called as D2EHPA-Na in cyclohexane was used as the organic 
phase while the model solution containing 1.0 mM of cobalt(II) and 1.0 mM of 
lithium ions was used as the aqueous phase. From the model solution, the extraction 
reagent extracted around 50 ppm of cobalt(II) and 1.5 ppm of lithium ions giving 
3.54 as the separation factor for cobalt(II) over lithium ions. Furthermore, using the 
droplet microfluidic device, 9896 purity of cobalt(II) solution was obtained within 
4.50 s extraction time [50]. 

In our previous study, we investigated the possibility to extract lithium ions with 
monoacetic acid derivative of calix[4]arene using T-type of droplet microfluidic 
system [51]. At first, the droplet microfluidic device was fabricated from Tempax 
glass and the microchannel was produced utilizing a diamond-coated needle through 
a micromilling process. The microchannel dimensions were 73.0 mm in length and 
200 um in width and depth as confirmed by laser reflection microscope. The top 
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and bottom parts of the T-type droplet microfluidic device were permanently bonds 
through a thermal fusion technique [52]. To prevent any disturbance of leaked 
sodium ions from the glass substrate, the microchannel was coated with SFE-X008 
coating agent (1% of m-bis(trifluoromethyl) benzene in 94% of the ethyl perfluo- 
robutyl ether). After coating, the wettability of the glass substrate was significantly 
changed from 20 to 120°, demonstrating a successful coating of the microchannel. 
From the preliminary investigation on lithium extraction using the monoacetic 
acid derivative of calix[4]arene (shown in Figure 4 as 1 Ac) on the individual metal 
extraction, it was found that 100% of lithium ions were quantitatively extracted just 
within 2.00 s as the extraction time. Moreover, the complete stripping of lithium ions 
from the organic phase was easily achieved using 0.10 M HCl as the stripping agent. 
Since the monoacetic acid derivative of calix[4]arene completely recovered 
(extraction and stripping efficiencies were 10096) the lithium ions, the competitive 
metal extraction was carried out using a mixture of metal ions as the aqueous phase. 
In a mixture of all alkali and alkaline earth metals (lithium, sodium, potassium, 
rubidium, cesium, magnesium, calcium, strontium, and barium ions) in 10.0 mM 
concentration each, the 5.00 mM of the monoacetic acid derivative of calix[4] 
arene extracted 4.99 mM of lithium ions showing a quantitative (99.8%) extraction 
percentage of lithium ions. By using 1.00 M of HCl, as much as 99.596 of lithium 
ions were successfully stripped from the ligand, which is remarkable. The results 
demonstrated that the lithium separation from seawater which mainly contains 
sodium, potassium, magnesium and calcium ions was successfully achieved using 
the monoacetic acid derivative of calix[4]arene in the droplet microfluidic system. 
Seawater was sampled from Minatohama, Karatsu, Japan and filtered using 
Millipore membrane to prevent any blockage on the microchannel during the 
extraction process. The lithium concentration in the seawater was 0.03 mM while 
the sodium, potassium, magnesium and calcium ions concentration was 495, 
10.1, 474 and 12.8 mM, respectively. Using a droplet microfluidic system, it was 
found that 31.696 of lithium (0.01 mM) and 1.0096 of sodium ions (4.95 mM) were 
extracted from seawater, however, other metal ions were not extracted at all and 
remained in the aqueous phase. In this case, sodium ions were extracted due to the 
huge amount of existing sodium (~19,000 times higher) compared to lithium ions. 
To achieve complete extraction of lithium ions from seawater, 20.0 mM of the 
monoacetic acid derivative of calix[4]arene in chloroform was used as the extraction 
reagent. As expected, 100% of lithium ions were successfully extracted. However, 
since the higher concentration of the extraction reagent was used, as much as 2.5696 of 
sodium ions were also extracted. Both metal ions were completely stripped using HCl 
to obtain an aqueous phase containing 0.03 mM of lithium and 12.2 mM of sodium 
ions. To obtain a pure lithium solution, another type of calix[4]arene, named as triace- 
tic acid derivative of calix[4]arene was used because it exhibits high sodium selectivity 
over lithium ions at the acidic condition. When the 10.0 mM of triacetic acid derivative 
of calix[4]arene (shown in Figure 4 as 3Ac) was employed, 100% of the remained 
sodium ions were completely extracted leaving pure lithium ions in the aqueous 
phase. Similar to monoacetic acid derivative, triacetic acid derivative of calix[4]arene 
was also recycled. The scheme of recovery of lithium ions from seawater is shown in 
Figure 5(a). From this experiment, a rapid and efficient lithium-ion recovery from 
the seawater was successfully achieved in the three-step process using calix[4]arene 
derivatives as the extraction reagents in a droplet microfluidic device [51]. 


4.2 Precious metals 


Precious metals, i.e. silver, gold, palladium, platinum, rhodium, iridium, osmium 
and ruthenium metals are indispensable metals because of their high economic value 
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Figure 5. 
(a) A process of lithium separation from seawater [51] and (b) precious metals separation from a real metal 
waste [52] using the droplet-based microfluidic system. 


and their application for smart materials [53]. Bizzo et al. reported that silver(I), 
palladium(II) and platinum(IV) ions were found in the electronic waste serving as 
the secondary resource for precious metal ions [40]. Stepwise recovery of precious 
metals from electronic waste, therefore, must be established as an untapped resource 
for future precious metals supply. Even though the recovery process is complicated 
because of the presence of other metal ions with similar physicochemical properties, 
a stepwise recovery of precious metals from real metal waste using neutral type of 
calix[4]arene derivatives in a T-type of droplet microfluidic device was established 
[52]. Silver(1), palladium(II) and platinum(IV) ions extraction were carried out 
employing tetramethylketonic, tetraquinoyl and tetrapyridyl derivatives of calix[4] 
arene, respectively, as the extraction reagents. The chemical structure of tetrameth- 
ylketonic, tetraquinoyl and tetrapyridyl derivatives of calix [4]arene were shown in 
Figure 4 as MKTOC, QuTOC and PyTOC, respectively. 

From the individual metal extraction experiment, it was found that as much 
as 99.9% of silver(I) was extracted at 2.00 s of extraction time from 0.10 mM of 
AgNO; solution in 0.10 M HNO; media using 5.0 mM of the tetramethylketonic 
derivative of calix[4]arene in a droplet microfluidic device. However, it required 
72h to achieve a similar extraction percentage of silver(I) using a batch-wise system 
showing that the droplet microfluidic device enhanced the required extraction 
time 130,000 times faster to reach the equilibrium stage of silver(I) ion extraction 
compared with the conventional batch-wise system. On the other hand, as much as 
96.2% of palladium (II) ions was extracted at 4.00 s extraction time from 0.10 mM 
of Pd(NOj), solution in 0.10 M HNO; media using a tetraquinolyl derivative of 
calix [4]arene, and as much as 93.196 of platinum(IV) ions were extracted at 4.00 s 
extraction time from 0.10 mM of H;PtCl; solution in 0.10 M HCI media using a tet- 
rapyridyl derivative of calix[4]arene. The stripping of silver(I), palladium(II) and 
platinum(IV) ions were obtained in medium percentages using 2.00 M ammonium 
thiocyanate in water, 1.00 M thiourea in 1.00 M HCl and 1.00 M thiourea in water, 
respectively, from the metal-laden organic phase within 2.00 s only. 
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The evaluated real metal waste, mainly consisted of silver(I), palladium (II), 
platinum(IV), iron (IID), copper(II), and nickel (II) ions at 1.32, 0.22, 0.03, 0.01, 
0.82, and 0.01 M, respectively. To prevent any metal precipitations, the real 
metal waste was diluted 1000 times using 1.0 M HNO3. Since the concentration 
of silver(I) was higher than palladium (II) and platinum(IV) ions, the silver(I) 
ions were extracted first and followed by palladium(II) and platinum(IV) ion 
extraction. The scheme of stepwise recovery of those precious metal ions from real 
metal waste is shown in Figure 5(b). All extractions and stripping processes using 
a droplet microfluidic device were conducted at 4.00 s extraction time. Silver(1) 
ion extraction was carried out in two-cycle process due to high existed concentra- 
tion. On the first cycle, 80.696 of silver(I) ions (1.064 mM) were extracted using 
20.0 mM of the tetramethylketonic derivative of calix[4]arene. A complete extrac- 
tion of silver(I) ions were achieved in the second cycle. Even though only 54.896 of 
silver(I) ions were stripped from the first cycle using 2.0 M ammonium thiocya- 
nate, as much as 94.4% of silver(I) ions were recovered on the second cycle. 

A complete palladium(II) ions recovery was achieved within a single cycle 
extraction process using 5.00 mM of the tetraquinolyl derivative of calix[4]arene, 
while the stripping percentages of palladium ions were 44.0 and 73.296 on the first 
and second cycle of stripping process using 1.00 M thiourea in 1.00 M HCl. On 
the other hand, 274 and 37.296 of platinum(IV) ions were extracted in the first 
and second cycles of extraction process using a tetrapyridyl derivative of calix[4] 
arene. Meanwhile, low stripping percentages of platinum(IV) ions, 12.1 and 19.596 
of the first and second cycles of the stripping process, respectively, were obtained 
using 1.00 M thiourea in the water as stripping reagent. Even though platinum(IV) 
recovery still needs to be optimized, a stepwise and complete recovery of silver(I) 
and palladium (II) ions from real waste was established using a droplet microfluidic 
system at 4.00 s extraction time [52]. 


4.3 Rare-earth metals 


Rare-earth metal separations leave a serious challenge due to astounding proper- 
ties as well as their extensive applications in electronic and advanced materials [1]. 
Rare-earth metals including lanthanide and actinide groups which are classified 
as 4f and 5f elements in the Periodic Table. Their physicochemical properties are 
quite similar to each other, thus challenging many researchers to find the optimum 
separation process for hydrometallurgy field [30, 54]. In this part, lanthanum (III), 
praseodymium(III), neodymium (III), samarium (III) and uranium(VI) extraction 
and separation processes are described and discussed. 

Lanthanum (III) extraction from a chloride solution was studied in a Y-type micro- 
fluidic device by Yin et al. By using phosphonic acid (2-ethylhexyl)-mono 2-ethyl- 
hexyl) ester (P507) as an extractant, high extraction percentage of lanthanum(III) 
ions was achieved within only 0.37 s [55]. Furthermore, the separation of 
lanthanum(III) and europium (III) ions using a T-type of droplet microfluidic device 
has been carried out. The separation factor of lanthanum (III) over europium (III) is in 
the range of 1.50 to 3.50 depends on the flow rate of the liquids [56]. 

Separation of praseodymium (III) and cerium(III) in a Y-type microfluidic 
device employing P507 as the extraction reagent was also investigated. At 12.0 s 
extraction time, as much as 75.0% of praseodymium (III) and 42.9% of cerium (III) 
ions were extracted into the organic phase yielding 3.25 as the separation factor. 
Even though the separation was not satisfied enough, the process could be recycled 
further to obtain pure cerium(III) ions as the remained aqueous phase [57]. 

On the other hand, neodymium(III) separation over praseodymium (III) was 
studied by He et al. using a Y-type droplet microfluidic device [58]. One and a 
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half molar of P507 in sulfonated kerosene was used as the organic phase while 
praseodymium (III) and neodymium (III) in HCl and lactic acid media were used as 
the aqueous phase. At 12.0 s extraction time, it was found that P507 extracted 88.0% 
of neodymium (III) and 75.0% of praseodymium (III) yielding 2.23 as the separation 
factor of neodymium(III) over praseodymium(III). The similar separation factor 
was obtained using microfluidic at 12.0 s extraction time while the batch-wise sys- 
tem required 560 s extraction times. These findings initiated a rapid enrichment of 
neodymium from the wastewater sample within less than 600 s. From the wastewa- 
ter containing 30-90 ppm of neodymium (III) ions, the neodymium (III) ions were 
enriched 200—450 times higher which is remarkable [58]. 

He et al. also investigated the extraction process of samarium (III) using P507 
as the extraction reagent in a Y-type of droplet microfluidic device. Within 45.0 s, 
as much as 75.196 of samarium(III) can be extracted compared with the batch-wise 
system that requires 600 s to reach a similar extraction percentage [59]. Meanwhile, 
using a T-type of droplet microfluidic device, samarium (III) was extracted in 82.5% 
within 14.8 s, shorter than Y-type device. This phenomenon was caused by shorter 
diffusion distance and a larger specific surface area of T-type than Y-type of droplet 
microfluidic device [60]. 

Uranium (VI) extraction process was studied in a microfluidic device using 
several extraction reagents, such as phosphorus reagent, ionic liquids, and so on 
[61-64]. Among them, tributyl phosphate (TBP) exhibited as the best extraction 
reagent for uranium(VI) extraction process due to stable fluid dynamics, as well as, 
a possibility for large scale process through numbering-up the device. Darekar and 
coworkers investigated the extraction of uranium (VI) in 1.0 M HNO; using 30% 
v/v TBP in dodecane as the organic phase. It was found that quantitative (10096) 
extraction percentage was achieved at 1.50 s extraction time, however, the separa- 
tion of uranium(VI) over other metal ions has not been reported yet [65]. 


4.4 Heavy metals 


Heavy metals, especially lead (II) ions have been used in several industrial 
processes over the past several years and they get released into environment causing 
serious effects on environment and human health. A serious control monitoring of 
heavy metal ion concentration at industrial wastewater effluents should be compre- 
hensively managed. So far, a lot of heavy metal removal techniques have been evalu- 
ated for wastewater treatment, however, they were unsatisfied because of heavy 
metals existed in trace amounts together with huge concentrations of base metal 
ions [5, 6]. In this part, two examples of a selective removal process of lead(II) 
removal using two derivatives of calix[4]arene, i.e. tetraacetic acid derivative (4Ac) 
and tetrakis(diethylamide) derivative of calix[4]arenes (EATOC) will be discussed. 

In a conventional batch-wise system, the equilibrium state of lead(II) extraction 
using 4Ac reached after 24 h, however, it was shortened to 6.00 s using a droplet 
microfluidic device. The significant enhancement was observed by shorter diffu- 
sion distance and larger specific surface area aforementioned above. The lead (II) 
extraction percentages were increased from 73, 89 and 10096 at 8.00 s for pH 2.00, 
2.25 and 2.50, respectively, however, the separation of lead(II) and other metal ions 
was unsatisfied. As shown in Figure 6(a), the extraction percentages of iron (III), 
copper(II), and zinc(II) were 60, 6.9 and 7596, respectively. This low selectivity 
of lead(II) extraction over Fe(III) was caused by the ion exchange mechanism in 
which the acetic acid groups of calix[4]arene were deprotonated and Fe(III) ions 
were easily extracted due to strong electrostatic interactions [66]. 

Another neutral calix [4]arene host molecule, EATOC was also evaluated to 
obtain a selective lead(II) separation from other metal ions mixture solution. 
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Figure 6. 
Lead (II) separation from other metal ions from a competitive metal mixture solution using (a) 4Ac [66] and 
(b) EATOC [8] as the extraction reagents in a droplet microfluidic systems. 


In this case, a neutral host molecule was used to achieve high metal selectiv- 
ity. In 1.00 M HNO; media, one molecule tetrakis(diethylamide) derivative 

of calix[4]arene formed a 1:1 complex with one lead(II) ion revealed from 

the Job’s plot experiment [8]. Using FTIR and ‘H-NMR titration studies, it 

was elucidated that the lead (II) ion was chelated by EATOC through phenoxy 
oxygen atoms and carbonyl oxygen moieties. Since the lead(II)-EATOC com- 
plex was stable enough, a selective lead(II) extraction (100%) was achieved 
from a competitive metal system containing copper(II), zinc(II), nickel(II), 
iron(III) and cobalt(II) ions at 0.10 mM concentration of each metal ion [8]. 
A complete lead(II) removal was also successfully performed from a simulated 
wastewater effluent containing 0.10 mM of lead (II), 45.0 mM of iron(III), 
while copper(II), zinc(II), nickel(II) and cobalt(II) was in 30.0 mM concentra- 
tion, respectively as shown in Figure 6(b). Furthermore, the lead(II) ions was 
easily stripped in 98% using distilled water as the stripping reagent, which is 

a convenient and low-price process. These finding demonstrated that droplet 
microfluidic system is a promising device for a complete and efficient lead (II) 
removal from environmental samples [8]. 


5. Conclusions 


Several successful reports on optimized separation and recovery of alkali and 
base metals, precious metals, rare-earth metals, and heavy metals using the droplet 
microfluidic system have been published. By employing the host-guest chemistry 
in a droplet microfluidic system offers a rapid and efficient metal separation from 
binary metal mixtures or real samples either from nature or industrial activities. 

In order to realize the true potential of droplet microfluidic devices for renewable 
metallurgical applications, the improved design of microfluidics addressing larger 
volumes with continuous operation along with effective utilization of host-guest 
chemistry are greatly needed to become a step forward in building a novel metallur- 
gical technique for selective metal separations from a real sample, which otherwise 
not technically feasible to extraction using the conventional batch-wise system. 
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with Off-the-Shelf Devices and 
Other Novel Designs 
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Abstract 


For the past three decades since the dawn of the concept of scaling down fluidic 
processes and systems, various microfabrication techniques have been significantly 
explored and developed. Glass and elastomers, as first-generation microfluidic 
device materials, are still widely used, while other alternatives have emerged. We 
have seen a rise in novel and innovative device designs and fabrication in droplet 
microfluidics which enable easier and more cost-effective approach of generating 
droplets, such as the use of commercially available “off-the-shelf” components, 
plug-and-play modular devices, and 3D-printed droplet generators. This chapter 
aims to review some of these facile approaches for droplet generation and discuss 
the versatility and functionability of these designs for possible commercial formu- 
lations. Discussions of rationales for and challenges of the development of these 
conceptual devices are included. 


Keywords: off-the-shelf, 3D printing, droplet generator, modular, parallelization 


1. Introduction 


Microfluidics defines the science and engineering of a small-scale fluid system. 
It is a study of design, fabrication, and operation of a system conducting fluids in 
microscopic channels of widths or diameters ranging from 10 to 500 micrometers 
(um). Gafian-Calvo first introduced the use of micron-scale capillaries later estab- 
lished by Thorsen et al. who demonstrated the use of a simple T-junction microflu- 
idic device to control the flow of immiscible liquids, forming monodisperse droplets 
[1, 2]. The manipulation of discrete fluid packets in microdroplets provides benefits 
in large reduction of reagent volume, size of sample, and the equipment [3]. This 
opened a completely new wave of interest in microfluidic systems in a broad range 
of fields such as biochemical, engineering, and pharmaceutical benefits from the 
advantages of small and precise. 

Droplet-based microfluidics, one subcategory of microfluidics, focuses on the 
discrete volume creation with the use of immiscible fluids and allows the handling of 
fluids under confined spatial and temporal control. Applications of droplet micro- 
fluidics are largely distinguished into two aspects: droplet reactors for biochemical 
reactions and analysis [4] and droplet-templated material synthesis. The proverbial 
use of lab-on-a-chip (LOC) in various biochemical analyses is driven by the time and 
cost-efficiency of performing automated high-throughput analysis with small volume 
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of reagents, especially in the areas of genomic study, point-of-care (POC) diagnostic, 
drug discovery, etc. To date, only a handful of selected microchips have been com- 
mercialized [5]. Compared to conventional bulk methods, microfluidic techniques 
provide a power platform that enables the creation of highly controllable emulsion 
droplets in which the size, shape, and composition of the droplet can be manipulated. 
This allows the synthesis of functional droplets such as microcapsules for drug 
delivery and cell encapsulation, microparticles, and various types of Janus particles. 

With the rapid development in microfluidic systems and its applications, differ- 
ent selections of techniques and materials are increasingly accessible for fabrication 
of microsystems (especially for those who are just starting out in the field). Some of 
these fabrication techniques require high-priced tools and machineries with skilled 
workers to produce intricate designs of microchannels and microchambers in which 
droplets are generated and manipulated, while simple assemblies of microdevices 
with commercially available components have also been reported. To produce 
emulsions, which are droplets in immiscible phase, the selection of materials with 
compatible surface wettability is important. 

The huge potential of droplet-based microfluidics has stimulated rapid develop- 
ment of flexible platforms to generate highly monodispersed droplets with diverse 
structures. The recent years see a soaring number of researches and publications on 
life sciences and material synthesis applications based on microfluidic approach. 
Likewise, there has also been a rise in novel and innovative device designs and 
fabrication in droplet microfluidics which supposedly enable easier and more cost- 
effective approach of generating droplets, such as the use of commercially available 
“off-the-shelf” components, plug-and-play modular devices, and 3D-printed 
droplet generators. While this signifies a positive outlook on the development of 
microfluidics, the ultimate reception of the technology depends primarily on the 
end users. Different applications where microfluidic approach is employed have dif- 
ferent prerequisites for commercialization, for instance, for material synthesis such 
as catalyst preparation or drug formulations, high throughput is expected to gener- 
ate sufficient materials, whereas for chemical/biological assays, throughput might 
not bea focal criterion, instead development of chip-to-world interfaces for stable 
connection to analytical units would be crucial. Up until now, the successfully com- 
mercialized microfluidic products are mainly for POC diagnostics, genotyping and 
sequencing, and other biomedical-related applications [5]. 


2. Summary of existing droplet-based microfluidic technology 


As the standard material interface used in the field of microelectronics which 
pioneered microfluidics, glass and silicon are naturally the first-generation materi- 
als for microfluidic devices [6]. Originally, glass and silicon substrates are processed 
via standard photolithography which mainly comprises three steps: (i) lithography 
or pattern transfer onto a substrate, (ii) etching of microchannels, and (iii) lastly 
bonding to create an enclosed structure [6]. The process of photolithography has 
been well documented. Devices fabricated via lithographic methods generally have 
quasi-two-dimensional (2D) planar flow [7]. Umbanhowar et al. and Utada et al. 
both introduced three-dimensional (3D) co-flow and hydrodynamic flow-focusing 
devices with tapered capillaries to produce monodisperse single and multiple 
emulsions, respectively [8, 9]. Both dispersed and continuous phase fluids meet in 
parallel streams in the co-flow geometry, while for the flow-focusing device, the 
dispersed and continuous phases are introduced in opposite directions, whereby 
the flow is eventually focused through an orifice. By combining both geometries, 
monodisperse multi-emulsions could be generated via single-step emulsification 
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(see Figure 1). 3D axisymmetric flow-focusing devices circumvent the wetting 

of channel walls by the dispersed phase as the effect of wall channel is of less 
concern and droplet formation frequency is higher [10-12]. Microgrooves can also 
be formed on glass substrates via deep reactive ion etching (DRIE). Nisisako and 
Torii have produced a planar synthesis silica glass chip with 256 parallel generators 
formed by DRIE which could produce emulsions up to a rate of 320 mL/h [13, 14]. 
However, micromachining of glass is very expensive, and it gets even more difficult 
and costly to fabricate elaborated and complex designs on a glass chip. 

Soft lithography, an extension from conventional photolithography, allows the 
lithographic master to be used for rapid prototyping of elastomeric material such as 
polydimethylsiloxane (PDMS). PDMS has been and is still a popular option because 
it is cheap and elastic and it is easy to work with high flexibility [15]. Its intrinsic 
properties such as high permeability to gases and optical transparency have even 
greatly extended its use in biomedical researches such as cell culturing and tissue 
analysis [16]. PDMS devices (see Figure 2) manufactured via soft lithography are 
inherently 2D; however, fabrication of 3D microfluidic devices is now possible 
with advances in 3D printing, whereby 3D PDMS microchannel can be casted from 
3D-printed molds [17]. 

With the vast selections available, researchers often establish their choices for the 
development of microfluidic systems based on the applications, material compat- 
ibility, ease of fabrication, as well as the start-up cost incurred; many methods/tech- 
niques remain underdeveloped or ceased to be developed presumably due to high 
start-up cost. While glass and PDMS have been extensively used in pioneering many 
microfluidic researches and applications, the known disadvantages associated with 
the fabrication and use of these materials have often been deliberately overlooked by 
the microfluidic community. As PDMS is inherently hydrophobic, small hydropho- 
bic nonpolar molecules tend to be absorbed into the polymer, which could signifi- 
cantly disrupt the microenvironment at which biological studies are carried out [15]. 
Economy-wise, it would be very costly to work with glass due to laborious proce- 
dures and difficulty in reproducing, while expansion of PDMS for industrial work 
will be less viable. Apart from glass and PDMS, plastic has also been a popular option 
and one that is believed to play a major role in translating microfluidic research into 
commercialized technologies [19]. There has been an increasing interest in develop- 
ing materials such as polymethymethacrylate (PMMA), polycarbonate, and cyclic 
olefin (co)polymers through microfabrication techniques such as hot embossing, 
injection molding, micromilling, and stereolithography [20, 21]. Plastics are more 
rigid and hence more collapse resistant than PDMS, and they generally have good 
compatibility for biological applications [19]. 

Hot embossing and injection molding are both replicating techniques akin to 
soft lithography, with the use of thermoplastic materials. With hot embossing, a 
thermoplastic film is patterned against a master, and a cast is formed as the molds 
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Figure 1. 
Glass capillary microfluidic device for double emulsion formation through coaxial jet [9]. 
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Figure2. 

Double emulsions produced from a PDMS device replicated via soft lithography. Part of the device is subjected 
to plasma treatment to render it hydrophilic. (a) Optical image of a double emulsion forming at the second 
junction of device. (b) Double emulsions collected with high uniformity. Reproduced from Ref. [18] with 
permission from the Royal Society of Chemistry. 


are heated and compressed [22]. Injection molding involves the heating and subse- 
quent injection of a thermoplastic melt into mold cavity, after which the cast will be 
removed from the mold once cooled down. While the resolution of molds and casts 
for both methods is subject to the methods used to make the mold inserts, especially 
complex three-dimensional structures, such replicating techniques are promising 
for rapid prototyping of microdevices. Micromilling, which has been a less popular 
alternative compared to mold replicating, is a conventional machining approach 
of incising features on workpieces. Advances in the milling technology now 
allow 3D computer-aided design (CAD) models to be prototyped on automation 
through computer numerical control (CNC) mills with microscale resolution [19]. 
Micromilling can also be used to fabricate molds to be used in replicating devices. 
While these microfabrication techniques are more established now than two 
decades ago, ongoing researches are directed at reducing the cost of the prototyping 
techniques with cheaper materials. Currently, hot embossing and injection molding 
are limited to thermoplastics only, while the main concern for micromilling would 
be whether the quality of micromilled devices would suffice for practical use. 


3. Some novel and innovative designs of droplet generators 
3.1 Off-the-shelf devices 


With rapid development in droplet emulsification, there has been much consid- 
eration in the use of “off-the-shelf” devices. While the term “off-the-shelf micro- 
device" might be ambiguous and may refer to already commercialized microchips, 
here it is defined as droplet generators assembled using commercially available 
components. Some of these components include medical dispensing needles, 
laboratory tubing, and tiny plastic fittings (such as ferrules and flanges) that might 
have been used as universal fluidic fittings or connectors to, for instance, laboratory 
analytical units. These easily assembled 3D devices are seen as cheap yet practical 
alternatives to conventional droplet generators requiring laborious effort to be 
fabricated. Simply said there could be no need for lithographical pattern transfers, 
microforging of fine glass capillaries, or any other hefty microfabrication units 
for the synthesis of these devices. Such robust and disposable devices from quick 
assembly of commercially available components are also useful and cost-saving 
when hazardous experiments are conducted [23]. Unsurprisingly, there have been 
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some publications that explore and incorporate the use of such droplet generators 
for production of functional materials such as microcapsules for drug delivery or 
gas sorption [24, 25]. Apart from being highly cost-effective and easy-to-make, 
some have reported that their off-the-shelf devices could also be easily disas- 
sembled for cleaning and there is also higher flexibility in droplet generation as the 
device can be reconfigured for desired formation [26]. 

Figure 3 shows the different devices assembled with components from commer- 
cial microfluidic sources, IDEX Health & Science. In Figure 3a, a 3D microfluidic 
nozzle is formed using a micrometering valve assembly and other fluidic fittings 
in an effort to focus and inject particulate samples into a downstream fused silica 
microfluidic device [27]. A silica capillary with a cone shape tip is used as the 
delivering nozzle and is sealed to the microferrule with a 5-min epoxy. Figure 3e 
and f shows a capillary microfluidic device with two polyetheretherketone (PEEK) 
chromatography tees, while Figure 3d shows a micro-cross droplet generator with 
high-pressure PEEK adapters [26, 28]. Microferrules and female nuts are used to 
secure the capillaries in place. According to Benson et al., the device can produce 
both oil-in-water (O/W) and water-in-oil (W/O) emulsions by just flushing with 
ethanol before the formation of the other, and the device is also reusable and could 
be disassembled for cleaning. The tips of two glass capillaries are flamed to create 
two orifices: one for the injection of the inner fluid and the other for the focusing 
of the flow. With the flamed-tip design, there is no need to modify the wettability 
of the capillary, allowing changeable formation of O/W and W/O. The distance 
between the orifices of the two aligned capillaries can also be adjusted, and the flow 
regime in the device changes even when the flow rates are maintained. Wu et al. 
compared the synthesis of droplets in the micro-cross droplet generator to that from 
a planar PDMS chip having a similar geometry. The former has a higher droplet 
formation frequency due to faster 3D fluid thread collapse than in the 2D flow. The 
micro-cross droplet generator is also able to operate under high flow pressure of up 
to 400 psi [28]. 

Steinbacher et al. demonstrated an assembly of a simplified mesofluidic 
device entirely from off-the-shelf components, small-diameter tubing, barbed 
tubing adapters, and needles, into droplet-forming devices and particle concen- 
trator [23]. The readily available components can be arranged into a T-junction 
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Figure 3. 

Microfluidic devices made from off-the-shelf components from commercial sources. (a) A 3D microfluidic 
nozzle for focusing and injecting of particulate samples. (b-d) Off-the-shelf components and the assembly of 
a micro-cross droplet generator. (e and f) Image and schematic of a reusable modular glass capillary device 
assembled from commercial materials. Scale bar = 1 cm. Reproduced from Refs. [26-28] with permission from 
the Royal Society of Chemistry. 
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device or a flow-focusing configuration, and the formation of monodisperse 
droplets, microcapsules, and Janus particles has been reported to be successful, 
with the results comparable to those synthesized with more complicated devices. 
Li et al. also reported the assembly of yet another entirely off-the-shelf micro- 
device with flexible designs [29]. The device is mainly composed of dispensing 
needles of different sizes (60-1550 um) and assembled with mini tee- and cross- 
links as shown in Figure 4. The stainless-steel dispensing needles are arranged 
and aligned to form microchannels for the flow of the immiscible phases. The 
assembling of the device is made using plastic tubing as well as UV curable 

glue. Monodisperse single and multiple droplets have been successfully formed 
using the needle-based microdevice by reconfiguring the device. The size of the 
droplets can be controlled by varying the needle used as well as the flow rates of 
each phase. 

The past few decades have seen a myriad of microfluidic device designs, yet the 
innovative assemblies of off-the-shelf devices have shown that the basic modalities 
of a droplet generator could be substituted by commercially available components, 
tee-link for a T-junction, cross-link for a cross-channel, dispensing needles and 
glass capillaries as the inlet and flow channels, and other structures that could also 
be replaced by commercially available parts, which greatly reduce the fabrication 
cost to as low as 3USD per device. Assemblies of such devices have also shown that 
without sealing the droplet generators onto rigid substrates, the devices are recy- 
clable and can be reconfigured for versatile droplet generation. The development 
of such microdevices has seemingly negated the need for complicated microfab- 
rication techniques such as those listed previously. However, the need for manual 
assembly as well as the difficulty in aligning glass capillaries remains an issue in 
the fabrication of these devices as most of them still adopt the use of glass capillar- 
ies, especially when nontransparent fittings are used, it poses an obstacle for clear 
observation of droplet formation. Alas, there has not been many further researches 
reported following successful demonstrations of the use of these devices, and it is 
questionable if these devices would stand any chance to be commercialized. To the 
best of the author's knowledge, there has yet to be any work reported on the paral- 
lelization of off-the-shelf devices. In the meantime, such devices remain a short- 
term solution for preliminary droplet generation studies and a conceptual proving 
tool for new droplet-based applications. Nonetheless they would be useful for 
applications in detection or quick analyses, given their versatility to be connected or 
integrated to other existing units. 


(8) middie phase 
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Figure 4. 

(a) Schematic illustration of generation of double emulsions from needle-based coaxial microfluidic device. 
(b) Clear depiction of production of W/O/W emulsions. (c) The assembly of the needle-based device with 
a coin for scale. (d) Extended configuration of needle-based microdevice for complex emulsion production. 
Reproduced from Ref. [29] with permission from the Royal Society of Chemistry. 
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3.2 3D-printed devices 


Additive manufacturing (3D printing) has been vastly adopted in fabrication of 
microfluidic devices with the advance of the technology in generating structures at 
increasingly high precision. 3D printing is a layer-by-layer manufacturing technology 
that now allows an extensive range of materials to be printed, such as various types 
of plastics and even glass [20, 30]. As is it largely automated and assembly-free, rapid 
prototyping can be achieved for microdevices. There are different types of 3D print- 
ing techniques, but they are primarily categorized into (i) extrusion-based 3D print- 
ing, (ii) stereolithography, and (iii) multijet (or polyjet) modeling printing [17]. 3D 
printing of microfluidic devices with elaborated and complex designs is now possible 
with the advances in the resolution and accessibility of 3D printers. In this section, 
an emphasis is given to 3D-printed modular devices, to introduce and highlight the 
versatility of droplet generators with “plug-and-play” functionality, which could bea 
good solution for chip-to-world integration and development. 

Vijayan and Hashimoto developed 3D-printed fittings to form axisymmetric flow- 
focusing droplet generators with other readily available materials such as needles and 
elastic tubes [31]. While monolithic 3D-printed device is now possible with the rapid 
development in 3D printing technology, the authors have opted to print the device in 
parts and assemble with other commercially available parts, as they claim that this gives 
more customizability to the devices configuration, allowing for the generation of higher 
order and complex emulsions. Akin to the off-the-shelf needle device by Li et al., com- 
mercially available dispensing needles of varying inner diameters (60-250 um) are used 
to produce droplets of a larger range of sizes, while for fully 3D-printed microchannels, 
the minimum channel features achievable are a few hundred micrometers [32]. Multiple 
emulsions and Janus particles can also be formed by connecting the device in series or 
incorporating a 3D-printed branched Y-channel such as that shown in Figure 5g. 

For stereolithography and polyjet 3D printing, post-processing is needed to 
remove support materials or nonpolymerized resins from narrow channels; and it 
is even more challenging and time-consuming to do so for monolithically printed 
devices [31, 32]. Modularly printed microfluidic devices, however, allow easy 
removal of uncured materials and surface treatment onto different modules to 
facilitate generation of complex emulsions, and external active components may 
also be incorporated to the device as individual modules. Such modular plug- 
and-play units would also be useful in chip-to-world interface to existing systems 
requiring microfluidic solutions. Many have produced comprehensive reviews on 
the advances of 3D-printed devices, including works on plug-and-play modular 
devices, and the works included in the following discussion merely demonstrate the 
potential of 3D printing in microfluidics and hence are far from exhaustive. 

Ji et al. presented a modular multimaterial 3D-printed device that allows both 
passive and active (pneumatic control) generations of various emulsions [11]. Three 
modules have been printed, a function module, a T-junction module, and a co-flow 
module, of which at least two modules are needed for the assembly of a functional 
device. There are three switchable modes for the function modules, namely, single- 
inlet and dual-inlet modules, depending on the type of emulsions to be generated or 
a pneumatic control unit (PCU), for active control of droplet generation by control- 
ling the deformation of the flexible channel. A multimaterial polyjet 3D printer 
is used to print the PCU module which consists of both elastic (orange) and rigid 
(gray) parts, as indicated in Figure 6a, while the T-junction and co-flow modules 
are printed with a stereolithography printer using a clear resin. Surface treatment is 
carried out on individual modules to render the surface hydrophobic/hydrophilic. 
Figure 6d-f shows passive generation of single and multiple emulsions in the differ- 
ently configured device. Likewise, single and double emulsions can be generated by 
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(a) 3D-printed modules and fittings for the fabrication of microfluidic axisymmetric flow-focusing device: 
needle (N), 3D-printed holder (H), 3D-printed Y-fitting (Y), and tube (T). Schematics and optical images 
of (b) a single flow-focusing device, (c) serially connected flow-focusing devices to produce double emulsions, 
and a (d) flow-focusing device with a Y-shaped parallel laminar flow to produce compartmented particles. 
Scale bar = 10 mm. (e) a schematic illustration of serially connected flow-focusing devices for producing 
double emulsions. (f) W/O/W double emulsions with varying numbers of inner droplets. Scale bar = 600 um. 
(g) Schematic illustration of production of bi-compartmented Janus particles with a Y-channel as inlet. (h) 
Optical micrographs of produced compartmented particles at different flow conditions. Scale bar = 300 um 
[31]. Published by the Royal Society of Chemistry. 


replacing the single—/dual-inlet module with a PCU, and under different applied air 
pressure, the flexible channel undergoes deformation through pneumatic excitation, 
and regimes of unstable and stable droplet generations are identified. The modules 
are connected through snap-fit joints printed to the modules with an O-ring, and no 
leakage is detected during the experiment even at a pressure of 4 bar. 

Song et al. also fabricated a modular microfluidic system that allows both passive 
and active generation and manipulation of complex emulsion droplets [33]. Various 
modules such as inlet/outlet, I-junction, Y-junction, co-flow droplet generator, mixer, 
and a module containing electrode are printed with an i3DP printer with clear trans- 
parent UV curable polymer; however, the visibility of the modules is slightly reduced 
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Figure 6. 

go ed modular device. (a) Schematic of a single droplet-generating pneumatic device consisting of a PCU 
and a T-junction. (b and c) Schematic and optical image of a pneumatic device for generation of double or 
multicore emulsions. (d and e) Schematic and optical image of modular device with dual inlet. (f) Frames 
from high-speed recording showing complex emulsions produced via dual-inlet device [11]. 
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Figure 7. 

Schematic of 3D-printed modular microfluidic system. (a) 3D-printed modules of different functions and 
structures. (b) A full assembly of the modular device. (c) The design of the interlocking male and female 
connectors. (d) Male and female connectors printed on each module [33]. 


due to surface roughness. The modules have interlocking male and female connectors 
that secure the assembly, and an inherently printed O-ring groove around the female 
connectors circumvents the leakage around the joint. An assembly of all the modules 
into an operative device is depicted in Figure 7. When the generated droplets reached 
the electrode module, they experience an electric field in a direction perpendicular 
to the flow direction. The shape of the droplet is distorted as an effect to the electro- 
hydrodynamic flow patterns within and outside the droplet. Again, by reconfiguring 
the device, versatile range of droplets has been formed such as alternating droplets 
and single- and dual-core double emulsions. Under the influence of the electric field, 
the droplets are subjected to electrocoalescence to form Janus structure. 

Unlike other methods, 3D printing is an additive technique, and 3D parts are 
formed layer by layer without requiring molds. A variety of 3D printing methods 
has been developed; likewise, a wider range of choices of polymer-based print- 
ing materials are available, although some are proprietary materials limited to 
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certain printing techniques. At present, 3D-printed devices are limited by the 
minimum achievable channel dimensions, subject to the printing resolution. The 
printing time of a device is also subject to the complexity of its design structure. 
Nonetheless, given the current development trajectory in additive manufacturing, 
there is high anticipation on 3D printing to resolve challenges with current existing 
microfabrication techniques such as high start-up cost, low channel resolution, and 
manual assembly of devices. The demonstration of plug-and-play modular devices 
is believed to be able to partly resolve the long-standing chip-to-world issues. 


4. Droplet-based microfluidics: now and then 


Droplet microfluidics is an enabling platform that is now widely recognized and 
adopted in many research disciplines. The past few decades have seen a rapid evolu- 
tion in microfluidic-based technologies, especially in contributing to life science- 
related researches and gradually into material sciences. Much could be benefitted 
and have been realized from the subscale world of micros and nanos, attributed to 
rapid analyses and fast reactions, high precision, and good control of droplet size 
and compositions. Progressively, droplet-based platforms have been used for the 
synthesis of functional materials such as drug delivery vehicles and nanomaterials 
[34]. Conversely, the low production rate of droplets at 0.1-10 mL/h hinders the 
extensive use of the technology for high-volume production [35]. The direct solu- 
tion to addressing low throughput of droplets is to increase the number of droplet 
generators to scale up production of droplets within a single chip. 


4.1 Microfluidic emulsification: toward parallelization 


Such effort to scale up the production of chip-based emulsions has been 
reported as early as 2008 by Nisisako and Torii who demonstrated parallelization 
of up to 128 and 256 droplet generators on planar glass chips to achieve throughput 
of 128.0 and 320.0 mL/h of Janus droplets and single emulsions, respectively [14]. 
Microgrooves are etched onto synthesis silica glass substrate by DRIE to form differ- 
ent droplet generator geometries such as cross-junction, Y-shape co-flow, and even 
three-consecutive cross-junctions, joint in a circular manner for the production of 
various emulsions [13, 14]. The etched chip is sealed with another glass substrate, 
and the microfluidic chip is mounted on a stainless-steel holder connected to the 
inlets of different liquids. Some of the designs are illustrated in Figure 8, while 
Figure 9 shows the arrangement of commercialized modules sharing a single set of 
infusion pumps, for generation of Janus microparticles of a few hundred kilograms 
per month [36]. 

Other parallelization devices have also been reported, especially for mass pro- 
duction of single-phase emulsions. The channels in the parallelization devices are 
commonly distributed in a ladderlike geometry, in a treelike branched geometry, or 
in a row, while different layers are needed for the introduction of all liquid phases 
through through holes to be distributed into the microchannels [35, 36]. Apart from 
DRIE, parallel channels can also be micromachined onto PMMA substrate, repli- 
cated onto PDMS via lithography, or 3D printed [37, 38]. A PDMS millipede device 
is also reported, with more than 500 nozzles arranged in 2 rows, and droplets are 
formed through the nozzles and break under static instability [39]. Jeong et al. dem- 
onstrated the mass production of W/O emulsions at kilo-scale using 1000 parallel 
flow-focusing drop generators arranged in a 20 x 50 dense array on a PDMS device 
which is currently the highest degree of parallelization ever achieved (see Figure 10) 
[40]. Currently only Nisisako’s team has demonstrated mass-production of Janus 
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Figure8. 

Designs of glass microfluidic chips with parallelized channels in circular arrangements. (a) 144 cross-flowing 
droplet generators on two concentric circles, (b) 72 droplet generators on three concentric circles, and (c) 32 
emulsion generators in four concentric circles producing triple emulsions. Reproduced from Ref. [13] with 
permission from the Royal Society of Chemistry. 


Figure 9. 
Parallelized devices are arranged in series to further scale up the production of Janus particles. Reproduced with 
permission [36]. 


particles for commercial applications; however, there is no other information on the 
status of the parallelized device, whether the technology has been commercialized. 


4.2 From “proof of concept” to commercialization 


Inkjet printing is a classic example of commercialized microfluidic application to 
date, having been developed in the 1950s and largely commercialized and manufac- 
tured in the 1970s for digital image printing, through impelling droplets onto papers. 
The inkjet technology is subsequently adopted as a patterning tool for material 
deposition onto a wider range of substrates, for applications such as direct printing of 
electronic and optical devices [41]. More recently, studies have been conducted on the 
synthesis of microparticles via emulsion solvent evaporation following inkjet print- 
ing of droplets on substrates [42]. Following the debut of inkjet printing, notably 
only a limited number of microfluidic applications successfully made their way into 
the market in the following years, mainly in the areas of microfluidic genotyping and 
POC diagnostics. Chin et al. and Volpatti et al. summarized the commercialization of 
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Figure 10. 

Sehetnati diagram of a PDMS parallelization device. (A) the channels are arranged in a 20 x 50 array. (B) 
The device has a ladderlike geometry, and solution is distributed evenly across the device. (C) Cross-section 
view of the 3D microchannel. (D) Schematic of the distribution of solution from the delivery channel to the 
through holes and then the droplet generator; droplets are formed and then exit through the outlet channel. 
Reproduced from Ref. [40] with permission from the Royal Society of Chemistry. 


microfluidic devices of which most are LOC-based POC diagnostic devices, report- 
edly commercialized by a few companies [5, 43]. These POC devices developed by 
individual companies cover a vast range of purposes, ranging from flu, cardiovas- 
cular diseases, malaria and E. coli detection to blood chemistry analysis. Devices for 
cancer detection and DNA and RNA signatures are still under development. These 
devices come in different forms such as disposable cards or cartridges containing 
preloaded test solutions and capillary-driven test strips, and the entire analysis could 
be performed on-chip or off-chip with handheld analyzers or benchtop instruments. 
More conventional POC devices are expected to adapt LOC concepts and functions in 
order to improve diagnosis and to meet greater clinical needs. 

Alas, the overwhelming accomplishments in microfluidic research and develop- 
ment are not reflected in the slow uptake of microfluidic technology in the market. 
The industrial realization of microfluidic innovations is most often held back by the 
inflexibility of microfluidic components for integration, economic non-viability, and 
the lack of standardization. A complete microfluidic system usually entails many 
different fluidic components, and a single manufacturing process is not possible at 
current stage without the intervention and support from manual labor. Automated 
and assembly-free prototyping of simple devices is perceived to be a goal that is not 
attainable in the near future. Most importantly, with many deemed “promising” inno- 
vations, the challenge is to ensure that such innovations can be adapted and integrated 
to other existing appliances, technologies, or systems. Chin et al., Volpatti et al., 
and Mohammed et al. have all echoed the same concern that the focus should be on 
practical, marketable devices that are readily integrable, instead of unending develop- 
ment of single microfluidic component with niche function. Such is the demerit of 
introducing microfluidic-formulated solutions; they are too unique to be adapted; 
hence, the issue of chip-to-world interface should also be addressed and resolved even 
from the initial design phase [44]. And unless a minimum volume demand is met 
and a reasonable margin is secured, it is both uneconomic and unsustainable for mass 
production of microdevices, a risk that no stakeholders would bear [45]. 

The translation of the technology seems more challenging as it appears that academic 
researchers are more interested in publishing their inventions rather than developing 
them into commercially viable products. Many successful inventions remain a “proof 
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of concept” and ceased to be developed. In fact, many of the microfluidic systems are 
limited to “one design one application” which are applicable only to niche research 
areas, rendering them less appealing for further development. The voluminous lateral 
development of microfluidic researches and lack of standardization on methods and 
materials instead become the stumbling block for a standard microfluidic solution to 
reach mainstream uses, hence slowing the commercial uptake of the technology [5, 44]. 


5. Conclusion 


Microfluidic-based researches are growing in an unprecedented pace with more 
research disciplines adopting the technology, techniques, and tools that are now much 
established. At present, novel microfluidic platforms offering diverse functionalities 
and solutions are surfacing each day, and there seems to be no limit to what the 
microscale technology could offer, be it in clinical technologies or the manufacturing 
sector. While researches are flourishing in the academic perspective, the ambiguity 
surrounding the commercialization of microfluidic applications did not go unno- 
ticed. The interest of much of the microfluidic community is in lateral development 
of microfluidic solutions with only a small amount ongoing researches targeting at 
resolving technical challenges for innovation-to-market translation: the technologies 
and economy for mass productions. Despite the continuous effort and enthusiasm for 
more microfluidic solutions to be realized industrially, the path to commercialization 
and industrialization for more microfluidic solutions depends potently on the joint 
effort and commitment from the researching teams and their industrial counterparts. 
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Abstract 


The advance of droplet-based microfluidics has enabled compartmentalization 
and controlled manipulation of monodispersed emulsions with high yield and 
incorporation efficiency. It has become a highly exotic platform in synthesizing 
functional material due to the presence of two immiscible liquids and the inter- 
face between them. With its intrinsic feature in high degree of product control, 
advanced emulsion-based synthesis of functional material is constituted as a 
template for effective water remediation and resource recovery. This chapter aims to 
provide an overview of recent advances in microfluidic technology for environmen- 
tal remediation. More specifically, the facility of microemulsion-based functional 
materials for water remediation is reviewed. Moreover, the removal and recovery 
of pollutants, such as heavy metal, dye, pharmaceuticals, etc., from aquatic envi- 
ronment by the applications of adsorption on functional micro/nanomaterials are 
unfolded with respect to its potential for wastewater purification. 


Keywords: microfluidics, remediation, functional sorbents, recovery, wastewater 


1. Introduction 


Over the past few years, tremendous growth in the manufacturing and wide- 
spread application of synthetic chemical compounds in industrial sectors has led 
to growing number of emerging contaminants in environmental matrices (air, 
wastewater, water, sediment, and soil), which poses a challenge for regulatory 
agencies. Apart from that, the environment has gradually suffered as a result of 
exposure to the emerging pollutants especially in water and wastewater. Various 
types of physical, chemical, and biochemical methodologies have been reported on 
effectively removing pollutants. These methods include membrane separation, bio- 
logical degradation, advanced oxidation process, and adsorption as well [1]. Recent 
developments have indicated that the adsorption method is favorable due to its wide 
availability, lower cost, as well as its recyclability. As in terms of cost effectiveness, 
simplicity of construction, and easy adaptation of operating conditions, particle 
adsorption materials become the key for realizing various adsorption applications in 
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environment remediation using fixed bed reactors, absorption columns, fluidized 
beds, and cyclone separators [2]. However, controlled shapes, size, and compart- 
ments are some of limitations for most conventional methods. Unlike conventional 
approaches, microfluidics is comparably conducive as it has improved and extended 
the possibilities to synthesize highly controlled size of microparticles with excellent 
adsorption capability and reusability [2]. 

Microfluidics is a multidisciplinary field stretch across engineering, physics, 
chemistry, microtechnology, and biotechnology. Microfluidic devices generally have 
two of the three geometric length scales in the order of microns. The micrometer 
length scale defines the most obvious but extremely important character of micro- 
fluidic devices especially their small size, which allows small sample volumes, low 
cost, and fast analysis, but with high resolution and sensitivity [3]. With the length 
scale associated with microfluidic devices, the flow within them tends to be laminar. 
Moreover, one characteristic of microscale miniaturization is the large surface 
area to volume ratio. Thus, this favorable aspect plays a major role in control and 
manipulation of fluid flow in microfluidics. Recent advances and innovations could 
make microfluidics technology ubiquitous and create microfluidic devices that are 
more functional, efficient, and cost effective than conventional techniques. 

The development of microfluidic systems that allow for the formation of 
microdroplets inside microfluidic devices has gained greatly attention over the 
past 20 years. The rise in interest is due to the utilization of microfluidic devices 
in a broad range of biological and biomedical application areas including disease 
diagnosis, cell treatment, drug screening, single-cell analysis, and drug delivery. 
Liquid droplets dispersed in a second immiscible fluid are useful, particularly when 
the sizes and the size distribution of droplet can be prescribed on a few hundred 
nanometers to a few millimeters [4]. Microfluidic emulsification approach offers an 
alternate and versatile route to produce emulsions that are highly monodispersed 
and have high formation frequencies in multiphase fluid systems [5]. Additionally, 
the geometrical attributes and flow characteristics within these microfluidics 
system constitutes flexibility in producing complex structured emulsions, such 
as double-emulsions and multi-emulsions [6, 7]. Highly monodispersed single, 
double, or multi-emulsions can be used as a template to prepare micro- and 
nanoparticles with various structures and morphologies [7], as shown in Figure 1. 
As seen in Figure 1, the presence of different particle shapes, compartments, and 
microstructures is formed based on the flow and geometrical attributes in micro- 
fluidics. Consequently, the superior properties of droplet-based microfluidic device 
have become extremely promising and attractive platform that enables the produc- 
tion of functionalized monodisperse microparticles. 

In order to fabricate micro- or nanoparticles, the analysis of droplet formation 
is imperative to understand the device operation and its process control to meet 
different application purposes. Additionally, other’s droplet manipulation, such as 
fusion, fission, mixing, and sorting with high precision and flexibility constitutes 
essential issue, at which extensive investigations have been directed. Based on the 
sources of the driving force involved, there are five approaches to droplet manipula- 
tion: hydrodynamic stress, electro-hydrodynamics, thermos-capillary, magnetism, 
and acoustics [8]. Hydrodynamic stress is a simple and effective approach to accom- 
plish droplet manipulation relating to the geometrical characteristics of microchan- 
nel [9]. In this mode of manipulation, various methods have been employed in the 
formation of droplets in microfluidics device, including co-flowing mechanisms, 
flow-focusing mechanisms as well as cross-flowing mechanisms [10]. These 
mechanisms enabled the formation of dispersions with highly attractive features, 
particularly the control over droplet and particle size distribution. 
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Figure. 
Classification of various structures of microparticles produced by droplet-based microfluidics. 


Over the years, microfluidic devices have been developed to synthesize particles 
for water remediation. For instance, Zhao et al. [11] synthesize graphene oxide 
microspheres using microfluidics technology for the removal of perfluorooctane 
sulfonate. Dong et al. studied the anionic dye adsorption using chitosan micropar- 
ticles [2]. In addition, the performance of microcapsules for the CO, adsorption 
and permeability was investigated by Stolaroff et al. [12]. Hitherto, there are still 
little attempts that have been considered to use microfluidic platform as a selection 
scheme for large-scale industrial wastewater treatment. Toward practical and high 
capacity, microfluidic platforms generally suffer from the high cost and limited 
capacity for high throughput production of microfluidic synthesized particles. 
Nevertheless, it is possible to operate the production in parallel in order to realize 
continuous processes [13]. Moreover, the fundamental microfluidic research is 
still highly demanded to bridge the gap between the functional material synthesis 
and industrial perspective on exploring the possibilities and potential benefits of 
microfluidic processes [2, 13]. 

In this chapter, we will discuss the current trend of employing microfluidic 
technologies for environmental remediation, specifically for wastewater treatment 
and water remediation. Apart from that, we will also provide a general overview 
of the facility of microfluidics emulsification for the fabrication of various mic- 
roparticles and nanoparticles as functional adsorbents. The sorption capacity 
and performance of the functional materials will be also evaluated in this review. 
Eventually, this chapter provides an impression of what are the consolidated fields 
of microfluidic formulation in functional material synthesis that will look like in 
about a decade from now. 
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2. Microfluidic technology for water remediation and resource recovery 


Due to excessive discharges of harmful wastes and by-products to the environ- 
ment, water contamination is the most prominent in which numerous organic and 
inorganic pollutants are found in the fresh water resources such as ponds, rivers, 
and underground water. Traditionally, there are various methods of water analysis 
including atomic adsorption, chemical analysis, chromatography, colorimetry, and 
spectrometry. Although these techniques have high sensitivity and accuracy, limita- 
tions such as expensive instrumentation, time consuming and requires manual 
operation in sampling process, causing the researchers to shift their interest toward 
microfluidics technology, which has a great potential to replace the traditional water 
analysis techniques [14]. 


2.1 Droplet microfluidic system 


Droplet-based microfluidic technology involves the formation and manipula- 
tion of discrete droplets inside microdevices [15]. Droplet microfluidics is being 
widely used in different applications, such as chemical reactions, therapeutic agent 
delivery, imaging, biomolecule synthesis, diagnostic chips, and drug delivery [16]. 
The major concerns when creating droplet microfluidic device are the type of 
microfluidic chip fabrication materials used and the fluids used for droplet genera- 
tion. Poly(dimethyl)siloxane (PDMS) is commonly used as fabrication material 
of microfluidic devices because PDMS is a relatively low cost and easily moldable 
elastomer. Nevertheless, PDMS has low solvent resistance, causing it to deform in 
the presence of strong organic solvents. Therefore, materials with higher solvent 
resistance, such as glass [17] and silicon, [18] are used. The pros and cons of the 
application of different microfluidic chip materials in droplet microfluidic system 
are shown in Table 1. The application of droplet microfluidics in the formation 
of functional materials for removal of pollutants will be discussed in Section 3 
subsequently. 


2.2 Microfluidic reactors 


Microfluidic reactors, also known as microreactors, have been widely used 
in wastewater treatment because the development of microfluidics technology 
in this area helps to overcome some existing problems in bulk reactors. The 
two essential issues of bulk reactor are photon transfer limitations and mass 
transfer limitations [20]. Thus, great attention and interest have been shown in 
microfluidic system as microreactor inherits the merit of microfluidics. Recently, 
the technologies of microfluidics in advanced oxidation processes (AOPs) for 
wastewater treatment were studied and focused. AOPs involve the utilization 


Application Glass/ Elastomers Thermoset Thermoplastics Hydrogel Paper 
silicon 

Droplets Excellent Moderate Good Good N/A N/A 

formation 

Production High Medium High Low Medium Low 

cost to high 

Reusability Yes No Yes Yes No No 
Table 1. 


Pros and cons of different microfluidic chip materials for different applications [19]. 
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of hydroxyl radicals (OH ) or sulfate radicals (SO4' ) as a major oxidizing 

agent to effect water purification [21]. This is because these powerful radicals 

are extremely effective to destruct the organic and inorganic contaminants in 
wastewater and transform them to less or even non-toxic products [22]. The most 
popular AOPs that employ the technologies of microfluidics are photocatalysis 
and Fenton processes. The development and application of several AOP mecha- 
nisms will be discussed in the next sub-section. Then, a brief example of the 
application of different types of microreactors in water treatment and resource 
recovery is listed in Table 2. 


2.2.1 Photocatalysis 


Heterogenous photoassisted catalysis, known as photocatalysis is one of 
the examples of hydroxyl radical-based AOPs which is of particular concern in 
wastewater treatment. It involves the utilization of light for decomposition or 
mineralization of organic pollutants into innocuous product, such as carbon 
dioxide and water, in the presence of catalysts [31]. The application of micro- 
fluidic technology offers a great number of advantages in photocatalytic water 
treatment. Microfluidic structures have larger surface-area-to-volume ratio, 
typically in the range of 10,000-50,000 m*/m? [32], compared to bulk reactor in 
which the surface-area-to-volume ratio is typically below 600 m?/m? [33]. The 
surface-area-to-volume ratio can be much larger if nanoporous photocatalyst 


Microfluidic device Application Result Ref. 
design 
Optofluidic planar reactor Degradation of 30% of dye degraded within 5 min 23] 
methylene blue with a reaction rate constant two 
orders higher than bulk reactor 
Microcapillary reactor Reduction of methylene Reduction rate of dye increased 24] 
blue by >150 times compared to batch 
system 
Tree-branched centimeter- Degradation of volatile 9596 of pollutants (benzene, 25] 


scale reactor organic compounds toluene, ethylbenzene, m-p xylenes 


and o-xylene) degraded in «5 s of 


residence time 


Jet-aerated microfluidic Degradation of Clopyralid is eliminated effectively [26] 
flow-through reactor clopyralid as model after 1 hour under several conditions 
organic pollutant 

Microfluidic atmospheric- Degradation of >97% of dye degraded [27] 
pressure plasma reactor methylene blue 
Droplet microfluidic Extraction of lead (II) Pb (II) ion was selectively and [28] 
reactor completely removed from the 

simulated wastewater effluent 

within 2.00 s 
Microfluidic chip with Adsorption of copper Adsorption capacity of 42.08 mg/g [29] 
polymethyl methacrylate (II) is achieved 
(PMMA) plates 
Simple cross microchannel Quantification of Bacteria were accurately enumerated [30] 


microfluidic device 


bacterial cells in potable 
water 


within 15 min after fluorescent 
staining 


Table 2. 
Application of microfluidic technology in water remediation and resource recovery. 
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film is used. Thus, higher heat transfer performance can be achieved by using 
microreactors. Furthermore, the rate of reaction is significantly increased and 
consequently favors having higher throughputs [34]. With the enhancement 

of reaction rate, the reaction time is reduced. The time taken for degradation 
process in a microreactor takes only several to tens of seconds [35], whereas 

bulk reactor requires several hours [36]. Besides, microfluidic layer has short 
diffusion length, typically 10-100 um to ease the diffusion of organic pollut- 
ants to the reaction surface [31]. In addition, microreactors usually contain an 
immobilized photocatalyst film under the thin layer of fluid. This can ensure a 
uniform irradiation on the reaction surface, resulting in higher photon efficiency 
[37]. Microreactor has self-refreshing effect as the running fluid can refresh 

the reaction surface naturally. This helps to move away the reaction products 

and stabilize the photocatalysts. According to journal by Wang et al. [35], the 
photocatalysts in microreactor can hold several hundred runs of photocatalytic 
reactions, whereas the activity of photocatalysts in bulk reactors starts to degrade 
after 10 runs of reactions [38]. 


2.2.2 Other mechanisms in AOPs 


Electro-Fenton process is an efficient AOP that involves activation of hydrogen 
peroxide by metal salts, typically iron, to produce hydroxyl radicals [39]. Electro- 
Fenton process is extremely effective in water remediation of the effluents, which 
cannot be efficiently treated using biological technologies [39-41]. Besides, AOPs 
with plasma-based water treatment (PWT) have been widely studied as PWT have 
the potential to reduce organic contaminants in wastewater. The application of 
microfluidic technology in PWT gives the benefits of large surface-area-to-volume 
ratio and flow control, in low-cost and portable devices [27]. However, more 
researches and development are needed to validate PWT performance at 
macro-scale [42]. 


3. Droplet microfluidics for the production of micro- or nanofunctional 
sorbents 


Droplet-based microfluidics is formed through fabricating emulsions of uniform 
size. There are two approaches to produce emulsions, which are active and passive. 
Unlike actively controlled microfluidic devices, the breakup of discrete phase in 
continuous phase driven is controlled in a fully passive manner, which is caused by 
flow instabilities and hydrodynamic pressure without external actuations, such as 
mechanical, electrical, thermal, and magnetic method [43]. In this section, passive 
formation of emulsion-based microparticles is mainly discussed. Emulsions can 
be produced and manipulated with micro device of different geometries forming 
different sizes and morphologies [44]. Monodisperse emulsions are produced in 
laminar flow region and generated drop-by-drop where the over size, shape, and 
morphologies of micro droplets can be control precisely [45]. 

The superior properties of droplet microfluidics are advantageous for precise 
microparticle manufacturing for wastewater treatment. Micro and nanofunctional 
particles with various morphologies are prepared using templates. In general, 
microfluidic device forms highly monodisperse emulsified droplets and forms 
microparticles via solidification, while nanoparticles are formed using photochemi- 
cal, chemicals or physical methods [46]. The space structure of the emulsions is 
controlled while preparing the adsorbent with microfluidic devices. By changing 
the device structure of microfluidic device, complex structure of droplets can be 
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produced such as single, double, and multiple emulsion [17, 47]. These emulsion 
droplets are classified according to the structures and the evolvement of the emul- 
sion into different structures and morphologies. Single emulsion can be a template 
to solid particles including spheres and non-spheres, while double and multiple 
emulsions can be templates for Janus particles, microcapsules, vesicles, hollow 
spheres or core-shell spheres [48]. 


3.1 Emulsion template: single, double or multiple 


Single emulsions are droplets of one phase fluid dispersed in another immisci- 
ble phase fluid. Flow focusing, cross-flow, co-flow are systems that are frequently 
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Figure 2. 

Wastewater treatment by adsorption method onto microfluidic synthesized microparticles. (I) The adsorption 
process of heavy metal copper (II) ions using solidified chitosan microspheres (CS-MS) synthesized via cross- 
flowing PMMA-based microfluidics [29]. (II) Absorption of fluorescent dye molecules using PEGDA-graphene 
oxide (GO) microspheres [11]. (III) Microfluidic-based synthesis of porous and non-porous PLA, PLGA, and 
PCL microspheres [54]. (I) is reproduced with permission from Zhu et al., 2017, Copyright 2017, Springer 
Nature. (II) is reproduced with permission from Zhao et al. [11], Copyright 2016, Springer Nature. (III) is 
reproduced with permission from Amoyav and Benny [54], Copyright 2019, MDPI. 
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used to form monodisperse droplets, and the coefficient of variation of droplets 
is usually less than 5% [3]. There are five breaking modes in passive generation, 
which are squeezing, dripping, jetting, tip-streaming, and tip-multi-breaking. 
These five modes have its own characteristics, for example, the structure and 
component of the droplets can be changed to produce inorganic nanoparticles, 
metal particles, and polymer particles [49]. Double or multi-emulsions are drop- 
lets with smaller droplets encapsulated in larger drops [48]. These emulsions are 
produced with capillary micro devices that involve three fluid streams in different 
capillaries. Initially, single droplets are formed when inner fluid is sheared by the 
middle fluid, then double or multiple droplets are formed when the outer fluid 
pinched off the single droplets containing in the middle fluid. There are dif- 
ficulties of precisely controlling the shell thickness, aggregation, and secondary 
nucleation. 


3.1.1 Special-shaped particles 


Non-spherical particles have unique properties, and they are usually fabricated 
with many strategies such as seeded emulsion polymerization [50], template 
molding [51], and self-assembly [52]. However, a high quality, monodisperse, non- 
spherical particles with tailored geometries and shapes yet still difficult to produce 
using these methods. Droplets with different sizes and shapes in microfluidic 
channels are confined with microfluidic technologies for fabrication of non- 
spherical particles. The droplet will be deformed into ellipsoid, a disk, or a rod if 
the largest sphere can accommodate in the channel of a larger volume of droplet. 
Non-spherical droplet will be formed after they are solidified in the confined 
channel [53]. 


3.1.2 Spherical particles 


There are different types of spherical particles, such as polymer microspheres, 
inorganic microspheres, noble metal nanospheres, and semiconductor nanospheres 
[48]. Polymer microspheres are usually prepared through spray-drying, coacerva- 
tion, and emulsification [48]. Inorganic microsphere, which is the composed of 
titanium, silica, and carbon, have potential application in biomolecules, sensor, 
catalyst, and drug deliver. Noble metal nanosphere, such as gold, silver, and plati- 
num, has shape and size dependent properties. It is hard to obtain desired size and 
size distribution as this individual nanoparticle tends to precipitate and coagulate 
to lower the surface energy. For the synthesis of semiconductor nanospheres, the 
microfluidic reactor should be chemically and thermally stable; thus, the droplets 
and carrier fluid could be stable, non-volatile, non-interacting, and immiscible 
from ambient to reaction temperatures. Hitherto, these microparticles are widely 
applicable in biological, pharmaceutical, medical (such as tumor treatment, drug 
controlled-release, and multi drug loading), optical, electrical applications, and 
researches. Moreover, polymer microspheres and inorganic microspheres are of 
great interest due to their potentials in adsorption separation as adsorbent in waste- 
water treatment, as shown in Figure 2 and Table 3 [29]. 


4. Sorption performance of different functional micro-sorbents for 
pollutants removal 


Pollution management is now one of the most challenging issue facing modern 
societies. Due to the increasing population as well as industrialization of most 
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countries, some pollutants are being discharged into aquatic environment without 
further treatment. This has a negative impact on the environment. The advance- 
ment of microfluidic technologies has allowed the synthesis of functional sorbents 
with greater sorption capacity. This is because the structure of the sorbent can be 
easily modified during emulsion. Thus, microsorbents with different functional 
groups can be synthesized to remove certain pollutants. Chitosan is a material 

that is widely used to make adsorbent for pollutant removal due to its affinity in 
removing heavy metals. In this subsection, the sorption performance of different 
functional microsorbents will be discussed as well as the kinetic model and adsorp- 
tion isotherm. 


4.1 Application and sorption performance of chitosan-based microsorbents 


Chitosan is a natural polymer material that is found in abundance. It is made 
from the chitin of crustaceans and shrimps. Due to its affinity with heavy met- 
als, it is a material with great potential for biosorbent synthesis. For instance, 
polyethylenimine-chitosan microspheres are used to remove methyl orange 
and Congo red dyes. Based on empirical observations, the uptake of methyl 
orange dye ranges around 88-9796, whereas the uptake of Congo red dye ranges 
around 86-9696 [55]. Zhai et al. [57] and Dong et al. [2] both have reported 
that the sorption performance of chitosan microspheres in the uptake of the 
common textile azo dyes. Apart from synthetic dyes, chitosan microspheres 
were also synthesized to remove copper (II) ions, the sorption performance was 
observed to be roughly 38.52 mg/g [29]. With the addition of polyacrylic acid, 
the sorption performance increased significantly [56]. Besides, Lv et al. [64] 
studied the sorption uptake of copper (II) ion with using polyethylenimine- 
poly(glycidylmethacrylate)-chitosan microsphere. Microfluidic synthesized 
ion-imprinted chitosan microspheres and thiourea-modified chitosan were also 
used to remove copper (II) ions. These studies have proved that microfluidic 
synthesized chitosan is a promising biosorbent for water remediation. Moreover, 
its mechanical intensity, sorption performance, and equilibrium adsorption 
amount of emerging contaminants are highly enhanced as compared to those 
conventional methods. 


4.2 Application and sorption performance of non-chitosan-based microsorbents 


Aside from chitosan, other microfluidic synthesized materials such as graphene 
oxide, silicon-based organic polymer, carbon, and silica were also used to be studied 
for wastewater treatment [11, 61-63, 65]. Copic et al. [61] reported that sodium 
dodecyl sulfate (SDS) and Congo red dye can be removed using carbon nanotubes 
microspheres. Li et al. [62] investigated that the synthesized hollow silica micro- 
spheres offer a much higher storage capacity as compared to conventional hollow 
nanospheres. Moreover, the drug detoxification capability of the hollow silica 
microspheres containing ethyl butyrate was tested on iodine removal [62]. The 
sorption capacity of iodine by silica microspheres goes up to 95% removal uptake. 
Lian et al. [63] studied the removal of toluene using polydimethylsiloxane micro- 
spheres synthesized via needle-based microfluidic devices. Ren et al. [65] synthe- 
sized anisotropic Janus microparticles loaded with Fe3O; and MnO, nanoparticles 
for the adsorption of basic dyes in wastewater. The sorption performance of the dye 
uptake using Janus micromotors can be varied from 47 to 94%. Thus, the current 
adsorption isotherm, kinetic modeling studies, and sorption performance of each 
microfluidics generated micro-sorbents for different water applications have been 
listed in Table 4. 
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5. Conclusion 


Certainly, microfluidic technologies are a relatively new research with great 
potential for development to enable more cost-effective synthesis of functional 
sorbents. Currently, the state-of-the-art microfluidic reactors for water remediation 
and resource recovery are being implemented in small-scale applications. Example 
of microfluidic reactors includes microfluidic atmospheric pressure plasma reac- 
tor, which is used to degrade methylene blue dyes. Furthermore, the advances in 
microfluidic technologies have improved the production of micro sorbents using 
microfluidic technology. Many efforts are also pouring into researching different 
functional sorbents for removal of different pollutants in wastewater. The results 
of such research have culminated in the discovery the adsorption capacity and 
isotherm model of different micro sorbents, which gives insight on the suitable 
sorbents required for different pollutants. 

At present, the need of highly effective water remediation and resource recovery 
has brought about mounting interest in the research of microfluidic technologies. 
Many breakthroughs had been achieved in such research which enabled highly 
specific and controlled synthesis of functional sorbents using microfluidic systems. 
However, there are still many barriers that prevent the implementation of microflu- 
idic technologies on an industrial scale. The greatest challenge against microfluidic 
technologies is the economical aspect. Due to the highly specific characteristics 
of sorbents used in wastewater treatment, most microfluidic systems are highly 
specific as well. Sorbents used in different applications would require different 
microfluidic systems to synthesize. Thus, more research will need to be done in 
order for such problems to be overcome. Of course, environmental problem such as 
water remediation and resource recovery cannot be solved by microfluidic technol- 
ogy alone and will need other technologies to complement it. 
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Chapter 4 


Application of Microfluidics 
in Biosensors 


Jing Wang, Yong Ren and Bei Zhang 


Abstract 


This chapter reviews the up-to-date researches in the field of biosensors inte- 
grated with microfluidic techniques, most of which are publications within the last 
5 years. The features of these biosensors, their applications, challenges, and possible 
future research interests in this field are also reviewed. 


Keywords: microfluidics, biosensor, bioreceptor, lab-on-a-chip, bioaffinity, 
PDMS, PAD, paper-based microfluidics, electrochemistry, optics, surface plasmon 
resonance, colorimetric, fluorescent, cell culture, food safety 


1. Introduction 


Biosensors are defined, by Tudos and Schasfoort [1], as “analytical devices 
comprised of a biological element (tissue, microorganism, organelle, cell recep- 
tor, enzyme, antibody) and a physicochemical transducer. Specific interaction 
between the target analyte and the biological material produces a physico-chemical 
change detected by the transducer. The transducer then yields an analog electronic 
signal proportional to the amount (concentration) of a specific analyte or group of 
analytes” [1]. The features of biosensors include the bio-recognition unit and the 
transduction mechanism from biological signals to measurable electronic signals, 
e.g., color, current, voltage, capacitance, light intensity, wavelength, and phase. 
Major parameters to assess the performances of biosensors include the following: 


* High sensitivity. The sensitivity of a biosensor is always the first and one of 
the most important parameters in assessing its performance. The efficiency in 
capturing analytes, the specific characterization of the analyte, the capability 
of converting biological signals into electronic signals (or the response of the 
system), and the systematic and environmental noises both determine the 
sensitivity of a biosensor. 


High stability and repeatability. The stability of a biosensor refers to the 
capability of a biosensor in performing consistently and reliably under desig- 
nated environments, and the stability of a biosensor is especially important 
when assessing portable or wearable biosensors that usually are applied in 
scenarios involving varied temperature, velocity, humidity, pressure, lighting 
conditions, etc. The repeatability of a biosensor mainly refers to the long-term 
performance of a biosensor under the same conditions and is usually tested 
regularly in commercial biosensors in order to recalibration. 
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* Quick response or real-time analysis and diagnosis. Real-time analysis usually 
delivers more information than providing a final result, e.g., binding rate, reac- 
tion time, kinetics, and saturation conditions, which can serve for the analysis 
in the applications of biological and chemical reactions and drug analysis. 
Response time required to bind sufficient molecules upon the sensing surface 
is typically determined by diffusion, which can extend to hours and even to 
days to generate a signal above the background noise level. This applies funda- 
mentally to all sensors that accumulate and concentrate target molecules onto 
a transducer, including fluorophore-tagged molecules in microarray spots, 
label-free optical biosensors, and impedance-based sensors [2]. 


* Low consumption of sample volume. The samples for biosensors are usually 
with low volume due to the nature, e.g., tissue, antibody, and some biologi- 
cal samples are with low concentration or small molecular weight, and this 
enforces biosensors to perform with minimum sample consumption. 


* Ease of operation. The application of biosensors goes from laboratory-based 
research to commercially available devices for at-home use. The operation of 
biosensors should eliminate professional operation or understanding of the 
device, but simply involves collecting samples and reading results. 


* High throughput. Single-function biosensors are fading away from stage even 
with extremely low cost. Biosensors should be able to integrate all the good 
qualities mentioned in above bullet points together with the capability of 
multiple-tasking. 


In traditional clinical healthcare, interests are in high-quality biosensor for the 
measurement of physiological indices. With the development in the requests of 
Internet of Things and self-service techniques, the interests and emphasis in health- 
care transfers from clinical healthcare to family healthcare, e.g., long-term monitor- 
ing of chronic disease [3], disease prevention and early detection, reachable clinical 
services for remote districts, etc. Under the precondition of high quality, biosensors 
miniaturized as portable or wearable are emerging to meet the new trend of era, 
which promise a bright future in health management and digital health; researches 
on the biological and instrumental parts of point-of-care (POC) and lab-on-a-chip 
(LOC) techniques belong to this area. The potential applications of biosensors 
include real-time health monitoring, remotely synchronizing health data with 
medical personnel, patient management, POC disease diagnosis, big data statistics 
in health management, etc. [4]. Real-time health monitoring in domestic applica- 
tions enables patients to monitor the health status by themselves at home without the 
assistance of professional personnel at minimum cost. Health data synchronized with 
medical personnel spares the patients from transporting and waiting for outpatient 
services, saving time and reducing medical expenses. Based on the collected health 
data, a potential service of clinical diagnosis is possible in an artificial intelligent 
health system in the near future [5]. Meanwhile, these new applications impose more 
requirements on the researches and developments of biosensors, as stated below. 


* Low cost. For at-home applications, lower cost of biosensors is vital, so there 
are emerging researches on adopting cheap materials, simplifying sensing 


systems and adopting smart phones in data processing, etc. 


* Noninvasive collection of samples. Biosensors for sensing human samples, 
noninvasive sensing, is preferable, especially for everyday or frequent 
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measurements. Researches using human samples as saliva, tear, sweat, and 
urine are one of the hottest topics. 


* Miniaturization of the biosensing systems [4, 5], including sample pre- 


processing unit, sensing unit, data collection/processing system, and data 
displaying unit. 


* The design of integrated sensor chip. A sensor chip with multiple functions is 


preferable especially for sensing human samples. 


2. Challenges in biosensing technologies 


All the qualities mentioned above, which both researchers and industry are seek- 


ing for, raise multiple challenges, and we try to summarize and list below: 
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* Specific binding. The recognition of analytes should be specific only to the 


analytes which is not affected by other chemicals, molecules, or cells. This 

is significantly more challenging when the sample components are complex 
and mixed with various kinds of molecules. For example, the detection of one 
specific antibody in human blood sample should eliminate the effects of all 
other antibodies, cells, electrolytes, etc., the detection reflects only the concen- 
tration of this antibody, and the detection of one specific heavy metal ion ina 
real polluted water sample should be able to distinguish the reaction induced 
by all other ions. 


* Non-specific binding, i.e., biofouling, in some cases significantly introduces 


signal noise, drift, or delay in biosensors [6]. The most common method to 
reduce non-specific binding is to completely wash the binding surface with 
buffer after the binding processing is finished; thus the weak binding induced 
by non-specific binding could be eliminated to the minimum extent. 


* Properties of bioreceptors, e.g., concentration and alignment. Plenty of 


papers [3, 7-10] presented the protocols of surface activation, modification, 
and functionalization, but the protocols greatly depend on operation and 
environment. Even following exactly the same protocol, the coverage rate 

of bioreceptors and the repeatability of operation may vary a lot, due to the 
immobilization of multiple layers on the sensor surface, which usually involves 
linking layers for stable sensor surfaces like gold and silicon, so the surface 
treatment protocols should be tested before operating the binding events. So 
far quite limited number of papers presented theoretical and/or experimental 
analysis on the effect of bioreceptor alignment/orientation on the performance 
of biosensors. The linking layer molecules are usually randomly polarized and 
oriented, which can induce destructive interference and dramatically reduce 
the collective charge polarization [11], and this means that even the surface 
treatment protocols could be repeated and the alignment of bioreceptors is 
another parameter that will highly affect the outcomes of the binding events. 
The detection becomes less sensitive. Chu [12] proposed a method to homog- 
enize the orientation of the chemical linker on nanowire-based field-effect 
sensor by applying an external voltage on a metal plate about 1 mm above the 
chip surface at certain frequency while grounding the back gate electrode; thus 
the molecular conformation can be maintained for hours or longer, and this 
method has been tested and proven by the detection of DNA hybridization 


Advances in Microfluidic Technologies for Energy and Environmental Applications 


reactions with poly-15 T ssDNA, showing that the alignment process promotes 
the sensitivity by 10-fold. 


* Design of biosensor assay matrix. The effects of specific and non-specific 
binding on signal was tested and analyzed by Schneider [13, 14], which proves 
that for all the binding events, proper design of the sensor assay should be 
optimized, especially when the sample components is complex, for example, 
proper reference binding sites should be included in order to eliminate non- 
specific binding from different components. But this in another way increased 
the requirement in both the imaging capability of the biosensor and the data 
processing capability. Meanwhile, the surface treatment, modification, and 
functionalization [4, 7, 9, 10, 15, 16] will be more complicated and need to be 
tested and verified, and the complexity in sensor surface properties (e.g., vari- 
ous refractive indices of bioreceptors) requests better system compatibility 
and responsivity. 


* Low concentration target molecule within a low-volume sample, i.e., extremely 
limited number of analytes available for detection. For the example, in the 
research of POC and LOC, 20-50 uL finger prick blood contains over 20,000 
kinds of biomarkers of clinical interest at concentrations as low as 10 pg./mL, 
meaning that only 106-107 available biomolecules for one target [2]. 


3. What is microfluidics? 


The definition of microfluidics, given by Whitesides from Harvard University, is: 

It is the science and technology of systems that process or manipulate small 
amounts (10? to 10 5 liters) of fluids, using channels with dimensions of tens to 
hundreds of micrometers. It offers fundamentally new capabilities in the control of 
concentrations of molecules in space and time [17]. 

The material most commonly adopted for the fabrication of microfluidic struc- 
tures is poly-dimethylsiloxane (PDMS), which is optically transparent and able to 
support important microfluidic components, e.g., pneumatic valves; meanwhile 
there are other materials with research interests, e.g., polycarbonate, polyolefin, 
silicon, and glass. Recently paper-based (reviews by [3, 18-21], and research by 
[22-28]) and cloth-based [29] microfluidics are drawing more attention because 
of the low cost, easy fabrication, and lightweight which are essential properties for 
POC applications. The major features of microfluidics are the small consumption of 
liquid sample and tiny dimensions of structures, which have significant impact on 
the development of biosensors, so the integration of microfluidics into biosensing 
techniques complies with the development of the era and generates unique features 
in biosensors, e.g., trace level of sample at high sensitivity. 


4. Advantages of microfluidic-integrated biosensors 


Microfluidics provide a closed and stable biosensing environment so to 
improve sensitivity. For on-site portable biosensors, the effect of an open envi- 
ronment on sensing results hugely lowers the biosensor performance. By integrat- 
ing the microfluidic structures, sample processing and biosensing reactions are 
carried out within a closed and relatively stable environment, thus promising 
better sensitivity and reliability [30]. Taking the example of the application of 
solid-phase polymerase chain reaction (SP-PCR) in online molecular diagnosis, the 
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development of this technique is hindered by lack of sensitive and portable on-chip 
optical detection technology. Hung [24] proposed a LOC device which combined 
the solid-phase polymerase chain reaction with supercritical angle fluorescence 
(SAF) microlens array embedded in a microchip. He demonstrated a high sensitiv- 
ity of 1.6 copies/puL and showed comparable detection limit and linear range to 
off-chip detection using conventional laser scanner, and he stated this device as 

an on-chip highly sensitive and multiplexed pathogen detection with low-cost and 
compact optical components. 

Microfluidic channel can efficiently, accurately, and significantly reduce the 
sensing area. Simulations and experiment results have shown that reducing the 
sensing area could shorten sensing time and increase the sensitivity with smaller 
sample requirement, especially at lower target concentrations [7, 11]. An increase in 
sensitivity of two orders of magnitude has been reported by Li et al. [31]. Meanwhile, 
the distribution of binding events along the sensing surface could be heterogeneous, 
and this could be induced by the heterogeneous in bioreceptor coverage, different 
alignment of bioreceptors, nonuniform concentrations of targets within samples in 
laminar flow, and nonuniform temperature, pressure, or other physical parameters 
along the sensing surface, the heterogeneity can be eliminated to minimum. Reduced 
sensing area also means miniaturizing sample volume which is essentially valuable to 
low concentration targets and rare targets with limited access. 

Microfluidic structures are capable of integrating multiple functions within 
one device without introducing extra equipment or tools. For example, by 
designing and optimizing microfluidic channels, sample injection, pretreatment, 
and processing can be easily realized. Usually for biosensing, the modification of 
the biosensing surface is compulsory for specific binding of targets, and this is 
doable in microfluidic structures which are even more stable and more promising 
than manual operations. For the biosensing events, the volume and speed control 
of sample are achievable which provides more valuable information, e.g., binding 
affinity, binding rate, kinetics, etc. 

Microfluidic devices are capable of automation. With or without external 
pumping system/equipment/tools, microfluidics is capable of integrating sample 
pre- and post-processing, sensing, surface modification, temperature control, EM 
field control, etc. The automation of microfluidic-integrated sensors and structures 
is reviewed by [30, 32-34]. However, the automation of whole microfluidic-inte- 
grated biosensor as one device still seems quite challenging as the liquid handling 
in this field is usually more complex which could involve up to dozens of solutions 
and operations like filtering, centrifugation, etc., together with the activity of 
biological samples to be considered. Partially automized microfluidic-integrated 
DNA biosensors are reviewed by Ansari [35]. Joung [36] presented a novel lateral 
flow immunosensor (LFI) for microfluidic-integrated enzyme immunosorbent 
assay (EIA) in POC testing (POCT), a chemiluminescent LFI-based automatic EIA 
system, the operation of which does not require additional steps such as mechanical 
fluidic control, washing, or injecting. The key concept relies on a delayed-release 
effect of chemiluminescent substrates (luminol enhancer and hydrogen peroxide 
generator) by an asymmetric polysulfone membrane (ASPM). When the ASPM 
was placed between the nitrocellulose membrane and the substrate pad, substrates 
encapsulated in the substrate pad were released after 5.3 + 0.3 min. Asa proof of 
concept, the high-sensitivity C-reactive protein level in human serum was detected 
by this sensor. 

Microfluidics enables both separate and mutual processing of multiple 
binding assays with single or multiple samples simultaneously. For the detec- 
tion of single or multiple targets in a real or complex solution, the design of the 
binding assay usually involves more than one kind of bioreceptors, thus meaning 
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that the designed samples to flow over each bioreceptor spot could be different. 
The delivering of different kinds of samples in sequential orders can be realized by 
unique design of microfluidic channels, pneumatic valves [17, 37], and/or centrifu- 
gal forces [33, 34]. 

Microfluidic structures ensure the precise control over experimental condi- 
tions [38]. What can be precisely controlled by microfluidic structures include 
flow rate, sample volume, channel volume, channel height, reaction time, etc. 
Integration of sensors with microfluidic channels serves to reduce assay time by 
constraining the diffusion distance between the molecules in the sample and the 
sensor and to create laminar flow over the sensor to distribute target molecules 
broadly and uniformly [2]. 


5. The present of microfluidic-integrated biosensors 


Biosensors can be classified based on target recognition events and transduction 
mechanisms [4]. Based on the target recognition events, biosensor receptors are 
included. Based on the transduction mechanisms, biosensors can be classified into 
optical biosensor (Raman scattering [39-44], surface plasmon resonance (SPR) 

[6, 45-47], fiber Bragg grating [48-53], fluorescent [54-58], chemiluminescence 
[36, 59]), electrochemical biosensor [60-64], calorimetric biosensor [6, 22, 65-69], 
and piezoelectric biosensor [70-74]. 


5.1 Target recognition 


Biological targets to be detected by biosensors, especially for the detection of 
analytes holden by human beings/animals, could be divided into two kinds, i.e., 
physical parameters and physiological/biological targets. Physical parameters like 
the body temperature, blood pressure, heart rate, velocity, and location usually 
do not request a corresponding and unique bioreceptor on the biosensor, as these 
physical parameters usually can be detected directly by optical, electronic, and 
piezoelectric sensors. Analytes as physiological/biological targets, however, cannot 
be detected directly, because of the complex components in a real human sample, so 
bioreceptors are adopted for the specific recognition of these targets, including cell, 
antibodies, DNAs, aptamers, and molecularly imprinted polymers [4]. 

The most commonly adopted physiological fluids of human beings/animals are 
blood, which has to be collected in an invasive way, and fluids that can be collected 
in a noninvasive way, e.g., sweat, saliva, tears, and urine, can be used in the predic- 
tion and diagnosis of various diseases [75-77]. Comparing with other physiological 
fluids, saliva is the outstanding fluid with the advantages of easy accessing and large 
volume, but with a major disadvantage of large range of variability in components 
and concentrations depending on the extent of oral cleanliness; examples that have 
been experimentally verified are using human saliva for the detection of cytokine 
[78], dopamine [51], insulin [79], fetuin [80], bacterial load [81], cholesterol 
[25], and cortisol [82]; using tear for the detection of dopamine [83], proteomic, 
lipidomic, and metabolomic composition [77]; using sweat for the detection of 
cytokine [84] and proteomic [76]; and using urine for the detection of anticancer 
drugs [85], L-carnitine [86], Chlamydia trachomatis, and Neisseria gonorrhoeae [87]. 
Samples of sweat and tear have been significantly undeveloped until quite recent 
when flexible materials and flexible electronic techniques achieved some milestones 
[4]. Currently the most well-explored targets in human physiological fluids include 
electrolytes (e.g., K*, Ca”, Na*) and major metabolites (e.g., myocardial enzyme, 
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glucose, urea), which lack specification to diseases, indicating the general physi- 
ological conditions [4]. 


5.2 Transduction mechanism 


So far, optical biosensors deliver the best sensitivity among the three other 
kinds of biosensors; electrochemical biosensors are the most popular choice as 
commercial-potential biosensors because of the compact size, low cost, and accept- 
able sensitivity; colorimetric biosensors are with a distinguished advantage of easy 
operation at extremely low cost but with a major disadvantage of low sensitivity; 
while the researches on piezoelectric biosensors are quite limited comparing with 
three other kinds of biosensors. Some most up-to-date researches on all fours kinds 
of biosensors are presented below. 


5.2.1 Optical biosensors 


Surface-enhanced Raman spectroscopy (SERS) and surface plasmon reso- 
nance are the two powerful optical biosensors with a unique feature of label-free 
sensing, as the analytes need no pre-processing to be labeled before sensing 
events and thus eliminate the false-positive or false-negative biosensing results 
induced by the labels. The first commercial product of SPR biosensor appeared 
in 1990 by the company of Pharmacia (named Biacore afterwards). Since then, 
more than 1000 papers were published annually using commercial SPR biosen- 
sors [88]. Most of these commercial SPR biosensors are bulky and only laboratory 
based. The development of plasmonic-based biosensors in the field of POC was 
reviewed in [43] together with recent advances in surface chemistry, substrate 
fabrication, and microfluidic integration. Here we explore a bit wider which is not 
limited to POC but microfluidic-integrated biosensors. In most of the researches 
mentioned below, the microfluidic structure usually serves as the sample han- 
dling unit. 

Tunc I et al. [39] presented the molecular specificity of Raman spectroscopy 
together with self-assembled monolayer of metallic AuNPs to detect CA125 
antibody-antigen molecules. Highly enhanced electromagnetic fields localized 
around neighboring AuNPs provide hot-spot construction due to the spatial 
distribution of SERS enhancement on the CA125 proteins at nM concentration level. 

Carneiro M et al. [41] reported the detection of carcinoembryonic antigen (CEA) 
in SERS using two different bioreceptors for CEA, i.e., a molecularly imprinted 
polymer (MIP) and a natural antibody. The MIP acted as a pre-concentration 
scheme for the CEA, while the natural antibody signals the presence of CEA on the 
MIP platform. The MIP film was first incubated in the sample containing CEA and 
next incubated in SERS tag, which is gold nano-stars coupled to 4-aminothiophenol 
(4-ATP) as Raman reporter, so the MIP acted as a pre-concentration scheme for the 
CEA. Then the MIP was exposed to the natural CAE antibody. A sensitivity down to 
1.0 ng/mL was reported. 

Zhu JY et al. [89] presented a biosensor that can be used for clinical diagnosis. 
This biosensor is based on localized surface plasmon resonance integrated with a 
biomimetic microfluidic “adipose-tissue-on-chip” platform for an in situ label-free, 
high-throughput, and multiplexed cytokine secretion analysis of obese adipose 
tissue. It was stated that this system enables simultaneous measurements of 
pro-inflammatory (IL-6 and TNF-alpha) and anti-inflammatory (IL-10 and IL-4) 
cytokines secreted by the adipocytes and macrophages and identified stage-specific 
cytokine secretion profiles from a complex milieu during obesity progression. 
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In the research of [90], the plasmonic biosensor integrated the microfluidic 
unit for plasma separation, allows the in-line separation of plasma directly from 
the bloodstream without any pre-processing outside the device, and channels it 
to the active detection area, where inorganic cerium oxide nanoparticles function 
as local selective dopamine binding sites through strong surface redox reaction. A 
detection limit of dopamine was achieved at 100 fM concentration in simulated 
body fluid and 1 nM directly from blood without any prior sample preparation. 
This demonstration shows the feasibility of the practical implementation of the 
proposed plasmonic system in detection of a variety of biomarkers directly from 
the complex biological fluids. Li XK et al. [91] reported the plasmonic biosensor 
integrated a multifunctional microfluidic system with small-volume microchamber 
and regulation channels for reliable monitoring of cytokine secretion from indi- 
vidual cells for hours. 

Besides the traditional plasmonic materials, graphene has recently received 
more and more attention in the field of both labeled and label-free sensing, because 
of its ability to harness electromagnetic fields, strong light-matter interaction of 
graphene layer, and its highly tunable optical properties [40]. Liu HP et al. [40] 
simulated the detection capacity of the graphene plasmonic biosensor using three- 
dimensional finite difference time domain method. Numerical results showed that 
the maximum sensitivity and figure of merit of the biosensor are 333.3 nm/RIU and 
16.665 RIU, respectively. 

Fluorescence is the other powerful optical biosensor which labels analytes and 
promises high sensitivity and specificity in target recognition. Raducanu VS et al. 
[56] reported a direct fluorescent signal transducer embedded in a DNA aptamer for 
versatile metal-ion detection. This sensor embedded with guanine-rich DNA aptamer 
internally coupled with Cy3 fluorescent dye that measures directly the DNA confor- 
mational changes upon metal-ion binding. Our signal transducer is environmentally 
sensitive that is internally coupled to the DNA aptamer. Potassium ion concentration 
was successfully measured in a variety of aqueous and biological test samples. 


5.2.2 Electrochemical biosensors 


There are plenty of researches on electrochemical biosensors, and majority 
of the commercialized biosensors belong to this category. Here we only present 
the electrochemical biosensors integrated with microfluidics that possesses both 
miniaturized structure and high sensitivity. 

Electrode-based chemoelectrical biosensors are the most common ones. Usually 
a working electrode and a blank/reference electrode are designed in such biosensor, 
and the samples cover both electrodes and generate a measurable electrical signal. 
Mi SL et al. [92] reported a sensitivity up to 567 nA mM(-—1) mm(-2), and the limit 
of detection was 4.5 M (vs. Ag/AgCl as the reference electrode) in the detection of 
metabolic lactate concentrations in HepG2 cells cultured with cancer drugs. 

Evans D et al. [93] demonstrated a fully integrated microfluidic amperometric 
enzyme-linked immunosorbent assay prototype using a commercial interferon 
gamma release assay as a model antibody binding system. What is unique in this 
research is that the assay cell is based on a printed circuit board (PCB) and the micro- 
fluidic assay chemistry was engineered to take place on the Au-plated electrodes 
within the cell. All components were manufactured exclusively via standard 
commercial PCB fabrication processes. À detection limit comparable to high-end 
commercial systems and a short diagnosis time of 8 minutes were demonstrated. 

Silicon nanowire field-effect transistor is one of the most sensitive biosens- 
ing techniques, but it is limited to analytes that carry charges. Weakly charged or 
uncharged analytes can hardly be detected [11]. Evans D et al. [31] presented a 
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method of immobilizing bioreceptors on the silicon nanowire sensing surface only, 
comparing with the traditional methods in which a large surrounding substrate 

is also covered with bioreceptors, and it was proven that restricting the surface 
modification substantially improves the sensitivity. 

Besides silicon nanowire, copper nanowire is adopted in electrochemical biosen- 
sors [94]. In [94], microfluidic chip is coupled to copper nanowires for the fast 
diagnosis of galactosemia in precious newborn urine samples. Galactosemia is a rare 
disease that is diagnosed through the identification of different metabolite profiles. 
The specific detection of galactose 1-phosphate (Gal 1-P), galactose (Gal), and 
uridyl diphosphate galactose (UDP-Gal) confirms type I, II, and III galactosemia 
diseases. The detection is extremely fast which is less than 350 s, required negligible 
urine sample consumption, and displayed impressive signal-to-noise characteristics 
and excellent reproducibility. 

Oliveira MC et al. [95] presented an amperometric biosensor using a screen- 
printed electrode modified with carbon nanotubes and nickel ions for the detection 
of glucose, which is characterized by the chemical oxidation of carbohydrate by 
NiOOH. Under optimized conditions, a limit of detection 3.9 pmol/L and a limit 
of quantification of 13 pmol/L were reported. The effect of concomitant species 
such as ascorbic acid, dopamine, and uric acid was investigated, and this method 
was successfully applied for the determination of glucose in a commercial blood 
serum human (original and spiked) sample. What is unique in this research is that 
the microfluidic system was assembled on a 3D-printed platform constructed with 
acrylonitrile butadiene styrene and integrated with nine cotton threads, providing a 
stable flow rate. 


5.2.3 Colorimetric biosensors 


Plenty of reports are available on colorimetric biosensors integrated with 
microfluidics, e.g., [27, 28, 96-98]; most of the reports highlighted the features of 
cost-effectiveness and miniaturization. Different from three other kinds of biosen- 
sors, the materials adopted for the integrated microfluidic structures are usually not 
PDMS, but paper for the majority and cloth in some researches. Currently majority 
of the researches focus on the applications in food safety [22, 27, 96] and heavy 
metal detection [67, 99, 100] in aqueous environment. The researches in the appli- 
cation of biological analytes are quite limited, due to the natures of analytes and 
bioreceptors and the environment conditions in order to keep the activities of both 
analytes and bioreceptors. 

Fraser LA et al. [101] presented a malaria biosensor whereby aptamers are 
coated onto magnetic microbeads for magnet-guided capture, wash, and detection 
of the biomarker. A biosensor incorporating three separate microfluidic chambers 
was designed to enable such magnet-guided equipment-free colorimetric detection 
of PfLDH. The biosensor showed high sensitivity and specificity when detecting 
PfLDH using both in vitro cultured parasite samples and clinical samples from 
malaria patients. 


5.2.4 Piezoelectric biosensors 


The research of piezoelectric biosensor integrated with microfluidics is quite 
underdeveloped so far. Possible reasons could be the lower sensitivity, poor biocom- 
patibility, and complicated fabrication. 

Yamaguchi M. [82] proposed a mass sensor based on mechanical resonance that 
incorporates a disk-shaped mechanical resonator, a separate piezoelectric element 
used to excite vibrations in the resonator, and a microfluidic mechanism. Electrical 
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power is used to actuate the piezoelectric element, leaving the resonator free from 
power lines. This sensor was reported to be suitable to analyze the concentration of 
a salivary hormone, cortisol in human saliva samples. 


6. Future research interests 


Future possible research interests in the field of microfluidic-integrated biosen- 
sors are proposed as the following: 


* Exploration of materials for both microfluidics and nanofluidics in different 
application scenarios. Besides PDMS, the exploration of other materials, e.g., 
engineering polymers, traditional glass, silicon, or metal, in special applica- 
tions that requires high chemical stability, high thermal stability, unique 
optical properties, and/or special mechanical properties. 


* Fabrication of microfluidics and nanofluidics and structures (e.g., valves, 
mixers). 


* Microfluidics with high chemical and thermal stability for special applications. 


* Integration of microfluidics and nanofluidics with sensors to form complete 
and functional systems that require no professional operations and ease in 
applications, e.g., LOC, etc. 


* Integration of microfluidics with data processing algorithms. The application 
of machine learning in sensing data processing could enhance the performance 
of biosensors in specialized environments. 


* Integration of microfluidics with communication techniques. Synchronization 
of sensing data with relevant users, remote control of the biosensors, and big 
data analysis of special sensing networks can be realized. 


7. Conclusions 


The state-of-the-art advances in biosensor development based on microfluidic 
technology have been reviewed in the book chapter with focus on the applications, 
challenges, and possible future research interests for each type of biosensor. It can 
be envisioned that microfluidic-based biosensors will remain a hot topic of inves- 
tigations because of the ever-increasing demands in various applications ranging 
from industry to biomedical detection. The interests in microfluidic-integrated 
biosensors promise even more prospective future in these areas. It is applications in 
wearable biosensor; portable biosensor can be explored in the future with enhanced 
sensitivity, improved stability, and miniaturized structure. 
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Abstract 


Clinical neurosurgeons used micro grinding to remove bone tissues, and drip 
irrigation-type normal saline (NS) is used with low cooling efficiency. 
Osteonecrosis and irreversible thermal neural injury caused by excessively high 
grinding temperature are bottleneck problems in neurosurgery and have severely 
restricted the application of micro grinding in surgical procedures. Therefore, a 
nanoparticle jet mist cooling (NJMC) bio-bone micro grinding process is put for- 
ward in this chapter. The nanofluid convective heat transfer mechanism in the 
micro grinding zone is investigated, and heat transfer enhancement mechanism of 
solid nanoparticles and heat distribution mechanism in the micro grinding zone are 
revealed. On this basis, a temperature field model of NJMC bio-bone micro grinding 
is established. An experimental platform of NJMC bio-bone micro grinding is 
constructed, and bone micro grinding force and temperatures at different measur- 
ing points on the bone surface are measured. The results indicated that the model 
error of temperature field is 6.796, theoretical analysis basically accorded with 
experimental results, thus certifying the correctness of the dynamic temperature 
field in NJMC bio-bone micro grinding. 


Keywords: grinding, nanoparticle jet, mist cooling, biological bone, thermal injury, 
temperature field 


1. Introduction 


As a precise material removal method, grinding has been extensively applied to 
bone tissue removal and surface treatment in surgical procedures by conforming to 
the combining trend of machine science and biomedical science. In the machining 
field, the method of using miniature grinding tools (diameter is generally below 
1 mm) for machining of feature size of 1~500 um is called micro grinding [1]. As 
bone tissues processed in osteological surgery are viable tissues and diamond grind- 
ing tool has obviously higher specific energy in the grinding process than other 
cutting methods, human bone, nerves, and blood vessels can be easily influenced by 
high temperature, and when grinding temperature is higher than 50°C, bone tissues 
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will experience irreversible necrosis, and nervous tissues will start experiencing 
thermal injury at 43°C [2]. In the present neurosurgery, surgeons usually perform 
drip irrigation of normal saline (NS) in the grinding zone following the following 
principle: The heat is carried away through convective heat transfer so as to cool the 
operative region. However, the drip irrigation-type cooling efficiency is not ideal; 
thermal injury in the operation process leads to necrosis of surrounding tissues and 
nerves, so the operation fails. The irreversible thermal injury of bone tissues trig- 
gered by high temperature in micro grinding is a bottleneck problem in bone micro 
grinding operation. Moreover, a large quantity of cooling liquid needs to be dripped 
into the grinding zone under drip irrigation-type cooling, which can easily reduce 
the visibility of the operative region under an endoscope [3]. 

In the machining field, as flood cooling grinding harms worker’s health and the 
environment, dry grinding can easily result in serious burn of the specimen, and 
microdroplet jet grinding has insufficient cooling performance [4-7]; researchers 
have put forward nanoparticle jet mist cooling (NJMC) technique according to the 
theory of heat transfer enhancement, which cannot only enhance the convective 
heat transfer in the grinding zone by taking full advantages of high specific surface 
of nanoparticles and heat capacity but, meanwhile, can improve the tribological 
properties of the grinding zone and can reduce heat generation by virtue of superior 
antifriction and anti-wear characteristics and high bearing capacity of nanoparticles 
[8-10]. On this basis, as it is the fact that clinical surgical bone grinding has 
disadvantages of thermal injury and low visibility of operative region and that 
nanofluid microdroplet grinding has good cooling and lubrication effect with a 
small use level of cutting fluid in the machining field, the research team where the 
author joined proposed NJMC bio-bone micro grinding process: adding medical 
solid nanoparticles of a certain proportion into the NS and selecting the 
corresponding surface dispersant according to physical and chemical properties of 
nanoparticle material, supplemented by ultrasonic vibration, so medical 
nanoparticles can enjoy uniform and stable distribution in the NS so as to form a 
medical nanofluid under stable suspension, which is then sprayed into the focus 
grinding zone in jet form as carried by high-pressure gas via a miniature nozzle 
[11, 12]. And this is expected to lower the bone grinding temperature while 
improving the visibility of the operative region. 

In recent years, researches have initially explored into thermal injury problem 
existing in the bio-bone grinding process. Taking heat transfer problem and its 
inverse problem in the bone grinding process with a miniature spherical grinding 
tool, Zhang et al. [13] conducted an in-depth study combining numerical simulation 
and experiment and analyzed the transient temperature field in the bone grinding 
process through a numerical simulation based on a heat generation model of bone 
grinding. Directing at the complex structure of bone issues and taking grinding 
temperature, grinding force, and blocking of grinding tool in the bone tissue grind- 
ing process as the study objects, Zhu et al. [14] obtained the grinding tool blocking 
mechanism in the bone tissue grinding by analyzing action rules of grinding tem- 
peratures on grinding temperature, grinding force, and grinding tool blocking. 
Sasaki et al. [15] analyzed heat generation during grinding using diamond tools with 
a diameter of 5 mm at a rotational speed of 100,000 r/min. They found that 
grinding for 10 cm along the grinding temperature without cooling liquid will 
expand the area of higher than 37.7°C continuously. To prevent thermal injuries to 
surrounding nerve tissues, Enomoto et al. [16] invented a new diamond spherical 
tool that can effectively reduce the temperature at the grinding zone. The surface of 
the diamond tool contains adsorbed super-hydrophilic TiO? nanoparticles, which 
effectively prevent temperature rise during bone grinding. The grinding zone is 
cooled during neurosurgery commonly by dripping a room-temperature NS into it. 
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The NS dissipates heat through natural heat convection, cooling the high- 
temperature region. 

For the bottleneck problem in the current clinical neurosurgical bone operation, 
namely irreversible thermal injury, a NJMC bio-bone micro grinding process was 
put forward. The convective heat transfer mechanism of nanofluid in micro grind- 
ing zone was studied. A dynamic temperature field model in NJMC bio-bone 
grinding was constructed. By reference to the suppression measures of grinding 
thermal injury in the field of mechanical engineering, the bottleneck in orthopedic 
operation was solved so as to provide theoretical guidance and technical support for 
avoiding or reducing thermal injury in clinical orthopedic operation. 


2. Grinding temperature field 


Starting from the universal temperature field problem, the heat quantity Q at 
one spatial point takes place transiently, and the temperatures at other points 
nearby this point will be changed under the action of the transferred heat quantity, 
and meanwhile, they are variable due to the spatial-temporal change. Temperature 
field is a generic term of temperature distribution at spatial points at a certain time. 
Generally speaking, temperature field is a function of space and time [17]: 


T =f (x,y,z, ¢) (1) 


Equation (1) denotes a three-dimensional (3D) unsteady-state temperature field 
where the object temperature is changed both in direction x, y, and z and at time t. 
In this study, temperature field is divided into two types: 


1.It is divided into steady-state temperature field and transient-state 
temperature field according to whether temperature changes with time. As 
indicated by Eq. (1), if the temperature field does not change with time, 
namely 27 — 0, it is a steady-state temperature field, or otherwise it is a 
transient-state temperature field. As mobile heat source is loaded, the 
temperature at specific point on the bone surface presents time-dependent 
change. Hence, the temperature field studied in this chapter is transient-state 
temperature field. 


2.It is divided into temperature field in cut-in zone, that in steady-state zone and 
that in cut-out zone according to the grinding tool/specimen contact arc 
length. As shown in Figure 1, the effective cutting part of the grinding tool is 
totally within the length of the specimen material, the grinding tool/specimen 
contact arc length is le, and the temperature field is temperature field in 
steady-state zone. As the grinding process starts, the grinding tool/specimen 


Micro 
grinding tool 


Bio-bone 


Figure 1. 
Schematic diagram of abrasive/bone contact states. (a) Cut-in zone. (b) Steady-state zone. (c) Cut-out zone. 
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contact arc length is gradually enlarged, but not reaching le, and at the 
moment, the grinding tool/specimen material is in cut-in phase, so the 
temperature field is in cut-in zone. As the grinding tool starts moving out of 
the specimen length, the contact arc length is gradually reduced to 0, and at 
the time, the grinding tool cuts out the specimen material, so the temperature 
is in cut-out zone. 


2.1 The solution method of grinding temperature field 


The solving method of grinding temperature field includes analytical method 
and finite difference method. 


2.1.1 Solution of grinding temperature field by analytical method 


Following the heat transfer theory and energy conservation law, analytical 
method solves the temperature rise function based on all kinds of boundary condi- 
tions in actual grinding machining so as to obtain temperature values on specimen 
surface and each internal node. The advantage of the analytical method is that it 
cannot only obtain the functional relationship regarding temperature distribution 
but also analyze different influence factors related to temperature field and their 
influence laws on temperature field distribution. Based on the theory of mobile heat 
source proposed by Jaeger [18] in 1942, many researchers have established theoret- 
ical grinding heat models with the main idea being superposition method of heat 
source temperature field. In other words, grinding interface is regarded as surface 
heat source constituted by numerous linear heat sources while linear heat sources 
are treated as a combination of numerous tiny unit linear heat sources, and each 
tiny unit linear heat source is simplified into combined action of point heat sources. 
Therefore, the basis for superposition method of heat source temperature fields is 
the solution of this temperature field at any time after transient point heat sources 
in an infinite object instantaneously emit partial heat quantity. 


2.1.2 Solution of grinding temperature field by finite difference method 


As for solving the grinding temperature field, it can be complicated to use 
analytical method to solve even simple heat conduction problem. Grinding 
machining is itself of high complexity relative to other machining modes [19, 20], 
where input parameters of grinding temperature field are miscellaneous, abrasive 
particle distribution is irregular, grinding state of abrasive particle (plowing, slid- 
ing, cutting) is uncertain, cooling medium participates in convective heat transfer 
in the grinding zone, and surrounding airflow field of grinding tool has an effect on 
temperature. A large quantity of nonlinear coupling relations exists in the grinding 
process, so if any input parameter is changed, it will influence the follow-up deri- 
vation of expressions, and the solving of temperature field becomes even harder for 
the analytical method [21, 22]. Under this circumstance, the finite element method 
based on numerical method is a very effective method of solving heat conduction 
problem. It is only necessary to determine boundary conditions and initial condi- 
tions in order to conveniently calculate grinding temperature field, so it has been 
widely used by researchers at present. As the finite element method makes many 
hypotheses for grinding temperature field and boundary conditions, only specific 
built-in modules of specific software can be used [23], and the deviation of calcu- 
lated grinding temperature from actual temperature is large. The finite element 
method based on numerical method is another effective method of calculating 
grinding temperature field between analytical method and finite element method. 
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It can accurately calculate temperature field through a theoretical modeling of 
boundary conditions of temperature field (heat flux, heat distribution ratio, con- 
vective heat transfer coefficient, etc.) according to actual grinding conditions. 
There have been few reports on calculation of grinding temperature field via finite 
difference method, so this method will be hereby described in details. 


2.1.2.1 Basic principle of finite difference method 


The object is divided into finite grid cells. A difference equation is obtained by 
transforming the differential equation, and the temperature at each grid cell node is 
solved through numerical simulation. As shown in Figure 2, under the 2D heat 
conduction problem, the object is divided into rectangular grids along directions x 
and y according to the spacing between Ax and Ay. The node is defined as intersec- 
tion point of each grid line, p (i, j) denotes the position of each node, i is serial 
number at node along the direction x, j is serial number at node along the direction 
y, and the intersection point between object boundary and grid is defined as 
boundary node. The basic principle of this method is to replace differential quotient 
with finite difference quotient so as to transform the original differential equation 
into a difference Equation [24]. 


2.1.2.2 Establishment of differential equation of heat conduction 


In the grinding heat conduction problem, Fourier' law is the most fundamental 
heat conduction equation, namely the heat quantity passing through infinitesimal 
isothermal surface A within limited time interval t is Q, which is in direct propor- 
tion of temperature gradient “ and is contrary to the direction of temperature field: 


oT 
dQ On dAdt (2) 
For heat flux: 
oT 
qx = x (3) 


where qx is heat flux in direction x and T is temperature field in direction x. 

The differential equation of heat conduction can be obtained according to 
Fourier’s law and energy conservation law, and it is hypothesized that there is 
internal heat source in the object. Based on the above hypothesis, microelement 
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Figure 2. 
Schematic diagram of gridlines and nodes of finite difference method. (a) Grid element. (b) Grid node. 
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Figure 3. 
Schematic diagram of thermal conductive microelement. 


dV - dx is divided from the object under heat conduction. As shown in Figure 3, the 
three edges of this microelement are parallel to axes x, y, and z, respectively. 
Thermal equilibrium analysis of the microelement is conducted. It can be known 
from the energy conservation law that the net heat quantity transferred in and out 
of the microelement within time dt should be equal to increment of internal energy 
in the microelement, namely net heat quantity transferred in and out of the micro- 
element (I) and increment of internal energy in the microelement (II). 

Net heat quantity transferred in and out of the microelement and increment of 
internal energy in the microelement will be, respectively, calculated as follows. 

Energy equilibrium analysis is implemented for the microelement in Figure 3. 
The net quantity transferred in and out of the microelement can be obtained by 
adding the net heat quantities transferred in and out of the microelement from 
directions x, y, and z, respectively. Along the direction of axis x, the heat quantity 
transferred in the microelement within time dt through plane x is: 


dQ, = 4, dydedt (4) 
The heat quantity transferred out through plane x + dx: 


dQ i a. = q, a, dy dedt (5) 


zm 90. 
where q, 4, = x + 3: dx. 
Hence, the net heat quantity transferred in and out of the microelement along 
the direction of axis x within time dt is: 


à 
dQ, ediscere A dxdydzdr (6) 


Similarly, the net quantities transferred in and out of the microelement along the 
directions of axes y and z within this time are, respectively: 


oq 

dQ, - dQ,,a4, = — gore (7) 
O 

dQ, dy ice is dxdydzdi (8) 


The following is obtained by adding net quantities transferred in and out of the 
microelement in directions x, y, and z: 
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aq, Ay ð 
[25 (St qu a) dxdydzdt (9) 
Within time dt, the increment of internal energy in the microelement is: 
oT 
II = pc dxdydzdt (10) 


where p is density of point heat source heat conducting medium and c is the 
specific heat. 
The three-dimensional heat conduction model can be obtained by Eqs. (9) and (10): 


2 2 
PaaS oT a) (11) 


à ^"latab a 


where a; is thermal diffusion coefficient of point heat source heat conducting 
medium. 


2.1.2.3 Differential equation converted into difference equation 


This specimen is assumed as a rectangular plane, and it is discretely decomposed 
into a planar grid structure. Isometric spatial step length Ax = Az = Al is taken; two 
groups of equally spaced parallel lines are drawn to subdivide the rectangular 
specimen and the equation of parallel line: 


x =x; =iAl,i = 0,1, .., M,MAI = l, 
he (12) 
D = 2; =jAl,j = 0,1, .., N, NAIL = b, 

where x; and z; are coordinate value of transverse line i in direction x and 
coordinate value of vertical line j in direction z, respectively, and both lines consti- 
tute the difference grid; /,, and bw are specimen length and height, respectively; 
and M and N are natural numbers. 

The grid area of difference calculation is obtained through subdivision as 
shown in Figure 3. 

A set of difference equations is built based on second-order difference 
quotient, namely: 


OT... T(-c-Lj)-cT(i—1j)-2T(ij) 2 

a Ë j= AP + O(AI) 

OT, . T(ijr1)-T(ij-1)-2T(ij) 2 13 
zz (bl) = Af + O(Al*) (13) 
^ (i,j) = Trad T(i,j) 4 O(At) 


The difference equation of each node in the grid can be obtained: 


At{ky - [Taj +1) + Taj — 1) +k: [TE - 5,5) + Th — 5,5)]) 
PylwAl? 


Ty (i,j) (14) 


Tatli, j ) = 


2A z 
" i _ 2At(ka + ka) 
PwCwAL 
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2.2 Boundary condition 


Heat transfer means transferring heat quantity from a system to another 
system with three main forms: convection, heat conduction, and radiation 
[25, 26]. The grinding tool/specimen contact interface is cooled with grinding fluid 
in the grinding process, so the grinding fluid flowing at this interface will conduct 
heat transfer on the specimen surface, where the main heat transfer form is con- 
vective heat transfer. As the grinding fluid contacts the specimen surface, partial 
heat quantity on the specimen surface is brought away by cooling fluid via convec- 
tive heat transfer, and the residual heat stays on the specimen matrix and is trans- 
ferred inside it. Therefore, the cooling effect of the machining zone can be 
strengthened by cooling fluid or by enhancing its heat transfer performance [27]. 
The temperature of external medium and convective heat transfer coefficient on 
boundary are given as shown in Figure 4. The heat quantity at grinding interface is 
generated due to grinding action of grinding tool and brought away by external 
cooling medium. In the meantime, both cooling medium and adjacent node (i-1, 1) 
on the specimen surface contact this node (i, 1), so the heat quantity at this 
node will be transferred to cooling medium in convective way and also to adjacent 
node, and in the end, a stable temperature is reached at grinding tool/specimen 
interface [28]. 

The thermal equilibrium analysis of heat conduction in temperature field is to 
solve differential equation. The complete description of heat conduction process 
includes monodromy conditions and differential equation of heat conduction, 
where the former includes time condition, physical condition, geometrical condi- 
tion, and boundary condition [29]. The boundary condition specifies characteristics 
of heat transfer process on object boundary, namely reflecting conditions for inter- 
action between heat transfer process and surroundings, so initial condition (tem- 
perature at f = 0) and boundary conditions (set temperature in boundary region or 
input and output heat fluxes) must be set in the temperature field. Based on the heat 
transfer theory, thermodynamic boundary conditions are generally divided into the 
three following types [30]. 


(a) 


Figure 4. 
Schematic diagram of heat conduction model and convective heat transfer at grinding interface. (a) Heat 
conduction model. (b) Convective heat transfer at grinding interface. 
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2.2.1 The first kind of boundary conditions 


Boundary condition of this type is forced convective boundary condition, 
namely temperature value at boundary of the given object at any time, and it is also 
called Dirichlet condition: 

Tl, = Tw (15) 
where Tw is set temperature at boundary surface sp. If boundary temperature is 


kept unchanged, T,, is a fixed value. If boundary temperature changes with time, T 
is a time-related functional expression. 


2.2.2 The second kind of boundary conditions 


Boundary condition of this type refers to heat flux on object boundary surface in 
normal direction at any time, and it is also called Neumann condition. The relation- 
ship between temperature gradient and heat flux is obtained through Fouriers law, 
and this is equivalent to temperature change rate on boundary sfin normal direction 
at any time: 


oT Qs 
onse (16) 


Sf 


where Qw is heat quantity passing through boundary surface sf. Qw = 0 when the 
boundary is under thermal insulation. Q,, is a fixed value when heat conduction at 
the boundary is constant, and it is a time-related function when the heat conduction 
changes with time. 

Differential equation at the second type of boundary: 


ð oT ð oT ð oT 
er) rm) +5 (be T) = Oe (17) 


2.2.3 The third kind of boundary conditions 


Boundary condition of this type denotes convective heat transfer between 
boundary surface and surrounding medium, and it is also called Robin condition. It 
can be known from Newton’s law of cooling that convective heat transfer occurs 
between boundary layer of the specimen and cooling heat transfer medium: 


oT 
h(T|;-o — To) — bx P" Iw (18) 


where qw is convective heat flux at the boundary surface sc between cooling heat 
transfer medium and specimen and h is convective heat transfer coefficient at the 
boundary between cooling heat transfer medium and specimen. 

Differential equation at the third type of boundary: 


[7 oT ð oT [7 oT 
(n) +3 m) +2 (te Z) - A(TI,, — Tr) (19) 
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2.3 Heat distribution coefficients 


An important problem exists in the establishment process of temperature field 
model in the grinding zone. It is necessary to determine the proportion of grinding 
heat quantity transferred into the specimen, namely heat distribution coefficient 
Rw. Under non-dry grinding condition as shown in Figure 4, the total heat quantity 
generated in the grinding zone is [31]: 


Fws 
Qtotal — "Ns = Wy +F dg $ dc oT qf (20) 


where b, is width of grinding tool; qw is heat quantity staying on the specimen 
surface; q, is heat quantity transferred into abrasive particle of grinding tool; q< is 
heat quantity carried away by grinding chips; and qfis heat quantity transferred out 
by cooling medium. 

The heat distribution coefficient model proposed by Rowe in consideration of 
convective heat transfer in the grinding zone is generally applied. As indicated by 
the model, under high-efficiency deep grinding, total energy in the grinding zone is 
transferred into grinding tool, grinding chips, grinding fluid, and theoretical math- 
ematical model of specimen material, and furthermore, its feasibility has been 
proved by a large quantity of experimental results. 

Heat flux transferred into specimen, grinding tool, grinding fluid, and grinding 
chips is correlated with parameters like maximum contact temperature Tmax» boil- 
ing point Tp of grinding fluid, and melting point Tm of specimen (Tmax < Tp) as 
shown in the following formula [32]: 


Iw = hw ` Tmax 

d; = hg Tmax 65 
qr = hg - Tmax 

4, = hac Tm 


where hw, hg, he, and hg are heat transfer coefficients of specimen material, 
grinding tool, grinding fluid, and grinding chips, respectively. 

The proportion of heat quantity flowing into the specimen can be obtained as 
follows: 


Rw = hw ^ Jas ded (22) 


3. Experimental study on micro grinding of biological bone with 
nanoparticle jet mist cooling 


3.1 Nanoparticle jet mist cooling bone micro grinding experimental platform 


Figure 5 shows the established experimental platform of NJMC bio-bone micro 
grinding, including feed system, fixed system, cooling system, and measuring sys- 
tem. Both axes x and y of the 3D displacement device are Shinano stepper motors, 
and axis z is Panasonic servo braking motor. A 120# diamond grinding tool is used, 
where the diameter of grinding head is 1 mm. The mist cooling supply device is mist 
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Figure 5. 
Nanoparticle jet mist cooling bone micro grinding experimental platform. 


cooling supply system produced by Shanghai KINS Energy-Saving Technology Co., 
Ltd. The jet parameters are uniformly set in the experiment as shown in Table 1. 

The clamping mode of thermocouple: a blind hole drilled on the back face of the 
bone specimen has a certain distance from the bone surface (grinding depth), the 
thermocouple wire is placed into the blind hole, the bone specimen is fixed on the 
dynamometer, and the force and temperature in the micro bone grinding process 
are simultaneously measured. Thermocouple type is K (TT-K-30), and measuring 
frequencies of both dynamometer and thermocouple are 100 Hz. 

In the biomedical field, HA, SiO?, and Al,O3 nanoparticles feature nontoxicity 
and good biological compatibility and are commonly used drug carriers in nanodrug 
release system [33, 34]; NS is a common clinically used cooling medium as osmotic 
pressure is basically equal to osmotic pressure of human plasma. Therefore, HA, 
SiO2, and ALO; nanoparticles with a diameter of 50 nm were used as solid nanoscale 
additives and NS as a nanofluid base to prepare HA, SiO», and ALO; nanofluids. 
Polyethylene glycol 400 has been extensively applied to lubrication in colonoscopy 
and gastroscopy by virtue of superior lubricating property and nontoxicity, and its 
safety in human body has been clinically certified. In the meantime, PEG400 also 
has good dispersity. Therefore, PEG400 was used as dispersing agent in this 


Jet parameters Value 
Nozzle diameter do 1mm 
Compressed gas pressure Ap - 0.5 MPa 
Feed flow rate Qs = 50 mL/h 
Spray inclination angle p-15* 
Nozzle height H = 0.6 mm 
Spray cone angle dz2/7 
Table 1. 


Table of mist cooling jet parameters. 
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chapter. When the dispersing agent with volume fraction of 2 vol.% and 2 vol.% is 
used, the suspension stability of the nanofluid will be the best. Hence, the nanofluid 
preparation method in this chapter was “two-step method,” namely adding 2 mL of 
HA, SiO;, and Al,O3 nanoparticles and 0.2 mL of PE in 100 mL of NS supplemented 
with ultrasonic vibration for 15 min. 


3.2 Temperature field of bio-bone micro grinding with different cooling 
methods 


The dynamic temperature field model of NJMC bone grinding was verified. 
Different cooling modes were adopted: dry grinding, mist, and NJMC. As dynamic 
heat flux was loaded to theoretically calculate grinding temperature field, tempera- 
tures at different measuring points on the surface of the bone material were mea- 
sured. As shown in Figure 6, three groups of thermocouples were used to 
simultaneously measure the temperatures at three measuring points—T1, T2, and 
T3—with spacing of 5 mm. 

Figure 7(a) shows typical grinding force signals measured under dry grinding 
condition, where F, is the force used to generate heat quantity. 

The temperature field of bone micro grinding is solved under initial condition 
(room temperature) of To = 20°C. Figure 8(a) shows temperature curves measured 
at different measuring points under dry grinding condition and theoretical temper- 
ature curves. It can be known that in the bone grinding temperature curves, the 
peak values at three measuring points are 36.1°C, 38.2?C, and 36°C, respectively, 
and measured peak values at the three measuring points are 36.7°C, 38.5°C, and 
36.6°C, respectively, namely the temperature on bone surface is ever-changing. 

Figure 8(b) and (c) show peak values measured through the experiment at 
three measuring points under mist and NJMC conditions and those obtained 
through calculation. It can be known that in comparison with average temperature 
(36.8°C) on bone surface under dry grinding condition, the temperatures on the 
bone surface under mist and NJMC decline by 13.696 and 33.996, respectively, thus 
proving the superior cooling effect of NJMC mode. The temperature error is 6.796, 
and theoretical analysis basically accords with experimental result, which verifies 
the correctness of the dynamic temperature field in NJMC bio-bone micro 
grinding [35]. 


~ Micro-grinding tool 


Thermocouples 


Figure 6. 
Schematic diagram of thermocouple clamping at different measuring points on bone surface. 
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Figure 7. 
Typical force signals in micro grinding of biological bone materials. (a) Typical force signal of bio-bone micro 
grinding. (b) Grinding force in y direction. (c) Fitting curve of grinding force in y direction. 
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Figure 8. 
Temperature curves and values of different measuring points on bone surface under different cooling conditions. 
(a) Dry grinding temperature. (b) Mist cooling grinding temperature. (c) NJMC grinding temperature. 


In order to explore into the dynamic characteristics of temperature field in 
bio-bone micro grinding, temperature fields in cut-in zone, steady-state zone, and 
cut-out zone are, respectively, analyzed. Figures 9 and 10 show temperature fields 
and temperature curves in bone grinding under NJMC conditions. The detailed 
analysis process is as follows: 


1. Cut-in zone: As shown in Figures 9(a) and 10, as the grinding starts, the 
grinding tool starts contacting the material and gradually cuts into the 
material, undeformed cutting thickness is gradually increasing [36, 37], and 
the heat quantity generated at grinding interface starts migrating into the 
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Temperature curves of cut-in, steady-state, and cut-out zones in bone micro grinding. 
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specimen surface. The temperature is low in the grinding zone, because the 
volume of material participating in the grinding in the cut-in zone is small with 
a small energy consumption. 


2. Steady-state zone: In Figures 9(b), (c) and 10, undeformed cutting thickness 
is kept at average value, and surface temperature no longer continues to rise, 
namely the formed temperature field reaches a steady state. As for grinding 
temperature fields calculated by researchers previously, after constant heat 
flux is loaded, the temperature value under which a steady state is reached will 
not be changed. However, in the temperature field of bone micro grinding 
with loaded dynamic heat flux as shown in the figure, the temperature value is 
ever-changing after the steady state is reached [38, 39]. 


3.Cut-out zone: As shown in Figures 9(d) and 10, according to the heat transfer 
theory, undeformed thickness is gradually reduced when the grinding tool is in 
cut-out zone, and the volume of specimen material participating in the 
grinding action is continuously reduced. If the heat quantity generated at 
grinding interface remains unchanged, the volume of the material to which 
heat is migrated in the cut-out zone is continuously reduced. As heat 
conductivity coefficient of air is extremely low, more heat is aggregated in the 
grinding zone, and then grinding temperature rises to a great extent. 


3.3 Effect of nanoparticle size on bone micro grinding temperature 


To probe into the influence laws of nanoparticle size on bone micro grinding 
temperature, Al,O3 nanoparticles with particle sizes of 30, 50, 70, and 90 nm and 
NS were used to prepare nanofluids for bone grinding experiment and measure 
grinding force and temperature at measuring point T2 on the bone surface. 

Figure 11(a) shows theoretical temperature curves under different cooling con- 
ditions, where abscissa axis denotes dimensionless distance (x/l = 2vwt/l., 21 = le), 
namely the location of this point in the grinding arc zone [40, 41]. It can be known 
from the figure that at the entry end in the contact zone, grinding temperature 
abruptly increases, peak temperature deviates from central position of heat source 
on the curve and deflects toward exit end of the contact zone, and the temperature 
at the exit end in the contact zone declines slowly. Figure 11(b) shows bone 
grinding temperatures measured using nanoparticles with different particle sizes, 
and it can be known that the grinding temperature increases with nanoparticle size. 
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Figure 11. 
Micro grinding temperature with different nanoparticle size. (a) Theoretical calculated temperature value 
(b) theoretical calculated temperature value. 
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Figure 12. 
Temperature values measured with different concentration of nanoparticles. 


3.4 Effect of nanoparticle concentration on bone micro grinding temperature 


SiO; nanofluid with volume fraction of 0.5, 1, 1.5, 2, and 2.5 vol.% was 
prepared in the experiment to investigate the influence laws of nanoparticle 
concentration on bone micro grinding temperature. Mist cooling was taken for 
comparative experiment and to measure grinding force and temperature at 
measuring point T2. 

As shown in Figure 12, the temperature measured under mist cooling is 32.7°C. 
Mist cooling is taken for control, and the surface temperature obtained using 
nanofluid with volume fraction of 0.5, 1, 1.5, 2, and 2.5 vol.% declines by 14.1%, 
17.1%, 19.6%, 22.9%, and 33.3%, respectively, namely the surface temperature in 
micro grinding decreases as the nanoparticle volume fraction increases. 


4. Conclusions 


As irreversible thermal injury and poor visibility of operative region exist in the 
current clinical bone micro grinding operation, a NJMC bio-bone micro grinding 
process has been proposed to investigate the convective heat transfer mechanism of 
nanofluid in the grinding zone and reveal the heat distribution mechanism. On this 
basis, a dynamic temperature field model in NJMC bio-bone micro grinding process 
has been established. An experimental platform of NJMC bio-bone micro grinding 
process has been set up to realize an experimental verification of dynamic temper- 
ature fields in this grinding process, followed by an experimental study of influence 
of laws of nanoparticle size and concentration on bone grinding temperature. The 
following conclusions are drawn: 


1. Heat distribution coefficient of specimen material, that of abrasive particle, 
that of grinding chips, and that of cooling medium, characterizes the abilities 
of heat transfer media to contend for heat quantity in the unit area of the 
grinding zone within unit time, so the heat distribution coefficient of the 
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specimen material can be expressed as the function of heat transfer coefficient 
of each heat transfer medium in the micro grinding zone. 


2.Compared with average temperature (36.8°C) on bone surface under dry 
grinding condition, bone surface temperatures under mist cooling and NJMC 
conditions decrease by 13.696 and 33.996, respectively, thus verifying the 
superior cooling effect of NJMC. The temperature error is 6.796, theoretical 
analysis basically accords with experimental result, and thus the correctness of 
temperature field in NJMC bio-bone micro grinding is proved. 


3. Al,O3 nanoparticles with different particle sizes and NS were employed to 
prepare nanofluids for bone grinding experiment. According to the 
experimental results, the grinding temperature rose abruptly at entry end of 
the contact zone, the peak temperature value deviated from the central 
position of heat source on the curve and deflected toward the exit end of the 
contact zone, the temperature at the exit end of the contact zone declined 
slowly, and the grinding temperature increased with the nanoparticle size. 


4. The bone grinding experiment has been carried out with nanofluids prepared 
using SiO; nanoparticles with different volume fractions and NS. In comparison 
with mist cooling, the surface temperatures obtained through nanofluids with 
nanoparticle volume fraction of 0.5, 1, 1.5, 2, and 2.5 vol.% declined by 14.1%, 
17.196, 19.696, 22.996, and 33.396, respectively, namely the surface temperature 
declined as the nanoparticle volume fraction increased in grinding. 
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Appendices and nomenclature 


NJMC nanoparticle jet mist cooling 


Tey set temperature at boundary surface 

NS normal saline 

Qw heat quantity passing through boundary surface 

3D three-dimensional 

dw convective heat flux at the boundary surface 

T grinding temperature 

h convective heat transfer coefficient at the boundary between cooling 


heat transfer medium and specimen 
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X,Y,Z 
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spatial coordinates 

width of grinding tool 

time 

heat quantity staying on the specimen surface 

contact arc length 

heat quantity transferred into abrasive particle of grinding tool 
the heat quantity passing through infinitesimal isothermal surface 
heat quantity carried away by grinding chips 

heat flux in direction x 

heat quantity transferred out by cooling medium 

density of point heat source heat conducting medium 
maximum contact temperature 

specific heat of point heat source heat conducting medium 
boiling point of grinding fluid 

thermal diffusion coefficient of point heat source heat conducting 
medium 

melting point of specimen 

isometric spatial step length 

heat transfer coefficients of specimen material 

specimen length 

heat transfer coefficients of grinding tool 

specimen height 

heat transfer coefficients of grinding fluid 

natural numbers 

heat transfer coefficients of grinding chips 

convection heat transfer boundary surface 

the proportion of heat quantity flowing into the specimen 
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Abstract 


Numerous researchers have developed theoretical and experimental approaches to 
force prediction in surface grinding under dry conditions. Nevertheless, the combined 
effect of material removal and plastic stacking on grinding force model has not been 
investigated. In addition, predominant lubricating conditions, such as flood, minimum 
quantity lubrication (MQL), and nanofluid minimum quantity lubrication (NMQL), 
have not been considered in existing force models. In this study, material removal 
mechanism under different lubricating conditions was researched. An improved theo- 
retical force model that considers material removal and plastic stacking mechanisms 
was presented. Grain states, including cutting and ploughing, are determined by cut- 
ting efficiency (/). The influence of lubricating conditions was also considered in the 
proposed force model. Simulation was performed to obtain the cutting depth (ag) of 
each “dynamic active grain." Parameter J was introduced to represent the plastic 
stacking rate and determine the force algorithms of each grain. The aggregate force 
was derived through the synthesis of each single-grain force. Finally, pilot experiments 
were conducted to test the theoretical model. Findings show that the model’s predic- 
tions were consistent with the experimental results, with average errors of 4.19% and 
4.3196 for the normal and tangential force components, respectively. 


Keywords: grinding, minimum quantity lubrication, nanofluid, material removal 
mechanism, force model 


1. Introduction 


With need of green manufacturing, nanofluid minimum quantity lubrication 
(NMQL) was introduced to grinding [1, 2], which could effectively improve cooling 
and lubricating performance compared to conventional dry or flood conditions 
[3, 4]. Furthermore, with the deeper understanding of grinding mechanism and the 
development of process technology for difficult-to-grinding material, advancement 
in the experiment has been rapidly made [5-7]. Of course, with such a demand from 
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designing of grinding technological parameter, material removal mechanism and 
force model under NMQL condition are of key importance. 

Material removal mechanism, here, mostly refers to deformation, fracture pro- 
cess, and strain rate of material in grinding zone, as well as mechanical behavior of 
chips and furrow formation [8, 9], which relates to grinding parameters and 
cooling/lubricating conditions [10, 11]. For dry grinding, the experimental studies 
have been carried out, and force, temperature, shape of chips, and surface quality 
were observed, which showed special phenomenon when higher speed (i.e., speed 
effect) and lower cutting depth (i.e., size effect) were used. Nevertheless, material 
removal mechanism under NMQL condition is vacant and important, which should 
be researched deeply. 

Force model of grinding, in general, was vector sum of sliding force, cutting force, 
and ploughing force of all dynamic active grains in grinding zone. Many researches 
were carried out to establish grinding force model under dry condition. Werner [12] 
derived a grinding force equation with two parts (cutting force and sliding force), 
which was adopted by Malkin and Hwang [13], Li et al. [14], Younis et al. [15], and 
Tang et al. [16]. To obtain grain stages and numbers, Rayleigh s probability density 
distribution [17, 18] and normal distribution [19] were used to describe cutting depth 
of grain [20]. Except for steel material, the grinding force model was also researched 
for single crystal sapphire [21], SiCp/Al composites [22], complex optical mirrors 
[23], etc. However, previous models were limited at three aspects: 


i. The average cutting depths were used for each grain states (i.e., cutting, 
ploughing, sliding), which cannot represent the actual situation in the 
grinding zone [24]. 


ii. The critical depth of cutting and ploughing was determined by an 
experience-based evaluation method: “when the grain-cutting depth 
reached 0.05 times of the radius of cutting grains, cutting action occurred" 
[25-27]. It's imprecise for different material. 


iii. The number of "static active grains" was presented to develop force models, 
which are not considered interference effect of grains on each other [28-32]. 


iv. The lubricating condition was not considered. 


This chapter proposes an improved theoretical force model that takes the mate- 
rial removal and plastic stacking mechanism into consideration. The scratch tests 
reveal the relationship between cutting depth (a5) and cutting efficiency (/) to 
distinguish the cutting and ploughing grains and determine the force algorithms of 
each grain. The tribological tests reveal the friction coefficients under different 
grinding lubricating conditions to calculate the frictional force of each grain. In 
addition, the a, of each "dynamic active grain" in the grinding zone is obtained by 
simulation, to assist the development of the grinding force model. To verify the 
proposed model, grinding experiments were performed, and force values were 
obtained and compared with predictive values. 


2. Material removal mechanism of NMQL grinding 
2.1 Deformation and strain rate 


The material deformation mechanism in the grinding zone with single abrasive 
particle is shown in Figure 1. The material deformation mechanism in the chip 
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Figure 1. 
Deformation mechanism of material in the cutting zone. 


Figure 2. 
Stress analysis of the workpiece material in the cutting zone. (a) Primary deformation zone and (b) secondary 
deformation zone. 


formation zone during the grinding machining is identical with that during cutting 
machining, and primary deformation zone and secondary deformation zone 
simultaneously exist. 

Ding et al. [33] conducted an analog simulation of high-temperature cutting of 
nickel-based alloy workpiece material with single abrasive particle and also 
observed the abovementioned phenomenon. Different from cutting machining, the 
negative rake cutting form of abrasive particle leads to a different mechanical action 
mechanism between abrasive particle/grinding chips; the calculation formulas of 
shear angle o, frictional angle Jf, and strain rate are also changed, but the impact 
trends of grinding process parameters and lubricating conditions on material 
removal mechanism are totally the same. The force analysis of the grinding zone of 
single abrasive particle is shown in Figure 2, and material force bearing and 
deformation mechanism can be described as follows: 


i. Inthe primary deformation zone, Figure 2(a) is force analysis graph of point A, 
and the resultant force of frictional force F;and extrusion force F, is F, 
According to the mechanical theory of material, the included angle between 
rupture chip formation direction (direction F, of shear zone) and resultant force 
F, is 1/4, and shear angle ¢ is on horizontal line; at the time, frictional force Ff 
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exerts a promoting effect on rupture chip formation of the material upward 
along the abrasive particle/cutting chip interface. Hence, the expression of shear 
angle ọ is solved according to geometrical relations of several angles as below: 


T 
91 = Zt Pr — Yo (1) 


The ratio of thickness to length of grinding chips formed through cutting 
machining of abrasive particle is about 1/100, so negative rake is yo ~ 0/2 and 
frictional angle is f; = arctan(j) under abrasive particle cutting, where y is frictional 
coefficient between rake face of abrasive particle/cutting chips. The relational 
expression between yo and f is substituted into Eq. (1) to obtain: 


m 


g — 4 *arctan x- ; 2) 
ii. In the secondary deformation zone, Figure 2(b) is force analysis chart, 
where the resultant force of frictional force Fs and extrusion force Fp is F,. 
According to the mechanical theory of material, the included angle 
between rupture chip formation direction (direction F, of shear zone) and 
resultant force F, is 1/4, and shear angle q is below the horizontal line; at 
the time, frictional force F exerts an inhibitory effect on flow of grinding 
chips downward along the abrasive particle/cutting chip interface. After 
kinetic analysis, speed of movement on abrasive particle/grinding chip 
interface is low due to frictional force, but the speed on free face of 
grinding chips is high; under high-speed cutting condition, transient speed 
difference is generated between abrasive particle/grinding chip interface 
and free face of grinding chips due to high strain rate of grinding chips, 
thus generating a shear zone and forming the secondary deformation zone. 
Hence, the expression of shear angle gy can be solved according to 
geometrical relations of several angles as: 
¢, =f, +y -Z —andango o -* (3) 
2 071 
For high-speed grinding machining, the change of strain rate is the primary factor 
influencing the material removal mechanism. The strain rate reaches as high as 10^— 
10? s^! in the cutting machining process, but for low-speed grinding machining (V, < 
30 m/s), the strain rate of material removal reaches up to 10’-10° s ! because of 
negative rake characteristic of abrasive particle. Therefore, it's speculated that the 
strain rate under high-speed grinding machining condition is higher than that under 
low-speed grinding machining condition of grinding wheel by 1 or 2 orders of magni- 
tudes. Under high-speed grinding machining condition, the strain hardening effect 
and strain strengthening effect caused by strain rate of material deformation generate 
a remarkable effect on the material removal mechanism. Jin and Stephenson [34] 
established the equation of shear strain in the shear zone of abrasive particle and strain 
rate. Taking maximum undeformed cutting thickness and shear angle in the primary 
deformation zone and secondary deformation zone into account, strain rate of grind- 
ing chips in the two deformation zones could be calculated by the following equations: 


Av, cos C sinf, ! arctany- .. ) 
| (4) 
m, NEA D cost F arctan jo 


ü T 0 
1 
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rM x " 
aV. Jo, cosfarctan se f- 4 


where q is shear angle, 0 is vertex angle of abrasive particle, 4; is average length- 
width ratio in the shear zone (4, = 6-12), v is cutting speed of abrasive particle, and 
v & V, is for single abrasive particle grinding. 

It can be known from Eq. (4) that factors influencing deformational strain rate of 
workpiece material during the cutting process with single abrasive particle include 
grinding parameters, abrasive particle shape, and lubrication characteristics of 
abrasive particle/grinding chip interface. In this study, as V./Vw = 1 x 10* and 
abrasive particle shape are unchanged, the influence trends of grinding speed V,, 
grinding depth Vw, and grinding coefficient jJ on shear rate will be hereby discussed. 

As shown in Figure 3(a), when other parameters remain unchanged, the strain 
rates in both primary shear zone and secondary shear zone present a linear growth 
trend when V, increases from 30 to 120. In comparison with the primary deforma- 
tion zone, the strain rate in the secondary deformation zone has a higher growth rate. 

According to Figure 3(b), when the frictional coefficient j gradually increases, the 
strain rate in the primary deformation zone presents an increasing trend of quadratic 
function, because the frictional force in the primary deformation zone exerts a gaining 
effect on material removal and the strain rate increases with frictional force, which is 
better for chip formation. In the secondary deformation zone, as the frictional coeffi- 
cient y increases from small to great (characterizing lubricating condition), the strain 
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Strain rate under different grinding parameters and frictional coefficient. (a) Influence of grinding speed, 
(b) influence of frictional coefficient, (c) influence of grinding depth and (d) two-factor analysis. 
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rate in the secondary deformation zone presents a declining trend of quadratic func- 
tion, the frictional force in the primary deformation zone prevents grinding chips 
from flowing out, and flow direction of grinding chips is gradually turned from 
towards abrasive particle into along the surface of abrasive particle, so the strain rate 
gradually declines. As shown in Figure 3(c), the strain rate presents a linear declining 
trend due to increase of grinding depth. As this study focuses on speed effect and 
lubricating effect in the high-speed grinding machining, in consideration of impact 
trends of grinding speed and frictional coefficient on strain rate, Figure 3(d) shows 
the variation trend graphs of the primary deformation zone and secondary deforma- 
tion zone under the influence of grinding speed and frictional coefficient. 


2.2 Mechanism of chip formation 
2.2.1 Theoretical research 


Impact dynamics describes material deformation behaviors under high speed as 
an adiabatic shear process. Adiabatic shear effect refers to constitutive instability 
(thermal viscoplastic instability) of the material under impact load. Extremely 
strain rate exists in the material removal shear zone under the action of impact load, 
and the shear zone is a thermal insulating environment within a very short time. At 
the time, nonelastic energy in the material deformation process is converted into a 
large quantity of grinding heat, which leads to abrupt temperature rise of the 
material in the shear zone and declination of material hardness (softening) under 
temperature rise. The softening effect exceeds strain hardening effect and strain 
strengthening effect, and thus shear zone is formed and material instability is 
caused. Hence, the chip formation process under the effect of abrasive particle is 
dynamic stress and thermal force coupling action process under a high strain rate. 

The removal process of metal material derives from plastic deformation after the 
material strain increases. In the previous studies, scholars have conducted a large 
quantity of dynamic material tensile tests [35], proving strain hardening and strain 
strengthening phenomena. As shown in Figure 4, the material removal process 
with the change of strain rate can be divided into the three following types 
according to the stress-strain curve: 

As shown in Figure 5(a), under quasi-static conditions, the material deformation 
stress will increase obviously with the strain, indicating that notable strain harden- 
ing effect happens to the workpiece material. Grinding chips will be formed after 
the stress increases to material breaking limit ob1, and free face of cutting chips 
presents periodic upheaval (bamboo shape), which is caused by extruding defor- 
mation in the grinding chip formation process; however, the strain rate of material 
deformation is low because of low speed, and no obvious shear zone is formed. 

As shown in Figure 5(b), under high strain rate, assume that the shear zone of 
grinding chips is in an isothermal environment: the material deformation stress in 
the strain hardening phase is far higher than that under quasi-static conditions, and 
this phenomenon is called “strain strengthening effect” of plastic deformation of 
material, and its deformational resistance is enhanced under strain strengthening 
effect. At the time, as abrasive particle has a strong impact effect on the material, 
workpiece material forms shear layers under high strain rate so as to form grinding 
chips, and rupture stress limit between shear layers is breaking limit oj? of the 
material under high strain rate. Free face of grinding chips presents notable periodic 
shear slippage layer, and the slippage distance of shear layer under high strain 
condition is obviously higher than that under quasi-static state; and the higher the 
speed, the higher the strain rate, and the greater the slippage distance of shear layer; 
and even radical rupture and separation between shear layers is caused by bluff- 
type rupture process. 
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Figure 4. 
Stress-strain curve under "velocity effect." 
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Figure 5. 
Chip-forming mechanism under (a) quasi-static, (b) high strain isothermal, and (c) high strain condition. 


As shown in Figure 5(c), in consideration of temperature rise at abrasive parti- 
cle/grinding chip interface under high strain rate, the plastic deformation process of 
the material at the time is a comprehensive result of strain hardening effect and 
thermal softening effect, namely, adiabatic shear process. In Figure 4, the stress- 
strain curve is between quasi-static curve (curve 1) and high strain isothermal 
condition (curve 2) under the action of adiabatic shear, because a large quantity of 
grinding heat softens the material so as to reduce deformational resistance under 
the action of strengthening strain rate. As the strain rate increases, the heat trans- 
ferred to shear zone is also increasing, thermal softening effect becomes more 
obvious, and the breaking limit relation of the material under three boundary 
conditions is 6,2 > 653 > Op1. Material rupture and plastic flow after thermal soften- 
ing simultaneously exist between shear layers of grinding chips. 

Hence, the critical research hotpots regarding material removal mechanism 
under high-speed grinding conditions are to explore into influences of grinding 
parameters and lubricating conditions on strain strengthening effect and thermal 
softening effect and build an evaluation model. 


2.2.2 Experimental research 


SEM of chips were shown in Figure 6. 

The grinding chip shapes are similar at the grinding speed of 40 m/s under three 
different working conditions, namely, continuous banding chips, because the 
material strain rate is low and thermal softening effect in the grinding zone is small 
under low grinding speed. When the grinding speed is 80 m/s, obvious shear zones 
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Figure 6. 
Chip morphology at different lubricating conditions. 


are observed on the free faces of the three types of grinding chips, indicating that 
under high grinding speed, the material experiences adiabatic shear during the 
forming process. Differently, the shear slippage layer spacing obtained under dry 
grinding conditions is small, that under MQL condition is larger, and that under 
NMQL is the largest, where even rupture takes place. 

On the one hand, as the lubricating effect is enhanced and frictional coefficient 
is reduced in succession under the three lubricating conditions, the strain rate in the 
primary deformation zone at the same speed increases successively, so does shear 
slippage distance under the three conditions. On the other hand, the output of 
frictional heat is reduced as the lubricating effect is enhanced successively under the 
three conditions, and the heat quantity transferred out of the shear zone is increased 
as the cooling effect is enhanced successively under the three conditions. Therefore, 
the thermal softening effect during NMQL cutting chip formation process is lower 
than that under dry grinding condition. The shear layer slippage of grinding chips is 
mainly material rupture, followed by plastic flow, so shear layer spacing is larger, 
and grinding chips even rupture. Under the grinding speed of 120 m/s, the 
abovementioned rule is magnified again. Banding chips are still formed due to 
thermal softening effect under dry grinding, rupture of grinding chips induced by 
excessive shear layer spacing also appears under MQL condition, and the rupture 
under NMQL condition is further aggravated. 


2.3 Mechanism of furrow formation and cutting efficiency 
2.3.1 Plastic stacking effect 


Cutting depth a, (varies from 0 to dgmax) and critical value ag. were the key 
points for distinguishing ploughing and cutting stage. The value a5. could be gotten 
from cutting experiment of single grain for different materials. The change rule 
could be concluded from Figure 7. 
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Plastic stacking Material unit 


Figure 7. 
Material deformation behavior in cutting of a single grain. (a) Ploughing grain and (b) cutting grain. 


i. When ag < ag, (ploughing stage). Material of workpiece flows upheaved 
and stacked at two sides of the furrows because the stress between material 
units is lower than fracture stress. 


m. 
= 


i. When ag > ag. (cutting stage). Plastic stacking theory [36], which describes 
the coexistence phenomenon of material removal and elastic-plastic flow in 
grinding, is a good method to understand deformation behavior of material. 
Pink-colored material units in Figure 7(b) were the elastic-plastic flow 
part, and blue-colored units were the material removal part. However, a 
boundary was born between two regions. 


Therefore, chip formation region (0—a4) and elastic-plastic flow region (a4-2/2) 
were defined in this paper, in which two were differentiated by critical value of 
angle a4, as shown in Figure 8. a4 (could be calculated by cutting efficiency f [37]) 
was different when material changes: 


a, — arccos(4! P) (6) 


2.3.2 Model of P( ag) 


Parameter / will increase from 0 to almost 1 when a, increases [38]. Therefore 
the model of f/(a,) is the key point of grinding force prediction. Cutting experiment 
of single grain was carried out for stainless steel material (440 C, Ra = 0.04-0.05 
um). Surface topography of furrow is shown in Figure 9, and numerical fitting of 
f(ag) was shown in Figure 10. 

The f/(a5) curve present S-shaped (when a, varies from 0 to 4.5 pm) and further 
linear trend (when ag > 4.5 um) . Therefore, Gompertz growth equation and 
Gaussian fitting were considered simultaneously, and piecewise function 
(R? - 0.9945) is shown as Eq. (7). 
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It could be observed from Figure 10 that (i) three stages could be learned from 
curve slope of (aș); (ii) compared to study of Hahn [39], sliding stage (0—0.023 
um) is declined; and (iii) the transition stage exists between ploughing and cutting 
stage, which is different from traditional theory. 


Plastic 
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Grain 


Figure 8. 
Plastic stacking mechanism. 


Figure 9. 
Result of scratch tests. (a) Surface topography of furrow, (b) 2D profile curve of furrow, (c) 3D graph of furrow. 
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Equation of cutting efficiency (f). 


3. Force model of single grain 
3.1 Stress state of grain 


As a premise for force model, analysis of stress state was carried out in this 
section according to the above results, as shown in Figure 11. Plastic flow is the 
material behavior, and equal stress value (yields stress 6,) distributes on grain 
surface. The stresses at chip formation region (ôo) and elastic-plastic flow region 
(à) show characteristics as follows: (i) equal stress value (fracture stress ô) dis- 
tributes on chip formation region and (ii) linear formula of stress distributes on 
elastic-plastic flow region, which varies from 0 to 6p. 


3.2 Force model in chip formation region (a,-7/2) 


ĉo was composed of 69; and ôo2 in chip formation region. As shown in Figure 12, 
the integral unit (ds) can be expressed as: 


a, stand 


d. = do (8) 


5 


A-vus E 


The plastic flow force equation F;(01) (a5) /F,.(01) (ag) in chip formation region 
could be calculated as: 


Fyn (4,) = [78, «a, tan0 cosa da (9) 
Fola) = f: à, a," tan 8-dar- = ay), a, tan 8 (10) 


where @ is grain vertex angle and a4 could be obtained by Eqs. (6) and (7), 
respectively. 

The chip formation force equation F;;(2)(a4) in chip formation region could be 
calculated as: 
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Figure 11. 
Schematic diagram of the grain stress state. 


Figure 12. 
Calculation principle of cutting force. 


Fas.) = Ds a," tan cosa di — Ò, A, (11) 
where A, is area of the chip formation region and 62 could be obtained by Eq. (11): 


ôn = dme, “by (12) 
2-(l—sinez,) 


The normal force could be further calculated by Eqs. (11) and (12): 


: zE- " 
C t.) 2ó,.-a,^- tan? 9-da = — 8, -a - Aun? 
Fal) = [s &,,-a,  -tan A-du Cy petuum à, -a, ^ lan’ 8 

In summary, ôo can be expressed as: 


z-lan8 


5, = Jy, +S =S 2 ———————- 
03 7 Og, t Oy, 70, 2-ü0-sina.) * 


(14) 


112 


Material Removal Mechanism and Force Model of Nanofluid Minimum Quantity Lubrication... 
DOI: http://dx.dot.org/10.5772/intechopen.91141 


3.3 Force model in elastic-plastic flow region (0— a) 


The equation of 6, and force in elastic-plastic flow region show increased trend 
and could be calculated as: 


» [ à, mõ tag | 
à(a)- 4 : uu (15) 
Gg, uU c(l-sing,) 


Fa. i (ea) a, tan cosa: da (16) 
Fo (Ay) -[ ó. CORN tan? @-de (17) 
Force equation of cutting grain could be obtained by the above equation: 


Fa, J= Fay Gn) Lu le, d+ Fon (d, ) 


à; S. (18) 
- (^5 (aa, tang- cosa da- |? 8 «a, tan -oosa da + Ö-A, 
5 (a, j= Fa Aa, ) Foron ) Faxon (a, ) 
T(L-a i] (19) 


= [, 8 (e) de - 2. 6, tane - 06, 4, fan @ 
o 2-(1 -sin c. ) 2 g 


3.4 Force model of cutting and ploughing grain 


Material behavior in ploughing stage could be categorized as plastic flow, in 
which stress should reach at 6,. Therefore, force model of ploughing grain could be 
calculated as Eqs. (20) and (21). 


F,(a,)7 f 8, -a -tan@-cosa-de (20) 
$T : 
F(a,)= f: 6,-a,° -tan Oda = - 8,2, -tan’ 8 (21) 


3.5 Frictional force model 


The frictional force is composed of rake face on ploughing grains Fy ¢(ag), rake 
face on cutting grains F.¢(ag), and the wear plane of the grain (f, ft), as shown in 
Eqs. (22)-(30). The force magnitude is determined by stress and the lubricating 
state between the grains and the workpiece. 


= x i x > Lk 0 
Fula - 2: |? u-8, ds = fius, TUE «he (22) 
SE w " ^o pns 
fda, aL Cer, sind = (: u-Ó,-d ^ tan? 8- da (23) 
Fla, E eia, cos = E 108, -a tan Ode (24) 


F(a.) - 2-(ue 3, ds È" wd, ds) 


> land 


Ll > P T 5 
= f a 0, - 5^ edu: -du 
et . D ^ 


CUN 
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anf ; lung 


-da-[" pð a — da (26) 
a ^ c 


; Z " 2b 
f(a, ) =f, (a, ):sinO = [ 4-8, VÀ 2 
1 Uus 


- 


Piglet, 7 Fy (a, )- vos Ü = [246 a” - tan da [ Hà duy ‘tan@-da (27) 


— 4.B.S,-N,-V,  4-K,-N,-V, 


f, -N,S, 2 

eee ee V.D VD (28) 

: — 4-u- P. 6 NE, 4-u- K,- N,-F, 

f = ` N MS S pz D € a w D 1 Pd a 

i a V-D V.D (29) 
K,=P,-S,, (30) 


where K; is the physical quantity related to the grain’s wear state of the grinding 
wheel, which could be reversely solved through grinding experiment, and y is the 
friction coefficient between the workpiece and the grains, which could be obtained 
by tribological tests [40]. 


4. Grinding force model and prediction 
4.1 Procedure of modeling common grinding wheel 


The matrix of the location distribution of the grains in the grinding zone is 
established as G(z,). The matrix of the protrusion height of the grains can be 
calculated as [40]: 


4.2 Dynamic active grains in grinding zone 


For two continuous cutting grains, maximum undeformed chip thickness dgmax 
(ag = 0.545... [42]) could be solved as [41, 42]: 


CE Ru LAT (32) 


where 4 is the space between the continuous cutting grains, D is the diameter of 
the grinding wheel, Vw is the feed speed, and V, is the peripheral speed of the 
grinding wheel. 

For common grinding wheel, 4 value represented the space between the 
dynamic active grains. The equation of the nth grain is deduced as: 


d 


K , 
O, wt) = 2A. 1) YU B t (5o) 7 Cy n) (33) 


where @gmax(n) is the maximum undeformed chip thickness of the nth dynamic 
active grain, 4(5..5..1 is the space between the nth and the (n — 1)th dynamic active 
grains, 4p(n) is the protrusion height of the nth dynamic active grain, and a; (n — 1) 
is the protrusion height of the (n — 1)th dynamic active grain. 
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4.3 Grinding force model 


Concrete calculation of grinding force can be described in Figure 13; the grind- 
ing force can be expressed in Eqs. (34) and (35). 


F, = XF, (2,) 1 Fg (4,9): 315, Gand! Fy mt f, 84 


F, 2Y LE (aE y 0,)]: Y IE 5, Gg) f, 85) 
| | 


where F/F, are the tangential/normal grinding force, respectively, and f,/f, are 
the tangential/normal frictional force on the wear plane of grains, respectively. For 
the nth cutting grains (1 <n < Nj), ag, is the cutting depth; F,-(@gn)/Fnc(4gn) are the 
tangential/normal cutting force, respectively; and Fy¢(4gn)/Fnce(4gn) are the tan- 
gential/normal frictional force on the rake face of the cutting grain, respectively. 
For the mth ploughing grains (1 < n < Ne, Ne + Np = Na), 45, is the cutting depth; 
Ftp(agm)/Fnp(agm) are the tangential/normal, respectively; and Fip¢(@gm)/ Fnpt(4gm) 
are the tangential/normal frictional force on the rake face of ploughing grain, 
respectively. 


4.4 Experimental verification 


MoS,-palm oil nanofluid minimum quantity lubrication grinding experiment 
was carried out with K-P36 surface grinder. Prediction and experimental force 
values was shown in Figure 14. 

The average percentage of the deviation in the normal force is 4.1996, while in 
the tangential force is 4.31%. For a certain grinding condition (dry grinding, V, = 20 
m/s, Vw = 2 m/min, a, = 15 um), the contribution of the tangential frictional force in 
the total tangential grinding force is approximately 89.1796 and approximately 
90.7196 for normal force. However, it was decreased to 86.5296 and 89.4396 for 
tangential and normal direction, respectively. 
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Figure 13. 
Process of the grinding force predictive program. 
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Figure 14. 
Force prediction and experiment result. 


5. Conclusions 


To develop a model based on the stress state of a single grain, this study 
investigated the material removal mechanism and interference action in the 
grinding zone. The combined effect of the material removal and plastic stacking 
of a single grain has been considered. 

Cutting experiment of single grain was carried out to obtain the critical cutting 
depth of the cutting and ploughing states. Results indicated that the variation trend 
of the grain-cutting efficiency (f) presents an S-shaped trend as the cutting depth 
(ag) increases. 

This study also developed and verified an improved grinding force predictive 
model. The average percentage of the deviation in the normal force is 4.19%, 
whereas that in the tangential force is 4.31%. 
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Chapter 7 


Kerosene-Water Multiphase Flow 
in Vertical and Inclined Pipes 


Faik Hamad, Nadeem Ahmed Sheikh and Muzaffar Ali 


Abstract 


This chapter presents the volume fraction distribution of kerosene-water two- 
phase flow in vertical and inclined pipes. The study of liquid-liquid two-phase flow 
is very significant to oil industry and many other processes in industry where two 
liquids are mixed and flow together. Pitot tube and optical probes are used for the 
measurement of velocity of water and volume fraction. The experimental measure- 
ments of the local parameters demonstrate that the single-phase and two-phase 
flows reached the fully developed axisymmetric conditions at L/D > 48 (L, pipe 
length; D, pipe diameter). The results also showed the severe asymmetry distribu- 
tions of the volume fraction at the entrance region (L/D = 1) downstream the bend 
and in the inclined pipe. The comparison of volume fraction profiles with void 
fraction profiles indicated a significant difference in their shapes. The results also 
showed that the kerosene accumulated at the upper wall of the inclined pipe and the 
distribution improved by increasing the volumetric quality. 


Keywords: volume fraction, kerosene-water two-phase flow, vertical and inclined 
pipes, optical probe, pitot tube 


1. Introduction 


Multiphase flows are important for the design of steam/water flow in steam 
generators, jet engines, condensers, extraction and distillation processes, gas and 
oil mixture in pipelines, and refrigeration systems. The mixture of two immiscible 
liquids is characterized by the existence of interfaces between the two fluids, associ- 
ated with a discontinuity of properties across the interface. The single-phase flow 
is traditionally classified into laminar, transitional, and turbulent flows according 
to the flow Reynolds number. The two-phase flow in vertical pipe can be classified, 
according to the geometry of the interfaces. 

The primary condition for all two-phase flows is specified by the volumetric 
quality £, which is defined as: 


2 Q4 
Bots o 


where Q, is the flow rate of the dispersed phase and Q, is the flow rate of the 
continuous phase. 

For a pipe of radius R, the corresponding (area averaged) superficial velocities 
are defined as: 
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Continuous phase superficial velocity: 
U, -2 (2) 
Dispersed phase superficial velocity: 


Uu (3) 


zR? 

For both gas-liquid and liquid-liquid flow systems, the continuous phase is 
usually water. 

In spite of the large number of published work in multiphase flow area, the 
publication on using local probe measurements for liquid-liquid flow is very 
limited compared to the gas-liquid two-phase flow. The purpose of this chapter is to 
publish some data on volume fraction profiles for liquid-liquid flow in vertical and 
inclined pipes. The following data are presented in this chapter: (i) the void fraction 
distribution for gas-liquid two-phase flow in vertical pipe, (ii) the volume fraction 
distribution for flow development of kerosene-water flow in vertical pipe, and (iii) 
the volume fraction distribution for the fully developed kerosene-water flow in 
vertical and inclined pipes. 


2. Void fraction/volume fraction definition 


Most experimental results for the void fraction a (volume fraction for liquid- 
liquid flow) have been obtained by a point sensor, which was used to discriminate in 
time between the two phases. Experimentally the void fraction a has been evaluated 
from the time record from such a probe as: 


_ LAty 


a= 24 (4) 


where > At, is the time the probe is located in the dispersed phase and T is the 
total sampling time used for record. It should be noticed that in many investiga- 
tions, reference is made to the (average) void fraction a. The proper reference 
would have been the volumetric quality (f) obtained for gas-liquid systems by the 
quick closing valve method, X-ray or neutron techniques. For liquid-liquid systems, 
f can be obtained by measuring the two flow rates Q, and Q.. 

Experimental studies of the phase distributions in concurrent two-phase upflow 
in vertical pipes present a complex picture that has not yet been systematically 
evaluated. The common flow patterns for vertical upward flow, in which both 
phases flow upwards in a circular tube, are shown in Figure 1. As the volume flow 
rate of gas increased for constant water flow rate, the flow patterns would vary. The 
following types of flow patterns can be found in vertical pipes: 


i. Bubbly flow: bubbles of gas or liquid in a continuous liquid phase appear, 
and the size of the bubbles can be very small or large. 


ii. Slug flow: in this type a bullet-shaped plug of gas is formed from many 
bubbles concentrated in one part to make larger bubbles, which approach the 
diameter of the pipe. The liquid phase is in continuous flow. 


iii. Churn flow: braking down of large vapor bubbles in plug flows form the 
churn flow. This is a highly oscillatory flow, and there is tendency for each 


phase to be continuous with irregular interfaces. 
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Figure1. 
Flow patterns in vertical upward flow in a pipe [1]. 


iv. Annular flow: the liquid forms a film around the wall of the tube. The gas 
phase flows in the centre. 


3. Gas: liquid void fraction distribution in vertical pipes 


A significant number of measurements have been made for upflow in vertical 
pipes. Several investigators, e.g. Malnes [2], Serizawa et al. [3], Michiyoshi and 
Serizawa [4], Wang et al. [5] and Liu and Bankoff [6], have observed the peaking 
phenomenon of the local void fraction near the wall as shown in Figure 2. Some 
investigators, such as Van der Welle [7], Moujaes and Dougall [8] and Johnson and 
White [9], have observed a maximum void fraction at the centreline as shown in 
Figure 3. Other researchers, such as Nakoryakov et al. [10], Spindler et al. [11] and 
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Figure 2. 


Void fraction distributions for air-water upflow in a pipe with a continuous phase mean velocity = 1.03 m/s, 
Z/D = 30. A-D are bubbly flow, and E and F are slug flow [3]. 
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Figure 3. 
Void fraction distributions for air-water upflow in a pipe with different flow rates, Z/D = 42.5. A-D are 
bubbly flow, and E and F are transitional flow [9]. 


Liu [12], observed both wall and centreline peaking void fraction distributions 
for two-phase flow. The actual void fraction distribution configurations have been 
found to depend on the initial conditions: bubble size and flow rates, physical 
properties of the fluids, and the test section condition geometry. 


4. Liquid: liquid volume fraction 


Compared to the large number of publications on gas-liquid flows, less work 
have been published on liquid-liquid flows. 

Most of the papers on liquid-liquid mixture flow were published by research 
group at the University of Bradford ([13, 14]; Hamad et al. [15]; Hamad and Bruun 
[16]). Most of these papers focused on the development of optical techniques 
for kerosene-water upward flow in vertical pipes. However, Farrar and Bruun 
[13] highlighted the problem of the severe asymmetry, and the swirl generated 
upstream the inlet due the existence of 90° bend as part of the experimental 
facility. 

Zhao et al. [17] used a double-sensor conductivity probe to measure the local oil 
phase fraction distribution for flow in a vertical pipe at L/D - 72. They found that 
the volume fraction profiles were uniform for f < 9.2% and changed into wall peak 
for f > 9.2%. The local oil phase fraction profiles at (a) constant water flow rate 
(Jw = 0.33 m/s) and (b) constant oil flow rate (J, = 0.066 m/s) are given in Figure 4. 

A comprehensive experimental data on kerosene-water two-phase flow were 
published by Hamad et al. [18, 19] and Hamad et al. [20] in vertical and inclined 
pipes. A summary of the results from each paper is given in the following sections. 


4.1 Development of kerosene-water flow in vertical pipe 


Hamad et al. [18] studied the flow development in a vertical pipe of 77.8 mm 
inner diameter and 4500 mm length downstream of a 90? bend experimentally 
at L/D - 1, 16, 38 and 54 using the experimental facility in Figure 5. Single-phase 
(water) flow measurements were made to check the establishment of fully devel- 
oped symmetrical flow conditions. Figure 6 shows the radial distribution of axial 
velocity in the plane parallel to the bend at different L/D ratios. Two values of U,, 
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a 06 


Figure 4. 
Oil volume fraction profiles at (a) constant water flow rate (Jy = 0.33 m/s) and (b) constant oil flow rate 
(Ja = 0.066 m/s) [17]. 


Figure 5. 
Schematic diagram of two-phase flow experimental facility [18]. 


are used (0.44 m/s (Re = 33,800) and 0.77 m/s (Re = 60,000)). The results show 
that water velocity distributions become fully developed at L/D > 48. The empirical 
power law velocity distribution given in Eq. (5) for single-phase turbulent flow 

[21, 22] was also included in Figure 6 to confirm the accuracy of the measurements: 


U -Ua -r/R)"" (5) 


Then, the kerosene was introduced to perform volume fraction measurements 
using optical probe [14] at four different axial positions at L/D = 1, 16, 38 and 54 
downstream of the pipe bend. Three different flow conditions are considered: 
Case 1: water superficial velocity, U,, = 0.44 m/s, and volumetric quality, f = 9.2%; 
Case 2: U,, = 0.44 m/s and f = 18.696 and one high v,, condition; and Case 3: U,, = 
0.77 m/s and f = 18.6%. For Case 1, the axisymmetric distribution is very poor at 
L/D = 1 (Figure7(a)) with high volume fraction values of 20% near the inner side 
of the bend to the lower value near the outer wall of the bend of 496. For Case 2, 
increasing f to 18.6% for the same U, (Figure 7(a)) improves the axisymmetric 
distribution across the pipe. For Case 3, the axisymmetric distribution improved 
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Figure 6. 


Local single-phase velocity distribution at different L/D ratios for average water velocity = 0.44 m/s (solid 
symbols) and 0.77 m/s (open symbols). The power law velocity distribution (solid line) is also included [18]. 
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Volume fraction distributions at different (a) L/D for U „ = 0.44 m/s and f = 9.2% (solid symbols) and 
18.6% (open symbols) [18], (b) L/D for U „ = 0.77 m/s and f = 18.6% [18]. 
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further by increasing U, for the same f (Figure 7(b)). The main conclusion is 
that the axisymmetric becomes better for higher f and v,,. The change in volume 
fraction distribution may be attributed to the improvement of the mixing process 
of the kerosene with water which acts against the effect of the buoyancy force and 
the centrifugal force at the outlet of the bend. However, the distribution becomes 
nearly symmetrical at L/D = 16. There appear to be no significant differences 
between the distributions at L/D = 38 and 54, which suggests that fully developed, 
symmetrical condition was achieved. 


4.2 Fully developed flow of kerosene-water flow in vertical pipe 


Hamad et al. [19] studied the flow of kerosene-water upward flow in a vertical 
pipe at (L/D = 54) using optical probes. The effects of 7,,, and f on radial volume 
fraction distribution [a(r)] of two-phase flow parameters were investigated. 

The local volume fraction is calculated from the output of the leading sensor 
of the dual optical probe by determining the average drop residence time using 
the procedure described in Hamad et al. [14]. Comprehensive measurements 
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Figure 8. 


Volume fraction profiles for different (a) U „and B = 4.6, 9.2 and 18.6% [19], (b) U „and f = 28.2, 38 and 
47% [19]. 
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Figure 9. 
The effect of pipe inclination on volume fraction distribution (a) (Uy = 0.29 m/s, B = 9.296) [20], (b) 
(Uw = 0.29 m/s, P = 18.696) [20]. (c) the o (x) distribution from Vigneaux et al. [24]. 


of a(r) were performed for a number of f values in the range of 4.6-4796 and 
constant U,, of 0.29, 0.44, 0.59, 0.69 and 0.77 m/s. The a(r) profiles have 
been plotted together for various values of U,,, for each value of p as shown in 
Figure 8(a) and (b). 
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As the a(r) profile primarily reflects the kerosene content in the mixture flow, it 
follows that the related a(r) profile sets for f = 4.6%, 9.2%, 18.6%, 28.2%, 38% and 
47% are centred around these values. The graphs also show distinct variations, both 
within each f group and between groups with different / values. 

The results from Figure 8(a) and (b) show that increasing U, with low f 
(<20%) will change the a(r) profiles from convex shape with peak at the pipe 
centreline to flat shape and then to concave shape with peak near the wall. For 
moderate f (20-40%), the a (v) profiles have a concave shape for different U,, with 
peak near the wall which has high values for higher U,,. In the case of f ~ 50%, the 
a(r) profile shapes are flat for the two cases in Figure 8(b). 

The a(r) profiles from centreline which peaked to uniform to wall peaked and 
then to uniform can be attributed to the change in lift force due to the change in 
drop diameter, slip velocity and radial velocity distribution of both phases. The 
present finding is supported by the results for liquid-liquid flows from Zhao et al. 
[17] and Hua et al. [23] for the same range of U, and f. 


4.3 Kerosene-water flow in inclined pipe 


Hamad et al. [20] used an optical probe to study the kerosene-water flow 
inclined at 5? and 30? from vertical at L/D - 54. The volume fraction was measured 
for U,, = 0.29 m/s and f = 9.2% and 18.6%. 

Figure 9(a) shows the radial a) distributions of the volume fraction, a(r) for 0°, 
5? and 30° inclination angles at U,, of 0.29 m/s and two values of p = 9.2 and 18.6%. 

The results in Figure 9(a) and (b) show that the inclination has a significant 
influence on the distribution of a(r). The kerosene drops were separated from the 
water accumulated at the upper zone of the pipe due to the gravity effect. The effect 
of increasing f in an inclined pipe leads to dispersion of the drops to the lower zone 
of the pipe due to the recirculation cells of the moving droplet swarms. 

The present results are supported by the findings reported by Vigneaux et al. [24] 
and Flores et al. [25]. Figure 9(c) presents the two sets of experimental data reported 
by Vigneaux et al. [24] ina pipe inclined at 15? from vertical. In the first case, 

p = 2396, and U,w = 0.27 m/s, and in the second case, p = 40%, and Usw = 0.21 m/s. 


5. Conclusion 


The results on void fraction profiles from literature show the complexity of the 
flow behaviour. It is reflected in different types of profiles due to the local interac- 
tion between the bubbles and the continuous phase. This may be attributed to the 
various forces at interface between the phases including drag, lift and virtual force 
as well as the size of bubbles and compressibility effect. In contrast, the volume 
fraction profiles for liquid-liquid two-phase flow have similar shapes. This behav- 
iour may be attributed to smaller drops, smaller density ratio, smaller slip velocity 
and the incompressible nature of the liquids. 

The results show that fully developed condition for liquid-liquid flow can be 
achieved at lower L/D compared to gas-liquid flow. This is due to incompressible 
nature of liquid drops which have the same volume compared to the gas bubbles 
which expand continuously due to the pressure drop in flow direction. 
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Abstract 


Polydimethylsiloxane (PDMS), due to its remarkable properties such as optical 
transparency and ability to easily mold, is one of the most popular polymers used in 
micro- and nanofluidics. Furthermore, 3D printing technology due to its low cost 
and simplicity is also gaining a great interest among the microfluidic community. 

In this work, the potential of 3D printing is shown to produce microfluidic devices, 
their ability for studying flows and heat transfer of nanofluids, and their applica- 
bility as a heat sink device. The low-cost fused deposition modeling 3D printing 
technique was combined with a PDMS casting technique for the microfluidic device 
fabrication. The potential of this technique was experimentally demonstrated by 
fluid flow and heat transfer investigations using different fluids, such as distilled 
water-, alumina (Al,O3)-, and iron oxide (Fe3O4)-based nanofluids. The simplic- 
ity, low-cost, and unique features of the proposed heat sink device may provide a 
promising way to investigate nanofluids' flow and heat transfer phenomena that are 
not possible to be studied by the current traditional systems. 


Keywords: microfluidics, 3D printing, microchannel heat sinks, nanofluids, 
heat transfer, electronics cooling 


1. Introduction 


The continuous investigation on strategies for size reduction while maintaining, 
or even increasing, power of technological devices demands cooling systems with 
higher thermal efficiency and smaller sizes. One approach relied on the modifica- 
tion of heat sinks by incorporating microchannels to increase the heat exchange 
surface. Despite being a well-adopted strategy, its complex configurations proved 
to be difficult to manufacture [1, 2]. Other strategies focused on the type of the 
used fluids for the cooling process. Dielectric fluids [1, 3], two-phase fluids [3, 4], 
and nanofluids (NFs) [5-8] have been thoroughly investigated. Among these three, 
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the nanofluids, which are fluids comprised of particles, with size ranging from 1 
to 100 nm, suspended in a base fluid, have been reported as presenting a better 
thermal conductivity than the base fluid [8, 9]. Nevertheless, some challenges 
regarding their usage, such as agglomeration, long-term stability, and high-costs, 
still need further investigation [8, 10-12]. The most commonly used nanoparticles 
(NPs) for NFs are metallic, such as Cu, Ag, Au, and Fe, or non-metallic such as 
ALbO;, CuO, TiO; SiC, and carbon nanotubes. Many authors reported that the heat 
transfer of NFs is influenced by several factors such as shape, dimensions, volume 
fractions in the suspensions, and the thermal properties of the particle materials 
[11, 13, 14]. Therefore, enhancement in heat transfer was also reported, but only for 
small concentrations of NPs [8, 9]. 

Abareshi et al. [15] have evaluated the thermal conductivity of nanofluids with 
different volume fractions of Fe3O, NPs, at different temperatures. The thermal 
conductivity was reported to increase up to 11.596 for a volume fraction of 3 vol96 
at 40°C. Xia et al. [16] have investigated the heat transfer coefficient of nanofluids 
using TiO; and Al,O; NPs, with different volume fractions. For a volume fraction 
of 196, the heat transfer coefficient was significantly increased for both nanofluids, 
compared with deionized water. Gavili et al. [17] have studied the thermal con- 
ductivity of ferrofluids with FeO, particles of approximately 10 nm in diameter, 
suspended in deionized water. With the application of a magnetic field, the thermal 
conductivity was increased up to 200% for a 596 volume fraction. Kim et al. [18] 
have investigated the thermal conductivity of alumina- and distilled water-based 
nanofluids, with concentrations of 0.5, 1, and 2 wt%. The conductivity was found 
to increase with the increasing of the NPs concentration. Al-Rjoub et al. [19] have 
tested four cooling liquids: deionized water; distilled water; borax buffer; and AO; 
NPs solution, on a microscale heat exchanger. It was found that the deionized water 
has presented the lowest heat removal capacity, while the Al,O; solution showed the 
highest capacity, corresponding to about 6996 increase. 

The majority of microchannel heat sink devices that can be found in literature were 
fabricated in silicon, due to its thermal conductivity. However, the fabrication process 
of those devices can be laborious and needs extremely expensive facilities. Novel, 
fast, and low-cost fabrication techniques have been developed by means of different 
kinds of polymers [20-22]. PDMS isa silicone elastomer with a set of properties that 
make it suitable for many applications and is a popular choice for microfluidic devices 
fabrication [21, 22]. Besides being cheaper than the monocrystalline silicon, it presents 
alow elasticity change versus temperature, high thermal stability, chemical inertness, 
dielectric stability, shear stability, high compressibility, and hyperelasticity [23-27]. 
Moreover, it is non-toxic and biocompatible [25, 27-29]. PDMS devices can be manu- 
factured by simple techniques at room temperature, such as replica molding. 

The 3D printers are gaining an increased attention by both academic and 
industrial community to produce microdevices and models at an extremely low 
cost. Some successful applications can already be found in lab-on-a-chip tools [30], 
microfluidics [31, 32], and biomedical in vitro devices [27, 32-34]. There are differ- 
ent kinds of printing methods such as the Fused Deposition Modeling (FDM) and 
stereolithography [32]. Among those methods, the most popular, simple method 
with the lowest cost is the FDM technology [32]. For this reason, the FDM process 
was combined with a PDMS casting technique to produce a PDMS microfluidic 
device. Hence, the main objective of this work is to show the potential of a FDM 3D 
printer to produce microfluidic devices and their potential to be used to perform 
flows and heat transfer studies of nanofluids. To demonstrate the potential of this 
technique, fluid flow and heat transfer studies were performed by using different 
fluids such as distilled water-, alumina (Al, O3)-, and iron oxide (Fe404)-based 
nanofluids with concentrations of 1 and 2.596. 
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2. Materials and methods 


This section describes the experimental protocol to develop the nanofluids and the 
fabrication process to produce the heat sinks, as well as the used experimental setup. 


2.1 Nanoparticles preparation 


Two different types of NPs were used on the flow and heat experiments: alumina 
oxide (Al; O4) NPs acquired from Sigma Aldrich (ref. 702,129, Sigma Aldrich) and 
iron oxide (Fe3O4) NPs synthesized on our laboratory by co-precipitation. This last 
method allows to produce magnetic iron oxide NPs in a cost-effective way and is 
appropriate for mass production. The AL; O4 NPs had a size inferior to 50 nm, while 
the synthesized Fe; O4 NPs had an average size of 11 + 2 nm. 

The co-precipitation was initiated with the preparation of the precipitation 
agent by adding 0.01 g of cetrimonium bromide (CTAB), diluted in 3 mL of 
distilled water, to 20 mL of ammonium hydroxide (N H4OH). Hereinafter, a ferrous 
solution was prepared by diluting 7.78 g of iron(II) chloride (Fe Cl) and 4.06 g 
of iron(II) chloride (Fe CL) in 20 mL of distilled water, in an ultrasonic bath. The 
solution was subsequently mechanically stirred at, approximately, 1500 rpm. The 
precipitation solution was then added, dropwise, to the ferrous solution under 
stirring, on a laminar flow cabinet. The co-precipitation occurred according to the 
following equation: 


Fe CL.4 H, Ow + 2Fe Cl46 Hy Ow + 8NH, (OH) (aq) 
= Fe3O45 t 8N H4 Cla, + 20 H,O0 (1) 


To conclude the process, the NPs were washed several times with distilled water 
with the assistance of a strong magnet. 

Figure 1 shows the synthesized NPs and the representative transmission elec- 
tron microscopy (TEM) image. A detailed description of this process can be found 
at Cardoso et al. [35]. 

Despite the TEM images show the aggregates of Fe;0, NPs, in the NF, they were 
stable and non-aggregated for a period of time after 10 min of ultrasonic bath. 

The NPs show a normalized magnetization of ~69 emu.g 1 at ~10 kOe, which 
corresponds to the saturation magnetization, and also show a superparamagnetic 
behavior with an extremely low coercivity of 1.6 Oe [35-37]. 


Magnetic iron oxide 
NPs synthesized by 
co-precipitation method 


E * 3 
e—a 50 nm 


Figure 1. 
Magnetic iron oxide NPs and representative TEM image, adapted from [35, 36] . 
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Figure 2. 
Schematic representation of the main dimensions of the ABS master mold. 


2.2 Fabrication of the heat sink microchannel 


The heat sink microchannel device was produced based on a scaffold-removal 
technique [25]. First, the molds were drawn by using the Autodesk Inventor® 
software and then printed at the FDM 3D printer Ultimaker 2+ (Ultimaker, 
Netherlands). The first mold was printed with acrylonitrile butadiene styrene 
(ABS), whereas the second one was printed with polylactic acid (PLA). The fabrica- 
tion of the molds was performed with a nozzle with a diameter of 0.4 mm, whereas 
the layer resolution was about 100 um. The main dimensions of the ABS master 
mold can be found in Figure 2. 

Once the 3D models were printed, PDMS was prepared by adding a PDMS cur- 
ing agent into the pre-polymer with a mixing ratio of 1:10. The PDMS was poured 
onto the PLA mold with the ABS master mold inside it. Once the PLA mold was 
filled with PDMS, the elastomer was cured at room temperature for about 1 day. 
Finally, the PDMS was removed from the PLA mold and immersed in an acetone 
bath to remove the ABS for approximately 24 h. Figure 3 shows the schematic dia- 
gram of all the main steps to produce the PDMS heat sink device. The overall cost to 
fabricate the PDMS heat sink device is about 3.8 €. This cost includes the printing of 
the ABS master mold (~1 €) and PDMS casting process (~2.8 €). 

Notice that after the PDMS curing process, small holes were made below the 
inlet and outlet to insert the thermocouples (type K). Figure 4 shows the PDMS 
heat sink microfluidic device used in the flow and heat experiments. 


2.3 Experimental procedures 


The PDMS heat sink was placed on top of a hot plate controlled by a 9400- 
temperature controller (CAL Controls). The temperature of the plate was set 
to 60°C, whereas the flow rate of the fluids was controlled by a syringe pump 
(Harvard) connected to the inlet of the heat sink. The temperature at the entrance 
and exit of the device was acquired through a data acquisition instrument con- 
nected to the thermocouples of the device. Wood and polystyrene blocks were used 
to minimize the heat losses. Figure 5 shows a schematic diagram of the experi- 
mental setup. The flow of the Fe4O, nanofluid was analyzed by optical microscopy 
ata flow rate of 10 uL/min. Note that the temperature measurements presented 
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Figure 3. 
Schematic representation of main steps to fabricate the PDMS heat sink. 
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Figure 4. 
PDMS heat sink microfluidic device with the inserted thermocouples. 


an uncertainty of +1°C. The thermographic studies were performed with distilled 
water at a flow rate of 1, 5, 7.5, 10, and 15 mL/min and a thermographic camera 
Onca-MWIR-InSb (Xenics Infrared Solutions). The setup and calibration proce- 
dures of the camera were performed as in Teodori et al. [38]. Images, with a resolu- 
tion of 150px x 150px, were taken with a frame rate of 1000 fps. 

To evaluate the influence of the nanofluids properties in the heat sink microfluid 
device, the tests were performed using distilled water, Fe O4 at a concentration of 1 
and 2.5%, and Al; O; at the same concentrations. All the fluids were set to a flow rate 
of 1 up to 30 mL/min. 


2.4 Heat transfer calculations 


The properties of the nanofluids were obtained taken into account fundamental 
equations described on previous studies [39, 40]. The thermal conductivity of the 
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Figure 5. 
Schematic representation of the experimental setup. 


nanofluid was obtained according to the Maxwell model described by the following 
Equation [39, 40]: 


K, + 2Ky + 29(Ky - =) 


Ey = Koy (z + 2Kiy - 29(K, - Ky) 2 


where K,; is the nanofluid thermal conductivity, K, is the NPs thermal conduc- 
tivity, Ky is the base fluid thermal conductivity, and q is the NPs concentration. 

The nanofluid density and heat capacity were calculated through the weighted 
average of the individual properties of both the NPs and base fluid. The first is 
expressed by Eq. (3) [14, 41] and the latter by Eq. (4) [11, 14]. 


Prf = Pp® + por(1 - p), (3) 


CD np = P ppCp,, T (1 E p) Pfo Cp (4) 


In the abovementioned equations, p,¢ represents the nanofluid density, p, the 
particle density, py the base fluid density, and c,,, and c,,, the specific heat capacity 
of the NPs and of the base fluid, respectively. 

The nanofluid viscosity was determined through the equation proposed by 
Batchelor [42]: 


Ung = (1+ 259 + 629?) up, 0 < q < 1096 [5] 


where Lis the nanofluid viscosity and uj the base fluid viscosity. This equation 
brings in a quadratic dependence with the volume fraction, which provides a better 
representation of the interaction between the particles on the fluid. 

Within the microfluidic device, the heat transfer will occur by convection inside 
the microchannels and by conduction in the walls between them. Consequently, the 
mathematical approach that better allowed the evaluation of the heat transfer was 
described by Ma et al. [10]. The convection heat transfer coefficient was calculated 
by iterations using the following equations: 


: Q 
h D N(A, + 2nA)(T, = Tavgf) ` (5) 
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2hi4 
Mm; = — 6 
] | RppmsWp ( ) 
th (miz, 
i= aes, (7) 
where A, represents the area of microchannel bottom, A; the area of microchan- 
nel sidewall, Tif the average temperature of the fluid, n the fin efficiency, w, the 
average width of fin, and z, the height of the channel. The main dimensions of 
the heat sink are illustrated in Figure 2. The heat transfer rate, Q, represents the 
amount of heat energy taken by the fluid when it flows through the channels and is 
given by Eq. (9). This parameter was obtained for each mass flow rate, m, after the 
temperature at the inlet, Tn, and at the outlet, Tou, were measured. 


Q -q PC of Tout = Tin)» (8) 


where q, cp, and pris the volume flow rate, specific heat capacity, and density of 
the working fluid, respectively. 


3. Results and discussion 


In this section, the obtained results are presented and discussed. The tempera- 
ture measurements and the known properties of the materials were used to calcu- 
late the parameters described on the previous section. Using those parameters, the 
influence of the environment conditions and of the fluid properties in heat transfer 
was analyzed. 


3.1 Influence of the nanofluid properties 


As described previously, the prepared nanofluids with iron oxide (Fe;0,) and 
alumina (Al, O4) NPs were used to verify the thermal properties influence of the 
nanofluids in the heat transfer performance of the developed PDMS heat sink device. 

Figure 6 shows the temperature difference between the inlet and the outlet 
as the tested nanofluids flow through the proposed PDMS heat sink. Overall, it is 
possible to conclude that both tested nanofluids present a bigger temperature dif- 
ference between the inlet and the outlet in comparison with distilled water. Hence, 
these measurements indicate that the amount of heat energy absorbed by both 
nanofluids was bigger than that absorbed by the distilled water. In addition, the 
amount of heat absorbed by the nanofluids was found to be bigger for smaller flow 
rates, as shown in Figure 6. These results corroborate the measurements performed 
by Chein and Chuang [43], where they have investigated the heat performance of 
nanofluids with CuO NPs in a microchannel heat sink. These results also show that 
the flow rate affects the amount of heat absorbed by the nanofluids. 

In Figure 7, an increase of the convective heat transfer coefficient was regis- 
tered on both nanofluids, in comparison to the base fluid. The increase was more 
pronounced for the alumina nanofluid due to the greater stability of the nanofluid 
and bigger thermal conductivity of these kind of NPs [6]. By increasing the convec- 
tive heat transfer coefficient, an increase of the heat transfer rate is expected. In 
Figure 8, it is possible to observe an increase in the heat transfer rate for both tested 
nanofluids in comparison with distilled water. From these results, it is also possible 
to conclude that for both nanofluids, the convective heat transfer coefficient and 
heat transfer rate increase with the flow rate, which agrees with the results obtained 
by Wen and Ding [44]. 
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Temperature difference between the inlet and the outlet for the tested nanofluids at the proposed PDMS heat 
sink microfluidic device. 


- a Distilled water 2 
o 
Fe2Oa [1%] P 
704 4 AbOs [196] P. 


0 2 E 6 8 10 12 14 16 
Flow Rate (mL/min) 


Figure 7. 
Convective heat transfer coefficient of distilled water, alumina, and iron oxide nanofluids as the function of the 
flow rate. 


Figure 9 compares the convective heat transfer coefficient for two different con- 
centrations of NPs, i.e., 1 and 2.596 of both Fe3O4 and Al; Os. From a macroscopic 
view, it was noted a better dispersion for the nanofluids containing 196 of NPs. In 
addition, by increasing the concentration, the heat transfer was not enhanced. In 
fact, it was noted a decrease of the convective heat transfer coefficient when the 
concentration of NPs was increased to 2.596. Although these results are somewhat 
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Figure 8. 
Heat transfer vate of distilled water, alumina, and iron oxide nanofluids as the function of the flow rate. 
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Figure 9. 
Convective heat transfer coefficient of nanofluids as the function of the flow vate, for two different 
concentrations of alumina and iron oxide NPs. 


paradoxical, other researchers, such as Wen and Ding [45] and Putra et al. [46], 
have reported similar results. The main possible reasons for the seen heat transfer 
deterioration include both the aggregation and sedimentation of the NPs. However, 
the possible reasons and mechanism attributed to such phenomena require further 
research. Currently, our group is carrying out both experimental and numerical 
work to identify the exact causes for such phenomena. 


3.2 Optical and thermal imaging analyses 


The biggest advantage of the developed PDMS heat sink device is the ability to 
visualize the flow phenomena happening inside the microchannels. By using a high- 
speed video microscopy system, it was possible to visualize several flow phenomena 
of the nanofluids such as the formation, growing, and breakdown of NPs clusters 
(see Figure 10). From these observations, it was concluded that one of the main 
causes for the formation of the clusters was the high roughness of the PDMS surface 
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Figure 10. 
Optical images of the formation, growing, and breakdown of a cluster of NPs. 


Figure 11. 
Optical image of the surface roughness of the heat sink device used in this study. 


channels (Figure 11). This verified roughness was caused by the ABS master mold 
fabricated by the FDM 3D printer. In order to improve the surface roughness, the 
ABS master molds should undergo an acetone vapor treatment before performing 
the PDMS casting procedure. More detailed information about this method can be 
found elsewhere [27]. 

Another interesting advantage of this PDMS microfluidic device is the ability to 
visualize both the flow and thermal performance of the system by using a thermo- 
graphic camera, as shown in Figure 12. Notice that the temperatures acquired were 
from the surface of the heat sink device and not directly from the working fluid 
flowing in the microchannels. In fact, PDMS is transparent to visible radiation but 
partially opaque to the infrared (IR) radiation. Hence, to obtain the temperatures 
more closely related to the working fluids flowing through the microchannels, 
the thickness of the upper walls should be reduced in future experiments. A very 
interesting observation was the ability to detect bubbles that are likely to happen 
in microfluidic devices. In Figure 13, it is possible to visualize a bubble within the 
microchannel and the thermal performance of the heat sink device. 


3.3 Limitations and future directions 
In this study, a microchannel device was successfully manufactured and used 


in microfluidic essays. Nevertheless, some limitations arose throughout the work. 
Despite the advantages of using PDMS, some properties of the material, such as 
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Figure 12. 
Temperature gradient analyzed through the thermographic camera at three different instants: t = 0 s; t = 60 s, 
and t = 120 s. 


Figure 13. 
Formation of a bubble at the entrance of the microchannels, which affects the thermal performance of the 
microfluidic device. 


thelow conductivity and partial opacity to the infrared (IR) radiation, were not 
favorable for the experiments. Also, the walls of the heat sink were rough due to the 
ABS master mold used in the fabrication technique, and the thickness should be 
reduced. The stability of the nanofluids has yet to be optimized since the deposition 
of NPs was detected on the heat sink walls. Future works will aim to improve of 
those limitations. 


4. Conclusion 


The main objective of this work was to show the potential of a PDMS heat 
sink microfluidic device to perform flows and heat transfer studies of nanofluids. 
The PDMS heat sink device was produced by using the FDM 3D printing process, 
combined with a PDMS casting technique. This fabrication process allowed to 
manufacture devices in an easy, low-cost, and reasonable reproductively way. To 
demonstrate the potential of the produced PDMS heat sink device, fluid flow and 
heat transfer studies were performed by using two different nanofluids, i.e., alumina 
(Al; O3)- and iron oxide (Fe4 O4)-based nanofluids with concentrations of 1 and 
2.596. Overall, it was found that the thermal performance of the working nanofluids 
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is in good agreement with several past studies. For instance, it was noted that the 
heat energy absorbed by both nanofluids was higher than that absorbed by the 
distilled water. In addition, it was found that the flow rate affects the amount of heat 
absorbed by the nanofluids. However, the most interesting and unique results were 
the optical and thermal imaging results. These results were only possible due to the 
optical transparency of the PDMS heat sink device. Hence, by using this device, 

it was possible to visualize several flow phenomena of the nanofluids such as the 
formation, growing, and breakdown of NPs clusters. From these latter observations, 
it was possible to conclude that one of the main causes for the formation of the 
clusters was the high roughness of the PDMS surface channels caused by the surface 
roughness of the ABS master mold fabricated by the FDM 3D printer. This drawback 
can be overcome by performing an acetone vapor treatment before performing 

the PDMS casting. Overall, the simplicity, low-cost, and unique features of the 
proposed PDMS heat sink microfluidic device may prove a viable alternative tool to 
investigate nanofluids flow and heat transfer phenomena that are not possible to be 
performed by the current traditional systems. 
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Nomenclature 

A area 

cp heat capacity 

h convection heat transfer coefficient 
K thermal conductivity 

N number of microchannels 
q volume flow rate 

Q heat rate 

R resistance 

Rt thermic resistance 

T temperature 

wp average width of fin 

ZC height of the channel 
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Greek symbols 

"n fin efficiency 

p viscosity 

p nanoparticle concentration 
p density 

Subscript 

b microchannel bottom 
bf base fluid 

ext. exterior 

f fluid 

in inlet 

l microchannel sidewall 
nf nanofluid 

out outlet 

p nanoparticles 
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Chapter 9 


The Characteristics and 
Application of Nanofluids in MQL 


and MQCL for Sustainable Cutting 
Processes 


Tran The Long and Tran Minh Duc 


Abstract 


Recently, there has been growing attention to nanofluids, especially in industry. 
More and more people nowadays see nanoparticle applications in various fields such 
as automotive, agriculture, medicine, machining, and so on. The addition of differ- 
ent nanoparticles to fluids has shown enormous advantages, particularly for improv- 
ing the efficiency and therefore lowering the energy consumption of processes for 
addressing a wide range of global challenges related with energy and environmental 
problems. Nanoparticles are of great scientific interest as they are in that nanofluid 
with unusual effects, and ultra-small sizes will be a new area for researchers and 
definitely offer novel mechanisms and technologies in the future. In this chapter, the 
authors will mainly present the characteristics as well as latest advances in applica- 
tions of nanofluids in machining practices. Nanoparticle additives contribute to 
reduce friction coefficient, lower the energy consumption, and significantly extend 
tool life by lowering thermal stress, from which the surface quality of manufactured 
parts improves. Moreover, the nanoparticle application in some of the green tech- 
nologies as MQL and MQCL using vegetable oils not only brings out superior cooling 
and lubricating properties and minimizes the use of cutting fluids, but also creates 
new solutions for machining, especially for difficult-to-cut materials. 


Keywords: nanoparticles, nanofluid, sustainable cutting, MQL, MQCL, 
hard machining, vegetable oil, metal cutting 


1. Introduction 


Climate change has become the most growing concern of people around the 
world. The rapid increase of population consequently leads to the use of more 
natural resources and giving out of more waste. The pollution in air, water, and food 
causes many serious human diseases. Accordingly, environmental laws are continu- 
ously tightening up to protect our Earth. Being a part of the production chain, 
manufacturing engineers are demanded not only to produce the products to meet 
the growing demand for higher quality and productivity but also to be responsible 
for achieving the sustainability in manufacturing. In metal cutting industries, the 
used cutting fluids after using for cooling and lubricating the contact zone con- 
tribute the largest amount of disposal (around 30%), which finally ends up as the 
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contamination in the rivers leading to the water pollution [1]. Therefore, it is neces- 
sary to find the solutions to reduce or eliminate the usage of coolants. Over some last 
decades, there were numerous studies concerning the reduction of coolant usage in 
machining, and dry cutting processes, the truly environmental-friendly method, 
had drawn most attention and brought out the obvious cost benefits derived from 
the elimination and treatment of cutting fluids. However, the selection of the proper 
cutting tools or inserts plays a very important role to ensure the proper tool life and 
high precision and accuracy of machined parts [2], and it also causes a strong influ- 
ence on technological and economic characteristics. Recently, to meet the continu- 
ously increasing demand for cutting difficult-to-cut materials having high-graded 
mechanical properties and high hardness, the tools with geometrically defined cut- 
ting edges are directly used for machining the heat-treated materials, with the typi- 
cal hardness of 45-70 HRC [3]. These processes are called hard machining, which 
has become the research trend in mechanical applications due to high productivity 
and accuracy. Up to now, people have seen hard machining more in metal cutting 
field, and therefore many of traditional grinding processes have been replaced. The 
new approach not only provides the alternative solution for cutting hard materials 
but also improves the cutting performance, significantly reduces coolant usage, and 
has low machine tool investment. On the other hand, the thermal shock caused by 
the use of cutting fluids must be seriously considered to avoid the insert breakage, 
so the flood cooling is not usually used for hard machining processes, especially for 
interrupted cutting. Furthermore, the enormous heat and high forces arising from 
cutting zone are the most challenging problems of hard cutting processes, which 
always demand the appropriate uses of high-graded cutting tools like coated carbide, 
ceramic, polycrystalline cubic boron nitride (PCBN), and diamond inserts [2, 4, 5]. 
Accordingly, minimum quantity lubrication (MQL) technique was proposed and 
proven to use and exhibited the promising results in some last decades [1, 6, 7]. 

The cutting fluids in forms of oil mist are directly sprayed to cutting zone, so the 
lubricating effect is very high to decrease the friction coefficient, from which 
cutting forces, cutting temperature, and tool wear reduce significantly, and tool life 
is extended. Interestingly, the minimal use of cutting fluid makes MQL an environ- 
mental friendly technique, and the vegetable oils can be used for hard cutting, which 
contributes to protect environment [8]. The main drawback of MQL method is the 
low cooling effect, which limits the applicability and cutting performance of hard 
machining [9, 10]. In order to develop MQL technique, there have been many studies 
proposing the very promising solutions to enhance the cooling performance, which 
includes MQL using nanofluids, minimum quantity cooling lubrication (MQCL), 
and MQCL using nano additives. In this chapter, the authors mainly discuss the 
latest studies on those up-to-date techniques used in hard machining processes. 

This chapter is divided into five sections. Section 1 of the chapter provides the 
literature review of new development of MQL and MQCL technology using nano- 
fluids for sustainable cutting processes. Section 2 is dedicated to hard machining 
under MQL condition using nanofluid. Section 3 describes the application of MQCL 
condition based on the new approach for hard cutting processes. Section 4 contains 
the latest advances on the utilization of nano additives for improving MQCL hard 
machining performance. Finally, Section 5 draws out the conclusions and some 
suggestion for future work. 


2. Hard machining under MQL condition using nanofluid 


Using nano additives suspended in MQL based fluids has opened a new 
approach for machining difficult-to-cut materials and is also an up-to-date research 
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topic gaining the growing concerns, especially for encountering climate change. 
There are many types of nanoparticles, such as Al,O3, MoS), SiO2, ZrO, CuO, TiO, 
CNT, ND, and so on, proven to use for improving the tribological property, thermal 
conductivity, and viscosity [11]. 


2.1 The improvement of cutting performance 


In order to apply this technique in machining practice, the parameters of MQL 
(the based fluid, air pressure, flow rate) and nanofluid (the type, size, and con- 
centration of nanoparticles) are needed to study and optimize, because they have 
strong effects on the cutting process. If the inappropriate values of each parameter 
are chosen, the little effectiveness and even the negative influence may occur in 
machining responses. 

Li et al. [12] investigated MQL grinding process for Ni-based alloy using six 
different types of nanofluids. The results indicated that the viscosity and thermal 
conductivity of nanofluids significantly improve when compared to the base 
fluids. The authors also pointed out that CNT nano additives exhibit the highest 
heat transfer coefficient. Hence, the cutting temperature and forces decrease. 
Another observation done by Ali et al. [13] indicated that the viscosity of Al,O; and 
TiO; nano-lubricants increases while their kinematic viscosity decreases slightly. 
Through the experiments, the coefficient of friction, power consumption, and wear 
rate much reduced due to the rolling performance together with the formation of 
tribo-films created by Al,O3 nanoparticles. Moreover, the technical specification 
and concentration of nanoparticles play a very important role and strongly influ- 
ence on the machining responses. For finish cutting, the nanoparticles with smaller 
grain size and higher concentration should be used to improve the surface quality 
and reduce the cutting forces [11]. Pashmforoush and his co-authors [14] reported 
that the big improvement in surface roughness in grinding process of Inconel 738 
super alloy is about 62.16 and 36.3696 compared to those of dry and flood condi- 
tions, respectively. The enhancement of lubricating effect was also reported in 
milling under MQL using MoS; nanofluid, from which the friction coefficient 
reduced to extend the tool life and improve the surface quality [15]. The presence 
of nanoparticles in MQL based fluid not only improves the cooling and lubricating 
effects but also brings out the better cutting performance for machining difficult- 
to-cut materials. Moreover, this approach will successfully replace the dry and flood 
conditions, which fulfill the technological, economic, and environmental require- 
ments, suitable for modern manufacturing. The performance investigation of end 
milling of SKD 11 steel using HSS tools under nanofluid MQL was done in [16]. The 
experiments are set up and shown in Figure 1. The cutting condition includes the 
three values of cutting speeds of 18, 24, and 30 m/min, feed rate of 0.01 mm/tooth, 
and axial depth of cut of 3 mm. The diameter of end mill is 10 mm. Al,O; nano 
additives (0.5 wt96.) are enriched in emulsion and soybean-based oil. 

Through experimental results, the cutting force components Fy, F,, and F, when 
changing the cutting speeds and based fluid are shown in Figures 2-4. Under NFMQL 
with emulsion-based oil, the cutting forces F,, F, and F, reduce with the increase of 
cutting speed from 18 to 30 m/min. Ata cutting speed of 30 m/min, the comparison 
between emulsion and soybean oil is made to find out the effect of the based fluid on 
machining responses. It can be clearly seen that the cutting forces reduce due to the 
better lubricating performance of soybean oil compared to emulsion fluid. 

In addition, the tool life under NFMQL using soybean oil much improves and 
is over two times longer than that under the case using emulsion fluid (Figure 5). 
The investigation of tool wear is shown in Figures 6-8. They clearly reveal that 
notch wear and flank wear on HSS end mills increase with the rise of cutting 
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Figure 1. 
Experimental setup [16]. 
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Figure2. 
The relation of cutting speeds and nanofluids to the cutting force F, [16]. 


speeds from 18 to 30 m/min under MQL using emulsion. In addition to that, the 
burn marks caused by heat deterioration develop. It can be explained that SKD 11 
tool steel has extremely high wear-resistant properties due to a high carbon and 
chromium (12% chrome) in chemical composition, from which it is grouped in 
difficult-to-cut material. In increasing the cutting speed from 18 to 30 m/min, 
MQL with emulsion-based fluid did not provide sufficient lubricating effects, so 
cutting temperature increased rapidly to damage end mills. In contrast, soybean 
oil has higher viscosity than that of emulsion, and the presence of Al,O3 nanopar- 
ticle additives contributes to improve the cooling and lubricating performance. 
Therefore, soybean oil-based nanofluid exhibits the superior lubricating effects 
to reduce friction coefficient in the cutting zone due to the easier formation of oil 
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Figure 3. 
The relation of cutting speeds and nanofluids to the cutting force F, [16]. 
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Figure 4. 
The relation of cutting speeds and nanofluids to the cutting force F, [16]. 


mist. In this study of end milling process of SKD 11 steel before heat treatment, 
cutting heat did not exceed the ignition temperature of soybean oil. That is the 
main reason why MQL using soybean-based nanofluid is better in this situation. 
Furthermore, AlO; nanoparticles with nearly sphere morphology suspended in oil 
mist as "the rollers" play an important role in improving cooling and lubricating 
effects. From those reasons, notch wear and flank wear on HSS end mills at cutting 
speed V, = 30 m/min significantly reduce. Notch and flank wear lands after 85 min 
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Figure 5. 
The relation of cutting speeds and nano-cutting fluids to the tool life [16]. 


Figure 6. 
Notch wear and flank wear (Vc = 18 m/min, emulsion-based nanofluid) [16]. 


Figure 7. 
Notch wear and flank wear (Vc = 24 m/min, emulsion-based nanofluid) [16]. 


of cutting are even lower than those of the case with MQL using emulsion-based 
nanofluid after 40 min of cutting (Figures 8-9). The significant reduction of burn 
marks indicates that the cooling and lubricating performance of soybean-based 
nanofluid is better and also suits for sustainable production due to the use of 
vegetable oil. Hence, tool life of end mill increases to 85 min even at cutting speed 
of 30 m/min, which is also higher than the manufacturers’ recommendations [16]. 
According to ISO 8688-2:1989 (en) [17], the cutting speed for soft steels using 
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Figure 8. 
Notch wear and flank wear (Vc = 30 m/min, emulsion-based nanofluid) at 40 min [16]. 


Figure 9. 
Notch wear and flank wear (Vc = 30 m/min, soybean-based nanofluid) at 85 min [16]. 


normal HSS end mills is recommended about 30-35 m/min, but for difficult-to-cut 
steels like SKD 11, with hardness 200—250 HB, the cutting speed must be reduced 
to 14-18 m/min to ensure the proper cutting performance and tool life [18]. 
Moreover, the cutting speed also increases from 18 to 30 m/min by using MOL 
emulsion-based nanofluid, which reveals the better cooling and lubricating effects 
compared to the pure fluids. 


22 Theimportant parameter of MQL nanofluid 


The concentration parameter of nanoparticles enriched in based fluids is among 
the most influential on machining outputs and costs, so it had been much studied 
to find out the appropriate and optimized values. Garg et al. [19] investigated 
the concentration effect of nanoparticles on micro-drilling process under MQL 
condition. The experimental results indicated that this parameter caused the 
significant reduction of drilling torque and power consumption. In the study of 
Lee et al. [20], the proper concentration of diamond nanoparticles was found with 
0.05 wt%, from which the reduction of friction coefficient was observed by 23%. 
The authors concluded that diamond nanofluid provided the excellent anti-wear 
and lubricating effects. Zhang et al. [21] studied the concentration parameter of 
MoS; and CNT nano additives in MQL grinding. The improvement in lubricating 
performance contributes to increase surface quality. Furthermore, hybrid MoS,- 
CNT nanofluids provided the superior cooling lubrication compared to that of the 
fluid with a single type of nanoparticles. Luo et al. [22] studied Al,O; nanoparticles 
enriched in MQL based fluid and concluded that Al,O; nanofluid exhibited good 
resistant ability for high temperature. Then, the cutting temperature is not high to 
cause the reduction of wear rate, which is much smaller than that of dry condition. 
Yildirim et al. [23] had done the study of MQL turning process of the difficult-to- 
cut steel Inconel 625 using hBN nano additives. The better lubricating performance 
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and surface roughness are reported from the obtained results, which led to reduce 
friction coefficient and wear rate. The authors also concluded that the optimal 

hBN nanoparticle concentration was 0.5 wt%. The experimental study on Al,O; 
nanoparticle concentration used as MQL based fluid in hard milling had been done 
by using ANOVA analysis and response surface methodology (RSM), from which 
the research direction was made for optimizing the concentration variable [24]. 
Figures 10-11 show the response surface plots of the relation of surface roughness 
and cutting force versus nano concentration (np), cutting speed V,, and feed rate F. 
It can be clearly observed from Figure 10 that, for better surface roughness, the 
low value of nanoparticle concentration about 0.5 wt% is more preferable than the 
larger ones (1.0 and 1.5 wt%). In contrast, the larger concentration (about 1.0 and 
1.5 wt%) contributes to reduce the cutting forces and cutting temperature signifi- 
cantly when compared to the lower one (0.5 wt%.). From those, the wear rate much 
reduces by increasing the concentration of Al,O3 nanoparticles to 1.0-1.5 wt%, so 
the tool life prolongs (Figures 12-15). Accordingly, the nanoparticle concentration 
must be chosen not only to ensure the good tool life but also to maintain the high 
surface quality. 

From Figure 13, the chip colors are well matched with the reduction of cutting 
temperature. The dark purple and blue colors in the case of using np = 0.5 wt% 
change to brown and dark straw one in the case of using np = 1.0—1.5 wt%, which 
indicates that cutting temperature decreases [25]. It can be concluded that the 
decrease of coefficient of friction and generated heat is reported by increasing 
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Figure 10. 
Response surface plots of surface roughness versus nano concentration and cutting speed (a), and nano 
concentration and feed rate (b) [24]. 
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Figure 11. 
Response surface plots of cutting force F; versus nano concentration and cutting speed (a), and nano 
concentration and feed vate (b) [24]. 
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the concentration parameter of Al,O; nanofluid, leading to reduce cutting forces 
and wear rate and prolong the tool life (Figure 14). From Figure 15, it clearly 
reveals that during the first 40 min, the values of surface roughness are higher 
when utilizing the high nanoparticle concentration (1.0-1.5 wt%). After that, the 


Figure 12. 


Wear on flank face under MQL using soybean-based nanofluid at 80 min using different nano concentrations: 
(a) 0.5 wt%, (b) 1.0 wt%, and (c) 1.5 wt% [24]. 
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Figure 13. 
Chip colors and micrographs with different AL,O; nanoparticle concentrations (at 80 min): (a) 0.5 wt%, (b) 
1.0 wt%, and (c) 1.5 wt% [24]. 
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Figure 14. 
Tool life with different Al,O, nanoparticle concentrations [24]. 


rate of reduction of surface roughness values rapidly increases, which is contrary 
to the case of using 0.5 wt%. It is the novel observation obtained from the valida- 
tion experiments, which is conducted until the tool life ends to see the actual 
phenomena after receiving the ANOVA and RSM results. Interestingly, the tool 
life in the case of soybean-based nanofluid 1.5 wt% is equal to that of emulsion- 
based nanofluid 0.5 wt%. It provides an important technical guide to enlarge the 
applicability of vegetable oil-based nanofluid in hard cutting processes as well as 
maintains its environmental-friendly characteristics. 

However, each type of nanoparticles has its own specific morphology and 
property, so the appropriate or optimal concentration parameters are different [26]. 
This is the up-to-date research topic. Accordingly, more investigations are needed 
to make and build up the technical guides for manufacturers in machining practice, 
even though many studies have been done to optimize these variables. In order to 
develop MQL method, minimum quantity cooling lubrication (MQCL) has been 
considered as another promising approach to solve the low cooling performance, 
the main MQL drawback. It is also the newest research topic, which is discussed in 
Sections 3 and 4. 
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Figure 15. 
Surface roughness Ra with different AL,O, nanoparticle concentrations [24]. 


3. Hard machining under MQCL condition 


The large amount of heat generated from hard cutting is always the big chal- 
lenge for selecting cutting tools and cutting condition while ensuring the technical 
requirements, productivity, and proper manufacturing cost. Hence, providing 
appropriate cooling and lubricating effects to cutting zone plays a vital role in 
the development of hard machining processes. The applicability of those can be 
enlarged to some or all of the traditional grinding processes. MQCL technique has 
been proposed and developed to fulfill the cooling and lubricating requirements 
and is also a solution for improving MQL method. Up to now, MQCL has drawn 
much attention and has been studied in recent years. Maruda et al. [27, 28] made the 
study on MQCL parameters using emulsion-based fluid in hard turning process. 
The obtained results indicated that emulsion oil mist formed under MQCL condi- 
tion plays an important role for improving the cooling lubricating performance in 
the cutting zone and increasing the cutting condition. The formation of tribo-films 
tends to occur easily with the droplets with smaller size, which help to decrease 
the coefficient of friction, cutting forces, and wear rate [29, 30]. The better cool- 
ing and lubricating effects of MQCL technique also reflect through the chip shape 
and the reduction of chip thickening coefficient [31]. Pervaiz and his co-authors 
[32] studied MQCL performance in the turning process of difficult-to-cut mate- 
rial Ti6Al4V. The author concluded that cutting forces and tool wear reduced and 
surface quality improved when compared to dry and flood conditions. It reveals the 
better cooling and lubricating effects of MQCL technique. Krolczyk together with 
his co-authors [33] investigated the parametric and nonparametric description of 
the surface topography under dry and MQCL conditions using emulsion-base fluid. 
The study results showed that the nozzle distance causes the strongest influence on 
droplet diameter. The most outstanding finding of this research is that parameters 
can be chosen for oil mist formation in a certain time, which is enough for creating 
cooling and lubricating effects and then evaporating due to generated heat from 
the cutting zone. Based on a brief review, it can be clearly seen that there is little 
information of MQCL technique and most of the studies relied on the based fluid 
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having cooling effect like emulsion oil to form MQCL method. The use of a real 
cooling method assisted to MQL technique to form MQCL condition is a novel 
approach. In this section, the author presents the newest advances in using the 
principle of Ranque-Hilsch vortex tube for separating a compressed gas into hot and 
cold streams from ordinary air [34], in that the cold stream is used to create cooling 
effects combined with MQL method to form MQCL [35]. The deep study on hard 
milling of SKD 11 steel (52-60 HRC) in terms of surface quality under MQCL con- 
dition was done, and the results were compared to dry and MQL conditions. From 
Figure 16, hard milling under MQCL method brought out better surface roughness 
than those under dry and MQL conditions. The main reason is that MQCL tech- 
nique provides sufficient cooling and lubricating effects, especially cooling effect, 
which helps to reduce the cutting temperature and tool wear. 

KEYENCE VHX-6000 Digital Microscope (Keyence Corporation, Osaka, Japan) 
was utilized for studying surface microstructure and surface profile (Figure 17). The 


Dry 
MQL 
MQCL 
= MQCL with MoS2 0.2wt% 
mpe MQCL with MoS2 0.5wt% 
MQCL with MoS2 0.8wt% 
0.202 


0.00 


Figure 16. 
The average values of surface roughness Ra under different cooling and lubricating conditions (cutting speed 
Vc = 110 m/min, feed rate F = 0.012 mm/tooth, depth of cut d = 0.12 mm, hardness of 56 HRC) [35]. 


Figure 17. 
KEYENCE VHX-6000 digital microscope for studying machined surface topography. 
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machined surfaces under different cooling and lubricating conditions are investi- 
gated (Figures 18-22). The white layer and burn marks significantly reduced under 
MQL and MQCL conditions compared to dry cutting because of cooling and lubri- 
cating enhancement. The burn marks under MQCL condition are less than those 
under MQL method due to better cooling performance (Figures 19(a), 20(a)). In 
addition, compared to dry and MQL conditions, the compression of machined sur- 
face observed from the surface profile much reduces (Figures 18(b), 19(b), 20(b)). 


(b) 


Figure 18. 
Surface microstructure (a) and profile (b) under dry condition [35]. 


(b) 


Figure 19. 
Surface microstructure (a) and profile (b) under MQL condition [35]. 


(b) 


Figure 20. 
Surface microstructure (a) and profile (b) under MQCL condition using pure emulsion-based fluid [35]. 
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(b) 


Figure 21. 
Surface microstructure (a) and profile (b) under MQCL condition using emulsion-based nanofluid of MoS, 
0.2 wt% [35]. 


(b) 


Figure 22. 
Surface microstructure (a) and profile (b) under MQCL condition using emulsion-based nanofluid of MoS, 


0.5 wt% [35]. 


4. MQCL hard machining using nanofluids 


The nanofluids used for MQL are successfully proven to be an alternative 
solution for difficult-to-cut materials while maintaining its environmental 
friendly property. Based on this idea, the use of nano additives in MQCL 
method will bring out the promising solution to increase the hard machining 
performance. The study of surface quality under MQCL using MoS, nanofluid 
for hard milling is the first attempt to investigate the cooling and lubricating 
effects [35]. From the obtained results, the values of surface roughness R, under 
MQCL using nanofluids are lower than those of dry and MQL conditions. By 
using MoS, nanoparticle concentration of 0.2 and 0.5 wt%, surface roughness is 
even slightly better than that of MQCL with pure fluid, but the R, value rapidly 
increases when increasing the concentration to 0.8 wt%. It can be explained that 
the morphology of MoS; nanoparticles is ellipsoidal with the low friction coef- 
ficient up to 0.03-0.05 or even lower [36], by which the better lubricating effect 
contributes to improve the surface quality. On the other hand, nanoparticles pos- 
sess the large surface area, which remain on the machined surface to form a thin 
protective film, which amplifies with the increase of MoS, nanoparticle concen- 
tration [37]. Furthermore, it also contributes to form MoS, tribo-film easily [29], 
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which can be observed from the so-called microbubbles on the machined surface 
(Figures 21-22). The protective film reduces and disappears when increasing the 
concentration to 0.8%, which causes the negative effect on surface quality [37]. 
Moreover, the white layer and burn marks are much reduced due to superior 
cooling and lubricating performance under MQCL condition using nanofluid. 
From those, the hard machining ability of normal carbide tools improves sig- 
nificantly and is about 157% higher than manufacturer’s recommendations [38, 
39]. It is the most outstanding finding of this research, and also the proper MoS, 
nanoparticle concentration in emulsion-based fluid was reported about 0.2 and 
0.5 wt%, which provides a very important technical guide for further researches 
and manufacturers. More investigations are necessary to be made for building up 
technical guidelines and optimizing nanofluid parameters. 


5. Conclusion 


The application of nanofluids continues to receive growing attention in basic 
science and machining technology. As shown, nano additives in based fluid of 
MQL and MQCL methods improve the cooling and lubricating effects as well 
as tribological property, thus increasing the cutting performance, especially for 
difficult-to-cut materials. It brings out the alternative solutions for improving 
productivity and reducing manufacturing cost. From those, the applicability MQL 
technique having environmental friendly characteristic has been enlarged in hard 
machining. Furthermore, the use of different types of vegetable oils can fulfill 
the cooling and lubricating performance by suspending nanoparticles, which is 
an interesting research topic and exhibits very promising results. On the other 
hand, MQCL has been considered as another approach for MQL development to 
overcome the low cooling effect. In this chapter, Ranque-Hilsch vortex tube, a real 
cooling method, used for creating cooling effect from ordinary compressed air 
rather than other gas sources to form MQCL method, is the first attempt applied 
to hard cutting processes. Also, nano additives enriched in MQCL-based fluids 
are the latest advances in the field of studying MQL and MQCL techniques. The 
parameters of nanofluid, such as types and size of nanoparticles, concentra- 
tion, and based fluid, play a key role in successful applications in metal cutting 
practice, and more studies are needed to make further development and optimize 
those variables. Those superior cooling and lubricating methods presented in 
this chapter will contribute to the solutions to reduce/eliminate the cutting fluids 
and replace dry and wet conditions. It is suitable for protecting our environment 
and aims for a sustainable production. In the future work, more attention will be 
paid on other types of nanoparticles, concentration, and parameters of MQL and 
MQCL methods. 
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CNTs 
MQL 
MQCL 
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ND 

NF 
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NP 

NPs 
ANOVA 


Author details 


depth of cut (mm) 

cutting speed (m/min) 

feed rate (mm/tooth) 

cutting force (N) 

high speed steel 

hexagonal boron nitride 

carbon nanotubes 

minimum quantity lubrication 
minimum quantity cooling lubrication 
nanofluid minimum quantity lubrication 
nanodiamond 

nanofluid 

nanofluids 

nanoparticle 

nanoparticles 

analysis of variance 
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